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Introduction 






hysics as a general discipline has no Hmils, from | he very huge (galaxy* 

* wide) to the? very small (atoms and smaller). This hook Ls a hoot the very 
small side oE — tEiat's the specialty oE quantum physics. When yuii lyuon* 

tize suincthing, yuu can't go any smaller; you're dealing with discrete units. 

Classical physics is terrific at explaining things like heating cups of cudiec ur 
accelerating down ramps or cars colliding, as well as a million -other things, 
but It has problems when things get very small. Quantum physics usually 
deals with the micro world, such as what happens when you look at Individual 
electrons zipping around. For example, electrons can exhibit both particle and 
wave-dike properties, much to the consternation of experimenters — and It 
took quantum physics to figure out the full picture. 

Quantum physics also kirroduccd the uncertainty principle, which says you 
can't know a particle's exacl position and mumt-mtum at the same time. And 
the Held explains the way that the energy levels o( I he dec trails bound in 
atoms work. Figuring out those ideas all look quantum physics, as physicists 
probed ever deeper lor a way to model reality. Those topics are all coining 
up in this book. 



About This Book 

Because uncertainly and probability are so important In quantum physics, 
you can't fully appreciate the subject without getting into calculus. This book 
presents the need-to-know concepts, but you don't see much in the way of 
though! experiments liiaL deal with eats or parallel universes, I focus on the 
mailt and how it describes the quantum world. 

I've tauRht physics to many thousands of students at the university level, 
and from that experience, l know most of them share one common trait: 
Confusion as to what they did to deserve such torture. 

Quantum Physics far Dtimmiea largely maps to a college course, but this book 
is different from standard texts. Instead nl writing il from the physicist's or 
professor's point of view. I've tried lo write it from I he reader's point of view. 

In other words. I’ve designed this book to be crammed lull of the good stuff — 
and only the good stuff. Nol only that, but you can discover ways of looking al 
things that professors and teachers use to make figuring out problems simple. 
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Although 1 encourage you to read this book from start to finish, you can also 
leaf through this hook as you like, reading the topics that you Bind Interest- 
ing. hike other for fiummm hooks, this one lets you skip around as you Like 
as much as possible. You don't have to read the chapters in order If you 
don't want to. This is your hook, and quantum physics Is your oyster, 



Contentions Used in This Book 

Some books have a dozen dizzying conventions that you need to know before 
you can even start. Not this one, Here’s all you need to know; 

v* I put new terms in italics, like tffrs. the first rime they’re discussed; 
it Eullow them with a definition. 

v* Vectors — those items that have hoth a magnitude and a direction — 
arc. 1 given ill bold, Like this: 11. 

t* Web addresses appear in mcneEont. 



Foolish Assumptions 

I don't assume that you have any knowledge of quantum physics when you 
start To read this hook. However, I do make the ini lowing assumptions: 

V You’re taking a college course in quantum physics,, or you're interested 
in how math describes motion and energy on the atomic and subatomic 
scale . 

** You have some math prowess. In particular, you know some calculus. 
You don't need to he a math pro. hut you should know how to perform 
Integration and deal with differential equations. Ideally, you also have 
some experience with HLI b err space- 

v* You have some physics background as well. You've hart a year’s worth 
of culJege-level physics (or understand ail that's in Physics For Dummies) 
before you tackle this one. 



Haul This Book Is Organized 

Quantum physics — the study of very small objects — to actually a very big 
topic. Tn handle it, quantum physicists break the world down intn different: 
parts. Here are the various parts thafr are coming up in this hook. 
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Part 1: Small World, Huh? Essential 
Quantum Physics 

Part l is where you start your quantum physics journey, anti you get a goad 
overview of I he topic here. L survey quantum physics and telJ yuu what It's 
good for and what kinds- td problems it can solve. You also get a good foun- 
dation In l he math tJiat you need for I he rest of tile bunk, such as state vec- 
tors and quantum matrix manipulations. Knowing this stulf prepares you to 
handle the other parts. 

Part 11: Bound and Undetermined: 
Handling Particles in Bound States 

Particles can he trapped inside potentials; for instance, electrons can be 
bound in an atom. Quantum physics excels at predicting the energy levels ui 
particles bound in various potentials, and that's what Part II covers. You see 
how to handle particles hound in square wells and in harmonic oscillators. 

Part 111: Turning to Angular 
Momentum and Spin 

Quantum physics lets you work with the micro world in terms i>f the angu- 
Lar irinmerilum of particles, as well as the spin of electrons. Many tarn mis 
experiments — Such as the Stcrn-GeriaCh experiment, in which beams ui par- 
tides split in magnetic fields — are understandable only in terms of quantum 
physics, and you gel all the details here. 

Part IV: Multiple Dimensions: Going 3D 
With Quantum Physics 

In |he first three parts, all the quantum physics problems are iim-(li:iii'iiHi:i:i.d 
to make life a little easier while you're understanding how to solve those 
problems. In Part IV, you branch out lo working with three-dimensional prob- 
lems in both rectangular and spherical coordinate systems. Taking things 
(rom ID to 3D gives you a belter picture of what happens in the real world. 
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Part V: Group Dynamics: Introducing 
Multiple Particles 

In this. part, you work wit]] iEUiltfple^parlirle systems, sucfi as atoms and 
gases. You sec huw to tiaiadJe many electrons in atoms, particles interacting 
with other particles, anti particles that scatter otl other particles. 

Dealing Willi multiple particles is all another step in modeling reality — alter 
all, systems with only a single particle don’t take you very far in the real 
world, which Is built of mega, mega systems of particles. In Pari V, you see 
how quantum physics can handle the situation. 



Part Vh The Pan of Tens 

You see the Part of the Tens in all For Dummies hooks. This part is made 
up of last-paced lists of ten Items each. You get to see some of the ton best 
on line tutorials nil quantum physics and a discussion ol quantum physics' 
ten greatest triumphs. 



Icons Used in This Book 

You find a handful of icons in this hook, and here's what they mean: 

This iron flags particularly good advice, especially when you're solving 
problems. 



This icon marks something to remember, such as a law of physics or a particu- 
larly juicy equation. 



This icon means that what follows is technical, insider stuff. You don't have to 
read if if you rirrn'l want to, hut if you wans to become a quantum physics pro 
[and who doesn't?), take a look. 



This iron helps you avoid mathematical nr conceptual slip-ups „ 
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Where to Go from Here 

AIJ right, you're all set and ready to go. You i:aia furn |> in anywhere you 
like. For instance. ii you'tt lure electron spin is going n> be a big topic ol 
conversation a? a party lids weekend, check out Chapter 6. And ii your 
upcoming vacation to Geneva, Switzerland, includes a side trip to your new 
favorite particle accelerator — the Large Hadron Collider — you cart flip to 
Chapter 12 and read up on scattering theory. But If you want to get the full 
Story From the beginning. Jump into Chapter 1 first — that’s where the action 
starts. 
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Parti 

Small World, Huh? 
Essential Quantum 
Physics 



The 5 th Wave By Rich Tennant 




just like the regular stevr oidy ti > s got 
Some bits o£ m&tier iri it vre car^t identity/ 
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In this part . . . 



T hus pari is designed to give you an introduction to the 
ways of {[ixanlmn physics. Yon see the issues that 
gave rise to quantum physics anti the kinds ul solutions it 
provide!;. II also introduce you lo the kind ol math that 
quantum physics requires, including the notion of stale 
vectors. 
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Chapter 1 

Discoveries and Essential 
Quantum Physics 



In Thi$ Chapter 

Pulling lorlh theories of quantization and discrete units 
tj{ per iriRnring with waves acting as particles 
Experimenting with particles acting as waves 
Embracing uncertainty and probability 



AM ccording to classical physics, particles are particles and waves are 
¥ ¥ waves, and never the twain shall mix That Is, panicles have an energy 
E. and a momentum vector p, and that's the end of It . And waves, such as light 
waves, have an amplitude A and a wave vector k ( where the magnitude of k - 

y. where }. Is the wavelength) that points in the direction the wave Is trav- 
eling. And that's the cud ul that, too, according tu classical physics. 

But the reality 3s dfclf ureal — particles turn nut tu exhibit wave-like proper- 
ties, and waves exhibit particle like properties as well, The Idea that waves 
(like light) eon act as particles (like electrons) and vice versa was the major 
revelation that ushered In quantum physics os such an important part of the 
world of physics. This chapter takes a look at the challenges facing classical 
physics around the turn of the 20th century — and how quantum physics 
gradually came to the rescue, Up to that point, the classical way ol looking 
at physios was thought to expLaiu just about everything. Out as those pesky 
experimental physicists have a way ol do mg, they cdinc up with a bunch ol 
experiments that the theoretical physicists couldn't explain. 



That made the theoretical physicists mad. and they got on the tub. The prob- 
lem here was the microscopic world — the world that's too tiny to see. On 
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the larger scale, classical physics could still explain most of what was going 
on — but when it came to effects that depended on the micro-world, classi- 
cal physics began to break down. Taking a look at how classical physics col- 
lapsed gives you an introduction to quantum physics that shows why people 
needed It, 

Being discrete: The Trouble 
ufith Black-Body Radiation 

One oi the major ideas ol quantum physics is. well, quantization — measuring 
quantities in discrete, not continuous, units. The idea of quantized energies 
arose with one ol the earliest challenges to classical physics; the problem of 
black-body radiation. 

When you heat an object, it begins to glow. Uven before the glow is visible, 
it’s radiating in the infrared spectrum. The reason it glows is that os you heal 
it. the electrons on the surface ol Use mate rial are agitated thermally, and 
electrons being accelerated and decelerated radiate light. 

Physics in the Late 1 9 Lb and early 20th centuries was concerned with Lhe 
spectrum of Light being emitted by black bodies. A black body is a piece of 
material that radiates corresponding to its temperature — but it also absorbs 
and reflects light from its surroundings. To make matters easier, physics pos- 
tulated a black body that reflected nothing and absorbed alt the- light Falling 
on it (hence the term black foody because the object would appear perfectly 
black as it absorbed all light falling on it). When you heat a Mark body, it 
would radiate, emitting Sigbl. 

Well, it was hard to come up with a physical black body — after all. what 
material absorbs light 100 percent and doesn’t reflect anything? But the 
physicists were clever about this, and they came up with lhe hollow cavity 
you see in Figure l-l, with a hale in it. 

When you shine light on the hole, all that light would go Inside, where It 
would be reflected Again and again — until it got absorbed (a negligible 
amount of Sight would escaper through the hole). And when you heated the 
hollow cavity, the hole would begin to glow. Sn there you have is — a pretty 
good apprnetimaTinn ol a blank body. 
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You can see the spectrum of a black body [and attempts to model that spec- 
trum) in Figure 3-2, for two different temperatures, T. and T . The problem 
was that nobody was able to come up with a theoretical explanation for the 
spectrum of light generated by the black body. Everything classical physics 
could come up with went wrong.. 



Figure l-i: 

SlaC t:- bLidy 
radietion 
spe-Mrun). 
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First attempt; Wien's Formula 

The first one lo try to explain the spectrum of a black body was Wlllhelm 
Wien, In 1889. Using classical thermodynamics, he came op with this formula; 

tr(u.r)= AuV ' T ' 1 



where A and |1 arc constants you determine from your physical setup, u is Hie 
frequency oi the light. aitd T is the temperature of the black body. (The spec- 
trum is given by w[n, T], which is the energy density of the emitted light as a 
function ol frequency and temperature.) 

This equation. Wien’s formula, worked fine for high frequencies, as you can 
see in Figure 1-2; however, .i failed for low frequencies. 



Second attempt: Raleigh -Jeans LaW 

Next up in the attempt to explain the black-body spectrum was the HaleiRh- 
Jcaus Law. introduced around IfltHU. Tins law predicted that the spectrum of 
a black body was 

dVr)=§£R' tr 

■ c ' 



where h is Bulbiiann's constant (approximately ] .381)7 x |ii _3;; J K' 1 ). However, 
the Raleigh- Jeans Law hadiheoppositeptobLcmoi Wien's law: Although it 
worked well at Low frequencies (see Figure 8-2). if didn’t match the higher- 
frequency data at all — in fad. if diverged a! higher frequencies. This, was 
called the uti/aviotei catastrophe because the best predictions available 
diverged at: high frequencies (corresponding to ultraviolet light), it was time 
For quantum physics to take over. 

An intuitive (quantum ) leap: 

Max Planck's spectrum 

The black-body problem was a tough one to solve, and with it came the first begin- 
nings of quantum physics. Max Planck came up with a radical suggestion — 
what if She amount of energy that a light wave can exchange with matter 
wasn't continuous, as postulated by classical physics, butddscrefe^ln other 
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words, Planck postulated that the energy of the light emitted from the wall? 
of the h lack-body cavity came only in integer multiples like this, where h is a 
universal constant: 



E = jnhti, where jt = (), 1, 2 , 



With this theory, crazy as it sounded in the early ISOOs, Planck converted the 
continuous integrals used by Raleigh-Jeans to discrete sums over an infinite 
number of terms. Making that simple change gave Planck the following equa- 
tion for the apectmm of black-body radiation; 



tf(u,T}= 




( U".- _ i 



This equation got It right — it exactly describes the black-body spec tmm, 
both at low and high (and medium, tor that matter) frequencies. 

This idea was quite new. What Planck was saying was that tire energy of 
the radiating oscillators in the black body cooltlnT take uji just any level oi 
energy, as classical physics allows.: it could take on only specific, tfiumiizetl 
energies. In fact. Planck hypothesized that that was true lor any oscillator — 
that its energy was an integral multiple of hx<„ 

And so Planck’s equation came lo be known as Planck s quantization niie, and 
h became Pianck's constant: ft = &.ft£6 * 1 1> 14 Joule-seconds. Saying that the 
energy' or all oscillators was quantized was the birth oF quantum physics. 

Dtp has to wonder how Planck came up with his theory, because it's not 
an nbvinns hypothesis. Oscillators ran oscillate only at discrete energies? 
Where did that come from? In any case, the revolution was an — and there 
was no stopping it. 

The First Pieces: Seeing Light 
as Particles 

Light as particles? Isn't light marie up of waves? Light, it turns out, exhibits 
properties ai troth waves and particles. This section shows you some of the 
evidence. 
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Solving the photoelectric effect 

The photoelectric erFect was one of many ok per Imental results that made up 
a crisis for classical physics around the turn oF the 20th century, it was also 
one of Flnstein’s first successes, and It provides prool of Hie quanlUatlon of 
light. HerE H s what happened. 

When you shire light onto metal, as figure 1-3 shows, you get emitted 
electrons. The electrons absorb the light you shine, Slid if they o; l- 3 encmgli 
utterly, they're able to break I roe oi the metal's surface. According to clas- 
sical physics, liyjii is just a wave, and it can exchange any amount of energy 
with the metal, When you beam light on a piece of metal, the electrons in the 
metal should absorb the light and slowly gel up enough energy to be emit- 
ted from the metal. The Idea was that If you were to shine more light onto 
the metal, the electrons should he emitted with a higher kinetic energy. Arid 
very weak light shouldn’t be able to emit electrons at all, except in a matter 
of hours. 

But that's nut what happened — electrons were emitted as soon as someone 
thoue light oil the metal. In fact, no matter how weak the intensity of t he 
incident light (and researchers tried experiments with such weak light that it 
should have taken hours to get any electrons emitted), electrons luere emit- 
ted. Immediately. 



Light 




rial 



Experiments with the photoelectric eFFect showed that the kinetic energy, K, oF 
the emitted electrons depended only on the frequency — not the intensity — 
Of the incident light, as you can see in figure 1-4. 
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Figure 1-4: 
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In Figure M, x\, Is called the thteshoid ftctfueiKv, and if you shine light with a (re* 
quency below this threshold on the metal, no electrons are emitted. The emitted 
electrons come from the pool of free electrons In the metal (all metals have a poo! 
of free electrons), and you need to supply these electrons with an energy equiva- 
lent to the metal’s work function, W, to emif the electron from the metal’s surface. 

■J'Jic results were hard to explain, classically, so enter Einstein. This was the 
beginning ol lib heyday, around I9(i5. Encouraged by 1'Janck‘s success (sec 
the preceding section), Einstein pustulated tiiat not only were uscillaturs 
quaiith ed but so was light — into discrete units called pltaiom. Light. lie sug- 
gested. acted like particles as well as waves. 

So In this scheme, when light hits a metal surface, photons hit the free elec- 
trons, and an electroa completely absorbs each photon. When the energy, 
dn, of the photon Is greater than the work Function of the metal, the electron 
Is emitted- That is, 

Ad = W * K 

w lurre W is the in ulars work Junction and K is the kinetic energy of the emit- 
ted alec iron. Solving for K gives you the following: 

K=Au-W 

Yon can also write this In terms ol the threshold Frequency this way; 
K=A<u-w„) 
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So apparently, Ll.|*ht Isn't Just a wave: yon can also view It as a particle, line 
photon, In other words, light is quantised.. 

That was also quite m unexpected piece of work fry Ffnstein. although it was 
based obi the earlier work of Planck, light Quantized? Light coming in discrete 
energy packets'.' What next? 

Scattering tight off electrons: 

The Compton effect 

To a world that still had trouble comprehending light as particles (see- 
the preceding section), Arthur Compton supplied the Final blow with the 
Compton effect. His experiment involved scattering photons off electrons, as 
Figure 1-& shows. 



Figure 1 -5: 

Light 

ncldentun 

an electron Photon 
n; rest \ 







El&ctron at rast 



incident light comes in with a. wavelength of f. and hits the electron at rest, 
After that happens, the light is scattered, as you see tn Figure 14. 
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Classically* here's what should’ve happened: The electron should've absorbed 
the Incident light, oscillated. and emitted It— with the same wavelength but 
with an intensity depending on the intensity of the Incident light. But that’s 
not what happened — In tact, the wavelength of the light is actually changed 
by a?-, called the mwkwgih shift The scattered light has a wavelength of ?- 
t aa — in other words* Its wavelength has increased, which means the light 
has lust energy. And Aa depends oh the scattering angle. 6, not on the inten- 
sity ul the incident ligtit. 

Arthur Compton coukl explain the results, ol Ills experiment only by making 
rite assumption Lhai he was actually dealing with two particles — a photon and 
an electron. That is. lie treated ligliL as a discrete parLide. not a wave. And. he 
made the assumption that the photon and the electron collided elastically — 
that is, that both total energy and momentum were conserved. 

Making the assumption that both the light and the electron were particles, 
Compton then derived (his formula for the wavelength shift (it’s an easy cal- 
culation if you assume That the litjhs is represented by a photon with energy 
E = hu and that its momentum is p = %): 

■ —— i ] - ccs n ) 

mx' ■ 

where h is Planck’s constant. m„ is ihe mass of an dec iron, c is the speed ol 
light, and e is the scattering angle of the light. 

You also see this equation in the equivalent form: 

Aa = 4/1.1. sin "| ^ j 

where a., Is the Compton wavelength ol an electron. \ - where 
k = hfln. And experiment confirms this relation — both equations* 



Note that Jo derive the wavelength shift, Compton had to make the assump- 
tion that here, light was acting as a particle, not as a wave. That is p the par- 
ticle nature of light was the aspect of the light that was predominant. 



Proof positron ? Dirac and pair production 

In 1928. the physicist Paul Dirac posited the existence of a positively charged 
anti-electron, the pasd/on. He did this by taking the newly evolving field of 
quantum physics to new territory by combining relativity with quantum 
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mechanics to create relativistic quantum mechanics —and that was the 
theory that predicted, through a pkis/mlnns-sign interchange — the exis- 
tence of the positron. 

It was a hold prediction — an anti-particie of the electron? But just four years- 
later, physicists actually saw the positron. Today's high-powered elementary 
particle physics has alt lands ol synchrotrons arid other particle accelerators 
tu create all l lie elementary particles they need, but in the early 2lJth century, 
this wasn't always Sir. 

hi those days, physicists relied on cosmic rays — those particles and high- 
powered photons (called gamma rays) that strike the Earth firom outer space -= 
as their source of particles. They used ckmU'fmmiwn, which were filled 
with vapor from dry ice, to see the trails such particles left. They put their 
chambers into magnetic fiefd-s to be able to measure the momentum of the 
particles as they curved in those fields, 

In 1(132, a physicist noticed a surprising event. A pair of particles., oppositely 
charged (which could Im determined from the way they Curved in Site mag- 
netic field) appeared from apparently nowhere. Mo panicle trail led to the 
origin of Hie two particles (hat appeared. That was pair-production — the con- 
version ol a high-powered photon LnLo an electron and posiLr-on, which can 
happen when the photon passes near a heavy atomic nucleus. 

So experimentally, physicists bad now seen a photon turning into a pair of 
particles. Wow. As if everyone needed more evidence of the particle nature 

of light, l.iitfir oil, researchers also saw pnfr (innjJiiiaUnn: the conversion of an 
Electron arid positron into pure light. 

Pair production and annihilation turned out to be governed tty Einstein's 
newly introduced theory ol relativity — in particular, his most famous for- 
mula, E = rue , which gives the pure energy equivalent of mass. At this point, 
there was an abundance ol evidence of the particle like aspects -of light. 

A Dual Identity: Lookiny at 
Particles as Waves 

In 1923. the physicist Louis de Broglie suggested that not only did waves 
exhibit part idol ike as pects but the reverse was also true — all material par- 
ticles should display wave-like properties. 
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How dots tills work? For a photon, momentum p = where t- is the 

photon's Irotiucitcy arid A is its wavelength. And the wave vector, k, is equal 
to ft = pifi. where A = ftffa. Dc Broglie said that tlie same relation should. hold 
lor all material particles. That is. 




P 






De Broglie presented these apparently surprising suggestions In his Phi), 
thesis. Researchers pul these suggestions to the test by sending a beam 
through a dual-silt apparatus to see whether the electron beam would act like 
it was made up of particles or waves, in Figure 1-7. you can see the setup and. 
the results. 




twn slits. 

a 




b 




c 



In Figure l-7a. you can see a beam oi electrons passing through a single slit 
and the resulting pattern on a screen. In Figure I -7b. the electrons are pass- 
ing through a second slit. Classically* you'd expect the intensities of Figure 
l-7a mid Mb simply to add when both silts are open; 

I - I, * I. 

But that's not what happened. What actually appeared was an interference 
pattern when both slits were open (Figure l-7e j, not just a sum nl the two 
slits" electron intensities. 
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The result was a validation of de Broglie's invention of matter waves. 
Experiment tore out the relation that fc - and de Broglie was a success. 

The idea of matter waves is a big part of what's coming up in the rest of the 
hook, in particular, the existence ol matter waves says that you add the 
waves' amplitude, ^,0; t) and xp .(>: 0, nor their Intensities, to sum them; 

¥(/, !") = V,0. !} - V,f/- .'j 

You square the amplitude to get the intensity, and the phase difference 
beiween w,<r. Qand v-.O- 0 is wliat actually creates the interference pattern 
that's observed. 

J hu Can't Know Everything (But lfou 
Can Figure the Odds) 

So particles apparently exhibit wave-like properties, and waves exhibit 
particle-like properties. But if you have an electron, which is it — a wave or a 
particle? The truth is that physically, an electron is gust an electron, and you 
can't actually say whether it's a wave or a particle. The act of mensurertrerrf is 
what brings out the wave or particle properties. You see more about this idea 
throughout the book. 

Quantum mechanics lives with an uncertain picture quite happily. That view 
nffended many eminent physicists of the time — notably Albert Einstein, who 
said, famously, “find does not play dice." In this section, i discuss the idea of 
uncertainty and how quantum physicists wort in probabilities instead. 





The Heisenberg uncertainty principle 

The (act that mailer exhibits wave-like properties gives rise to more trouble — 
waves aren't localized in space. And knowing that inspired Werner Heisenberg, 
In 1927. to come up with Ills celebrated uncertainty principle. 

You can completely describe objects In classical physics by their momentum 
and position, both of which you can monsure exactly, in other words, classi- 
cal physics is completely deterministic. 

On the atomic level, however, quantum physics paints a different picture. 
Here, tiie Heisenberg uncertainty principle says rfiat there’s an inherent uncer- 
tainty in the relation between position and momentum. In the jr direction , for 
example, that looks like this: 
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when* djr is ihe measurement uncertainty in the panicle's x position. S/3, is 
its measurement uncertainty in its momentum in t hex direction and ft = ft/ltir. 

That is to s ay, the more accurately you know the posit ion oE a particle, the 
loss accurately you know the momentum, and vice versa. Tills relation holds 
for all three dimensions: 



4vap, ^ | 
& \ 



And the Heisenberg uncertainty principle is a direct consequence ol the 
wauHIke nature of matter, because you can’t completely pin down a wave. 




Quantum physics, unlike classical physics* is completely undeterntmistir. 
Van can never know the precise position and momentum ol a particle at any 
tnit ticnc. You can give only proha bililics Eorlhnsc linked 
measurements. 



Rotting the dice: Quantum physics 
and probability 

In quantum physios, the state of a particle is described by a wave function, 
qif ?; t). The wave function describes the de Broglie wave of a particle* giving 
Its amplitude as a function of position and. time. {.Se^ the earlier section H A 
Dual Ideality: Looking al Particles as Waves" lor snore an tic Broglie.) 




N’ote that the wave lunction gives a particle's amplitude. not intensity:, if you 
want io lintl the intensity ol Hie wave f miction, you have to square it: t ^(r; fyi*. 
Tlie intensity ol a wave is what’s equal to Hie probability that the particle will 
beat that position al that time. 



That's bow quantum physics converts Issues ol momentum and position Into 
probabilities; by using a wave function, whose square tells you the prolwiw'% 
rfrnji/y that a particle Will occupy « particular position or have a particular 
momentum. In other words, li|i(r, I) l = tfris the probability that the particle 
will be found in the volume element ei'r, located at posit inn rat lime /. 
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Besides the position's pace wave function ijnr 0, there's also a mcinsentinn- 
space version of I he ware Junction: I). 

This book, is largely a study of the wave function — the wave Junctions of free 
particles, the wave functions of particles, trapped inside potentials, of identi- 
cal particles hitting' each other, of particles in harmonic oscillation, of light 
scattering from particles, and more. Using this kind of physics, you can pre- 
dict the behavior of all kinds of physical systems. 
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Entering the Matrix: Welcome 
to State Vectors 



In This Chapter 

Creating state vectors 
iJsing Dirac notation fur state vectors 
Working with br.;is gnci kcls 
Understanding matrix mechanics 
Getting to wave mechanics 



■ Quantum physics istiT just nbuoit playing arouml wilh your particle accet- 
^ er&lor while trying not to destroy lire universe. Sometimes, you get to 
do things that are a little more mundane, like turn lights off and on,, perform a 
bit of calculus, or play with dice. 

If you're actually doing physics with those dice {beyond hurling them across 
the room), the lab director won't even, get mad at you. in quantum physics, 
absolute measurements are replaced by probabilities, so you may use dice to 
calculate the probabilities that various [lumbers will come up. You can then 
assembly those values into a vector (sinylc-cuiunui matrix) m Hilbert space 
(a type ol infinitely dimensional vector space with some properties that are 
especially valuable in quantum physics). 

This chapter introduces how you deal with probabilities in quantum phys- 
ics, starting by viewing the various possible states a particle can occupy as 
a vector — a vector of probability slates. From there. I help you familiarize 
yourself with some mathematical notations common, in quantum physics, 
Including bras. ketS. matrices, and wave functions. Along the way, you also 
p|E?t to work with some important operators, 




24 



Part I: Small World, Huh? Essential Quantum Physics 

Creating \trnr Otfn Vectors 
in Hiibert Space 

In quail Lima physics, probabilities take ihu place ul absolute Measurements. 
Say you've been es p^erijneu licty; with rolling a pair of dice and are trying 10 
figure ihe relative probability that tlie dice will show various values. You 
come up with a list indicating the relative probability of rolling a 2, 3, 4 . and 
so on, all the way up to 12: 



Stt M nf ihtf DU*r 
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7 
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In other words, you're twice as likely to roll a 3 than a 3, you're four times as 
likely to roll a 5 than a 2, and so on, You can assemble these relative prob- 
abilities into a vector (If you're thinking of <1 ‘‘Vector’' from physics, think in 
terms or a column oE the vector's components, not a magnitude end direction} 
t« kf-r-jj track of thpm easily: 
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L 

> 

3 

4 
r> 
6 
r, 
4 
A 
2 
L 



Okay, now you're getting closer to the way quantum physics works. You 
have a vector of the probabilities that the dice win occupy various states. 
However, quantum physics doesn’t deal directly with probabilities but rather 
with probability amplitudes, which are tbe square roots of the probabilities. 
To hud the actual probability that □ particle will lie it: a certain state, yon add 
wave Juuttiuiis — which arc going to be represented liy these vectors — anil 
then square them (see Chapter 9 for inlo un why). So take the square rout ul 
;YL these cttlricG to get the probability amplitudes: 

1 

»' = 

3 < 

2 

S'* 

s'- 

5' 

2 

3 

** 

l 



TlClh' 



Brria 




Parr I: Small World, Huh? Essential Quantum Physics 



That's better, but adding the squares of all these should add up to a total 
probability oJ 1 ; as It is now. the sum of the squares oi these numbers is 36, 
so divide each entry by 36'- s , or 6; 
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2':/ 
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jV 
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/6 

r»^/ 

/6 

% 

z' J / 
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$0 m>w ynu can I’d the probability of rolling, airy combination (mm % to 12 by 
reading rfnwn the vector — the probability nF mllin^ a 2 is '6, ol railing a 3 is 
h 

-r~t arid so on, 



Making Life Easier With Dirac Notation 

When yon have a state vector that gives the probability amplitude that a pair 
of dice will be in their various possible states, you basically have a vector in 

dice space — all the possible states that a pair ol dice can take, which is an 
1 1-dlmensional space. (See the preceding section for more on stale vectors.) 
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But in moat quantum physics problems, the vectors can be infinitely large — for 
example, a moving particle can be In an infinite number of states. Handling 
large arrays of states isn’t easy using vector notation, so Instead of explicitly 
writing out the whole vector each time, quantum physics usually uses the 
notation developed by physicist Paul Dirac — the Dirac or (mr-toer nafttffon. 



Abbreviating state Vectors as kets 

Dirac notation abbreviates the slate vector as a ke J, Like ibis: ip =■. So Ln the 
dice example, you can write the state vector as a kut this way: 
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Here, the components of the state vector are represented by numbers in 
1 1 -dimensional dice space. More commonly, however, each component rep- 
resents a function, something like this; 
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You can Lise functions as components of a state vector ns long as they’re 
Linearly independent functions (and so can he 1 rented as independent axes 
in Hilbert space). In general, a set of vectors $ v in Hll her! spare is linearly 
independent if the only solution to rhe follnwirii* equation Ls that all the coef- 
ficients Q a - 0: 



y (i o = o 

i ■ ■ 



That is, as long as you can't write any one vector as a linear combination ol 
the others, the vectors are Linearly independent and so form a valid basis in 
Hilbert space. 
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Writing the Hermitim conjugate as a bra 

For every k:et„ there's a corresponding dro. (The terms come from b*a k?i. or 
frmcftef, which should be clearer In the upcoming section titled “Grooving 
With Operator*.*) A tut? is the Hern I dan conjugate of the corresponding ket 



SL]|j]]use you start with this ket: 
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The ' symbol means the complex conjugate. (A complex conjugate Hips the 
sign connecting the real and Imaginary parts of a complex number.) So the 
corresponding bra, which you write as ■■■ d I , equals l y> r '. The bra is this row 
vector: 
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Note that it any ol the elements of the ket are complex numbers, you have to 
take their complex conlugale when creating the associated hra. For instance, 
it your complex number in the feet i* fl t to. its complex conjugate in the bra is 
a — to. 
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Multiplying bras anti kets: 

A probability of 1 

You ran taki» the product al your kel arid lira, denoted as -=t|j ! ip>, like tliis: 
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This is just matrix multiplication, and tlie result is the same as taking the sum 
of the squares of the elements' 




< ur u , > = - 1 - + A + A + J . + - r> - - A + A + A . + A + A +. .L = t 
:«i :u; :tn 3G :w :ik mi :u; :« hg :w 

And that's the way it should be, because the total probability should add up to 
8. Therefore, in general, the product of the bra and ket equals I: 

C IV IV > = 1 



If this relation holds, the ket I ip =- is said to he najmiilizad. 
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Covering alt your (rases: Bras and kets 
as basis-less state Vectors 

The reason ket notation, ip, Is- so popular In quantum physios is that It 
allows you to work with state vectors in a basis-free way. hi other words, 
you're not stuck in the position basis, the momentum basis, or the energy 
basis, That’s helpful, because most of the work in quantum physics takes 
place In abstract calculations, and you don't want to ti&vfi (0 dray: all the com- 
ponents oi stale vectors through ‘hose calculalious (oltcu you can't — there 
may be icifiriiu- possible .states in the problem yuu're dealing with). 

For example, say that you’re representing your states using position vectors 
in a throe-dimensional Hilbert space — that Is, you have *. y, and z axes, 
form Inga position tors for your space. That's fine, but not all your calcula- 
tions have to be done using that position basis. 

You may want to, for example, represent your states in a three-dimensional 
momentum spate, with three- axes in Hilbert space, p . p , and p.. Now you’d 
have to- change all ymir position vectors to [in.:iii witimi vectors, adjusl- 
uig each component, and keep track of what happens to every component 
through all your calculations. 

So Dirac's bra-ket notation comes to the rescue here — you use U to- perform 
all the math and then plug in the various components of your state vectors as 
needed at the end. That Is. you can perform your calculations in purely sym- 
bolic terms, without being tied to a basis. 

And when you need to deal with the components ul a kut. such as when you want 
to get physical answers, you can also convert, kets to a different basis by taking the 
kefs components along the axes of dial basis. For example, if you want to convert 
the ket to the position basis, as represented hyf,/, and ft, which are posilioa- 
unit vectors along the x r y. and z axes, you can just find the three components of 
( ip along ij, and k for the new version of the ket, l+>. Here’s how that looks in 
general, where ^ are unit vectors in the bmm you’re switching to; 

Understanding some relationships 
using kets 

Ket notation makes the math easier than it is in matrix form because you can 
take advantage of a few mathematical relationships. For example, here’s the 
so-called Sehwara inequality lor slate vectors: 



>pyris 
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|-- ij |$ >| < < 1 / 1 cpr ijt > 

This says tlizit rhc square of the absolute value oE the product of two state 
vectors, ! is Less than or equal to *i|i I qr This turns out the he 

the analog ol the vector inequality: 

A B 's, A B|'' 

So why is the Schwarz inequality so useful? It turns out that you ran derive 
the Heisenberg uncertainly principle from it (see Chapter l for more on this 
principle). 

Other ket relationships can also simplily your calculations. For instance, two 
kets, I and art: said 1,6 be OrtlrGgunai it 

^ \y j > = 0 



And two kets are said to hr: mtimntmnal if they meet the following conditions: 
C •• I I 

v* = l 

^ = J 



With this Informal ion In mind, you're now ready to Start working with operators 



Grooving With Operators 

What about all the calculations dial you' re supposed to be able to perform 
with kets^ Taking Lite producl of a bra and a kel, ^ip I -p.’, is line as lar as it 
goes, hut what about extracting some physical quantities yon can measure? 
That's, where operators come in. 



Hello, operator: HoW operators Work 

Here's the general definition of an operator A in quantum physics: An opern- 
for is a mathematical rule that, when operating on, a f:r-T, hp>, transforms that 
ket into a new ket. I yV in the same space (which could just be She old ket 
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multiplied by a scalar). So when you have an operator A, It transforms kets 
like this: 

A ■/> '" > 



For that matter, the same operator can aLsit tnuisiurm bras: 

< i; 1 ' A ==< y 

Here are several exam pies of the kinds of operators you'll see: 

v* Hamiltonian (H): Applying the Hamiltonian operator (which looks diller- 
ent tor every different physical situation) ylves you E, the enerijy of the 
particle represented by the ket E is a scalar quantity: 

1 1 w > - E it -j 



t* Unity os* identity (I)* The unity or Identity operator leaves fcets 
line handed: 

I v > * y '> 

t* Gradient (?); The gradient Operator works like this: 

u i p>j + 4 <?>* 

1 ibr *Jy 1 Jz 1 

** Linear momentum (P): Tin? linear monietUum operator looks like this in 
quantum mechanics: 



F ip > - - iTi V y > 



V Laplaclan (A): You use the l-aplacfen operator, which is much like g sge- 
nnd-arder Radiant. to create the finarpy-iindin^ Hamiltonian operator: 

A w > a A y > = y > ! ■ Jp 1 . 1 ty > ■ JO. y j> 

ax ■ 1 oy ‘ 1 ilz" 1 
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In general, multiplying operators together is nut the same independent of 
order, sn Eor the operators A and E, Al-l ± BA. 

And an operator A is said rn be linear LF 

A ' i ■ (y > +c. ^ > J = c, A v- r > h-t A 1 / > 



/ expected that: Finding 
expectation Values 




Given that everything in quantum physics is done in terms ul pri:h abilities, 
making predictions becomes very important. And the biggest such prediction 
is list: expectation value. The expectation value ol an operator is the aver- 
age value that you would measure it you performed the measurement many 
times. For example, the expectation value ol the Hamiltonian operator (see 
the preceding section) is the average energy of the system you're studying. 

The expectation value is a weighted average of the probabilities ol the sys- 
tem’s being in its various possible states. Here’s hew you find the expectation 
value of an operator A; 

Expectation value =<^ A ip> 



Note that because you can express l as a row operator and fhp as a 
column vector, you can express the operator A as a square matrix. 

For example, suppose you’re working with a pair of dice and the probabilities 
of all the possible sums (see the earlier section "Creating Your Own Vectors 
in Hiibert Space"]. In this dice example, the expectation value is a sum of 
terms, and each term is a value that Can be displayed by ihe dice, multiplied, 
by the probability that that value will appear. 

The bra and ket will handle the probabilities, so it s up to Ihe operator that 
you create for this — call it the Roll Qf/erotor, Ft — to store Ihe dice values 
{2 through 12) for each, probability. Therefore, the operator R looks like this; 
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R =1 

21)0 0 (3 0 □ 0 0 I) 0 D 

0 300 0 0 0 Cl 0 0 0 0 

0 0 400 0 D 0 (3 I) 0 U 

000 509 000000 
0 Cl 0 0 600 0 0 0 0 0 
0 0 0 0 0 700 0 0 0 0 
0 0000 0 800 0 0 0 
0 0 □ (3 0 D 0 0. 1) 0 Q 

0 O 0 0 0 0 0 0 10 0 0 

0 0 0 (3 0 LI 0 0 1)11 U 

0900 0 0 000 013 



So to find the expectation value of R, you need to calculate <i?IRi ip. Spelling 
that out In terms of components gives you the following: 
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Dumg the math, >'iJu get 
< iv R > = 7 
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So the expectation value of a roll oi the dice Is 7. hew you can see where the terms 
fern and. Aw cons? from — they 'bracket' an operator to give you expectation 
values, Irt fact, the expectation value Is such a common thing to f ind that 
you’ll olten find <u I R hp> abbreviated as *R:<, so 

--.Rv - 7 



Looking at Umar operators 

An operator A is said m he tinew it it meets the following condition: 

A(c. I i|r> * Cj ht>) = £',Al V 3 " r tAl 

Fur instance, the expression lo^yi is act unify a linear operator. To see that, 
apply l<i':..-yr to a feet, lx-*: 

!*-■< v\?r* 



You can also write this as 
< v x > $> 

The expression <nrlx> b always fl complex number (which could be purely 
real), so this breaks down to 

fjo> 



where c is a complex number. Thus. loxf is Indeed a linear operator. 



Going Her mi ti an With Hermitian 
Operators and Adjoints 

The NermMQn adjoint — also called the adjoint orHemitUm conjugate — of 
an operator A is denoted A r . To find the hermetian adjoint, lollow these 
steps: 

I- Replace complex cotutonl* with Ihoir complex conjugate*. 

The Hermitian adjoint oEa complex number is the complex conjugate of 
that niimbsr: 



>pyr qhte 

r J O - 
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2- Replace kefs with (heir corresponding bnts, and replace bra* with 
their corresponding kets- 

You have to exchange the bras and keti when finding the KermJHan 
adjoint of an operator, so Finding the Ftermltlan adjoint of an operator fa 
not just the same as mathematically Finding Us complex conjugate. 

3. Replace operators with their I lenni than operators. 

In (jui 1 1 1 1 u]ii mechanics, Operators l Fiat are equal to Ifmir HermUia.ii 
adjoin ts are called ttermilian opemion. In other words, an operator is 
l-lcnnillan if 



A 1 = A 



HermitLan operators appear throughout; tJle hook, and they have spe- 
cial properties. For instance, the malm that represents them may be 
diagonalized — that is, written so that the only nonzero elements appear 
along Hie matrix's, diagonal. Alsu, the expectation value of a I termilian 
operator is guaranteed to be a real number, not complex (see the earlier 
section "l expected that; Finding expectation values"), 

4 - Write your final equation. 




< fff A 0 5 = c rt A 



Here are some relationships cnnceming HermitLan. aiijnints: 
^ |cdA| sfl'A 1 



^(A')' = A 
^ (A < Q}' ■A 1 +B r 



V* (AB|^>y =<tpp|B" A f 



Forward and Backward: 
Finding the Commutator 




The measure of how different It Is to apply operator A and then p, versus B 
and then A, is rated the operators' rommjjfrrfnr. Here's hmv you defiisf! the 
commutator nf operators A and 13 : 



l A. 11] = Ad - UA 
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Commuting 

Two operators fonrnjwte with each other f their' commutator is equal to zero. 
That is, it doesn't make any difference In what order you apply them: 



I A. BJ =o 



Note in particular that any operator coninmfes with itself: 



[A. A] =0 



And It's easy to show that the commutator of A, B is the negative of the com- 
mutator of B, A: 



f-V B1 - -[ B f A I 



It’s also true that commutators are linear — that is. A(c, i qf> - C, l'4'>) = C| A I y» 
+■ c.,Ai 



[A. B + e * D - ...] = [A. B] * [A. C\ ♦ [A. Dj + ... 



And the Hermit ian adjoint of a commutator works Lins way: 

| a,b[ -[b'.a 1 ] 



Vcm can alaolind the anticomimitator, [A, I4J-: 

(A. llj = AB + BA 

Finding anti-Hermitian operators 

Here's another one: Wli&t can you say about the I tenuitian adjoint oi the com- 
mutator of two I lermitian operators.? Here's the answer. First, write the adjoin L 

La. b]' 



The deEinitiun oE romiriLLtatnrs tfllls you Ihe FolLnwin^: 
| A, IS ]' =(AB-BA} r 



You know r AB)' ■ BW (see the earlier section "Going Hei-m Ilian with 
Hei-mltian Operators and Ad joints” for properties of adjoin is). Therefore, 

[ A.B]' =( AB - BA)' =B'A t " A'B 1 



TlClh' 
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But [or Hermltian operators, A - A r , so remove the ’symbols: 
| A, B | = ( All - BA | ' = B 1 A r - A 'B r = BA - AB 
But BA - AB is gust -I A. B], su you have the lulJuwiqjf 

[a.b| - -[a, bJ 




A and R here are Hermltian operators- When you take the HermLtJan adjoint of 
m expression and ii|et the same thing back with a negative sigp in Iront oi it. 
the expression is called anti-Hemition, so the commutator of two Hermltian 
operators is gntk-Heinritian. {And hy the wav, the expectation value of on anti- 
Honniliiin operator is guaranteed to be completely imaginary.) 



Starting from Scratch and Ending Up 
ufith Heisenberg 

If you've read through the last Tew sections, yoti're now armed with all this 
new technology: Hermit ion operators and commutators. Hr>w can you put it 
to work/ Ynu can rn>me up with the Heisenberg nncnrlninly relation starting 
virtually from scratch. 

Here's a calculation that takes you from, a few basic definitions to the 
Heisenberg uncertainty relation. This kind of calculation shows how much 
easier it is to use the basis- less bra and kel notation than the full matrix ver- 
sion of state vectors. This isn’t the kind of calculation that you’ll need to do 
In class, but follow it through — knowing how to use kels, bras, commuta- 
tors, and Hermltian operators is vital In the coming chapters. 

Thfi unrerlainty in n measurement of the Hermltian operator named A Ls for- 
mally given by 

AA = (< A’ > — < A > " I 

That is, AA is equal to the square root of the expectation value of A- minus 
the squared, expectation value of A. ff you’ve taken any math classes that 
dealt with statistics, this formula may be Familiar to you. Similarly, the uncer- 
tainty In a measurement using Hermltian operator B Is 

AR= ;<B : >- 
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Now consider rtip operators aA anti AB (not the uncertainties aA amt AtS any- 
more), ami assume that applying aA and. Ati as operators gives you measure- 
mml values like this: 

,-\A = A - < A > 

AC 4= C - B > 



Like any operator, using A A and AES can result in new Jfeets: 

AA y:>= x > 

AC y>=|p> 

tlert-fi the key: Tlie Schwarz inequiiility (from the earlier section 
"Understanding some relationships using tots") gives yon 



So you can sen; tliat the inequali ty sign, 2, which plays a big part In the 
Heisenberg unocrtaiEity reLaliun. has already crept inio the calculation. 

Because aA and aB are Hermitian. <% \%> is equal to njj l aA'-' and i$> is 
equal to m(j iaIB' ! Because AA' = AA (the delinitioti ol a Hermitian operalot), 
you can see that 

■i ;> = ■: u/AA AA ^ 

TJi is means that 

/I y > = < A A AA w > — < uAA iy > 



That is, *y_ is equal to ;AA-> and -^o lq> > is equal to So you can 
rewrite the Schwarz inequality like this,: 

< AA* x AB' ^>k AAAB 



Okay, where has, this gotten you? it's time to he clever. Note that you can 
write .‘YAaB as 

AAAll = ■' [ iA,Ali]+|{AA,AB} 
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Here, [iA, Ah] = aAaB + A IMA is the anticcmuiiurator of the operators &A and 
Alt. Because [aA, AJlj = [A, 13 J (tlie constants '~A> and -Lis subtract out), you 
can rewrite this equation: 

AAAB = -±[a,b] .[] AA. AB| 

Here's where the math gets intense. Take a look at what you know so far: 

t* The commutator of two Hermltian operators, I A, R], is anti-Hermittan. 
p* The expectation value of an anti-Hermitian is imaginary. 
p* -aA, AJJJ is hermit ian. 
m* The expectation value of a HumiitLan is real. 

AN this means that you am view the expectation value of the equation as the 
sum of real ({AA, AB)) and imaginary f [ A. R|j) parts, so 

And because the second term mi the right is positive or zero, you can say 
that the tallowing is true: 

|< VAR -■ 1 A'||-=:[A B ]:=■ 



Whew] Rut now compare this equation to the relationship From the earlier 
use of the Schwarz inequality: 

< AA* sc AB J > > |< BAAA 

Combining the two equations gives you this: 

■; AA* AB J >■> I|< A_ B 

•I I I 

This has the look of the Heisenberg uncertainty relation, except ior the pesky 
expectation value brackets, -- >, and the Fact teat aA anti aB appear squared 
here. You want to reproduce tee Heisenberg uncertainty relation here, which 
Lunks like this: 

AxAp 2 - 

2 
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Okay, so how (]<] you get The lelt sliLe> of the equation Irom <aA- to 

aA^U? Because an earlier equation tells you rtiat A A = A - «A*, von know the 
following:: 

■■■: AA ‘ > = ■=: A : - A > 2 - 2 A< A > -• 

Taking tin: expectation value the last term in this eq nation. you get this 
result:: 

< AA' i=<A'KA?'- 2A<A>* = <A'j-cA> ; 



Square Hie earlier equation A A = («A‘* - 'A-»-) lc to get the following: 
AA* =e A J >-< A>' 

And comparing that equation to the before it, you conclude that 



< AV > = AA-' 

Cool. Thai resii 1 1 means tliat '' vV >=: AR > i * [ A ' R 1 
AA ; AB ; > | -v f A. 1:5 ] >| 



becomes 



This inequality at last means that 
AAAB a l -:| \ b;>| 

Well, well, well. So the prod net ol two uncertainties is greater than Or equal 
to L /i the absolute value ol the commutator of their respective operators? 
Wow. Is that the Heisenberg uncertainty relation? Well, take a look. In quan- 
tum mechanics, the momentum operator looks like this: 

i" 1 ■ itiV 

And the operator for the momentum in the x direction is 




So what's the rnmnuiLtror of the X operator (which just returns the x posi- 
tion of a panicle) and I s ? |X. P, | = —/ft. so from AAAB 5 | < | .A. B J>|. you gr?t 

u 



'r- 
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this [text equation (remember. if arid i\p, here are tlie uncertainties in x and 
.ip . not I Ilc operators): 

A*A P„ 5 ^ 

Hot dog! That h the Heisenberg uncertainly relation. (Notice that by deriving 
it Emm scratch, however, yon haven't actually constrained The physical world 
through the use of abstract mathematics — you've merely proved, llh big a 
few basic assumptions, that you can’t mausum the physical wurhl with per- 
fect accuracy.) 

Eigenvectors and Eigenvalues: 

They’re Naturally Eigentastie! 

As you know if you've been following along in this chapter, applying an oper- 
ator to a kjet can result in a new ket: 

A ^> = 1* > 



To make things easier, you can work with eigenvectors and eigenvalues few?. 1 ? 
is German for "innate" or "natural"). For example, i is an eigentwteroi the 
operator A il 

I p* The number fj i:i a corn plus constant 

Note w hat's hap petting here: App lying A to one of its eigenvectors , ip, gives 
you tp hack, multiplied by that eigenvector's etgenoalue. a. 

Although o can be a complex constant, the eigenvalues of Hetmitlan opera- 
tors are real numbers, and their eigenvectors are orthogonal {that is, 

= «)- 

Casting a problem in terms of eigenvectors anil eigenvalues can make life a 
lot easier because applying the operator to its eigenvectors merely gives you 
the same eigenvector back again, multiplied by its eigenvalue — there's no 
pesky change of state, so you don't have to deal wilh a different slate vector. 
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Take a look at this Idea, using the R operator from rolling the dice, which Is 
expressed this way in matrix form (see the earlier section M l expected that: 
Finding' expectation values" for more on this matrix); 



2 0 0 0 0 0 0 

0 3 0 0 0 0 0 

() 0 1 () 0 0 0 

0 0 0 5 0 0 0 

0 0 Cl ft 0 0 0 

0 0 0 0 0 7 0 

ft 0 0 0 0 0 A 

ft 0 0 0 0 0 0 

ft 0 0 ft 0 0 0 

ft o n ft o □ o 

0 0 0 0 0 0 0 



ft ft 0 0 
0 0 0 0 
ft ft ft n 
ft Cl 0 0 
ft 0 ft Q 
ft ft ft 0 
ft 0 ft 0 
ft 0 ft 0 
0 10 ft 0 
ft Oil □ 
ft 0 012 



Tile R operator works in 1 L-dimensional space and is Elermilian. so there'D he 
II orthogonal eigenvectors and 11 corresponding eigenvalues. 



Because R Is a diagonal matrix, finding the eigenvectors is easy. You can take 
unit vectors in the- eleven different directions as the eigenvectors. Here's 
what the first eigenvector, would Look like: 



£,= 

i 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 



righ 
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Aj]d here's what the second eiy;enver.tor. c.,, would look like: 

S= = 

0 
1 
0 
0 
0 
rj 
0 
0 
(I 
0 
v) 

And so on, ujj to 4 n ; 

b 
s-i 

l> 

t) 

l> 

0 

0 

0 

0 

0 

[) 

(I 

] 

Note that all the eigenvectors are orthogonal. 

And the eigenvalues? They're the numbers you get when you apply the R 
operator Lo on eigenvector. Because the eigenvectors are just unit vectors in 
all L I dimensions, the eigenvalues are the numbers on the diagonal of the R 
matrix: 2, 3. 4. and so on, up to 12, 
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Understanding how they Work 




The eigenvectors ol a Herm Ilian operator define a complete set of orthonorraal 
vectors — that Is, a complete basis for the state space. When viewed In this 
"eigenbasks, h which Is built of the eigenvectors, the operator In matrix format is 
diagonal and 1ln> elements along the diagonal of the matrix are the eigenvalue!!. 



Tliss arrangement is one of the main reasons working with eigenvectors is so 
useful; your original operator may have looked something Like tliis (Note: Bear iri 
utintl that tliC dementi in an Operator call also Lie functions, not just numbers'): 

R = 

0 0 1 0 6 0 0 3 0 0 00 

0 1 D 0 0 S (5 0 0 0 0 

4 0 0 0 120 0 2 0 0 0 
9 0 0 0 0 0 7 0 0 0 0 

0 0 6 0 6 0 0 0 8 0 0 

1 0 30 S 00 0 0 0 S 

oo omi ohoo n 

0 0 0 0 0 0 0 9 0 0 0 

0 0 □ 0 !) it 0 (1 0 o 77 

0 0 0 0 L10 0 0 0 I 0 
0 0 7 0 0 0 880 0 0 D 



By switching to the basis of eigenvectors for the operator, you diagonalize 
the matrix Into something more like whal you’ve seen, which is much easier 
tu work witll: 

R = 

200000 00000000 000000 00 
030000 0 0000 
0 0 10 0 1 ) € 0 0 0 0 
0006 o 0 o o ti o n 
0 0 0 0 6 0 0 0 0 0 0 

0 0 0 0 ((7 00 (I 0 0 

0 0 0 0 0 0 6 0 0 0 0 

0 0 0 0 0 0 0 9 0 0 0 

0 0 0 0 0 0 0 0 LOO 0 

000000000 MO 

o r> o a oooo o o n 
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Yon can see why the term ctgm ta applied to eigenvectors — they form a 
natural basis for the operator. 

If two o r more of the eigenvalues are the same, that eigenvalue is said to be 
degenei me. So for example, if three eigenvalues are equal to 6, then the eigen- 
value 6 is threefold degenerate. 

Here's another cool thing: If two Hermit ian operators, A. and fl, con mule. and 
ii A doesn't have any degenerate eigenvalues. then each eigenvector ui A is 
also an eigenvector of B. (See the earlier section “forward and Backward: 
tin ding the Commutator" lor more on com muling.) 



Finding eigenvectors and eigenvalues 

So given an operator in matrix form, how do you find Its eigenvectors and 
eigenvalues? This is the equation you want to solve: 

A y > = a v -■ 



And you can rewriLe this equation as the following: 



[A- at) ^> = 0 

t represents the identity matrix, with Is along its diagonal and (Ks otherwise: 



r = 

L 0 0 0 0 0 0 0 

oi on o o n o 
0 0 10 0 0 I] 0 
i> U IJ ] 0 tt u u 

0 0 Cl 0 I 0 0 0 

0 U U 0 0 1 Cl 1) 



0 0 ... 
0 (I.„ 
0 

0 (S... 

0 0-J 
0 0 ... 



The solution to (A - ol) u - - 0 exists only if the determinant of the matrix A 
ol Is t>: 



del(A-nl) =ti 





Parr Pi Small World, Huh? Essential Quantum Physics 



Finding eigettVafues 

Any values of o Lhat satisfy Hie equation det(A - d) * 0 are eigenvalues of the 
original equation. Try to lind the eigenvalues and eigenvectors ol the Colin w- 
ing matrix: 




First, convert the matrix Into the lorm A - til: 



A - a\ = 



-1- a -I 
2 -\~a 



Nest, find the determinant: 

det(A -cjll = (-1 - a X-4 - «} + 2 
det(A-al) = hj- * 5a - h 

And this can be (adored as follows: 



det(A - oi) - q- + 3o * 6 - (a - 2 )(Yj + 3) 



You know that det(A- rtl) ■ 0, so the eigenvalues ol A are the roots of this 
Equation; namely. o. = -2 and a = -3. 



Fiitding eigenvectors 

How about Ming the eigenvectors? To find the eigenvector corresponding 
to cr (see the preceding section), substitute a, — the llrst eigenvalue, -2 — 
Into (he matrix In the form A - air 

. . l-l-fl -I 

-4 -a 

-3 

-2 




I A- of | = 0 



So you have 



l -1 




mi 


0 


2 -1 


Vi 




0 
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Because every raw of this matrix equal icm must he true, you know that qr, = 
i.5,i . And teat moans I hat, up to ail arbitral) 11 constant, the; eigenvector corre- 
sponding to a Is the lolJowi[LL£ 




Drop the arbitrary constant, and Just write this as a matrix: 

L 

1 

How about the eigenvector corresponding to cl? Plugging «l, -3, into the 
matrix In A -trl form, you get the follow lug; 

\-ul=2 -ll 
2 ~L 



Then you have 



2 


-1 




= 


0 


2 


-L 


V, 




0 



So = 0, and V, = 1H-, * 2, And that means that, up to 50 arbitrary con- 

stant, t ho iiicjenvec tor corresponding too, i$ 



r 



] 

2 



Drop the arbitrary constant: 
I 

■t 



So the eigenvalues ol Lius next matrix operator 




-t 

-I 



r no h t 



materia 
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are a, * *»2 and a., * -3. Anti she eigenvector corresponding to a, is 

1 

1 



The eigenvector corresponding to a, is 

f 

2 

Preparing for the Inversion: Simplifying 
(Pith Unitary Operators 

Applying the inverse of an operator undoes, the work the operator did: 

A- 1 A . AA-* = I 



Sometimes, finding the inverse ol an operator Is helpful, such as when you 
want to solve equations tike Aj- = y. Solving for x is easy if you can find the 
I nve rs e of A: jr = A y. 

However, finding the inverse of a large matrix often isn't easy, so quantum 
physics calculations are sometimes limited to working with unitary opera- 
tors, |j, where the operator's inverse is equal to its adjoint, U 1 = 11 f . (To find 
the adjoint ol an operator. A. you find ihe transpose by interchanging the 
rows and columns. A 1 '. Then lake the complex conjugate, A 1 " = A 1 .) This gives 
you the following equation: 

UTJ. UU r . I 

The product of two unitary operators. U and V. is also unitary because 



f u v }( vn )' =( y v }( v 'y' ) - u ( w 1 ) u T = uu ' - 1 



When ynuustt unitary operators, kels and liras transform this way: 
y/ > - < vr| L" 1 
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And you can transform other operators using unitary operators Like this: 
,V = UAH ' 



Note that the preceding equations, also mean the following: 



^ |y> = SJ' Y> 
v\ = < V 




v* A = VI ’A t? 

Here are some properties of unitary transformations: 

I t* If an operator Is Hermitlftn, then its unitary transformed version. A* = 
U AUt is als O Hermirlan , 



t* The eigenvalues Of A and its unitary transformed version, A = UAU 1 . are 
the same. 



t* Commutators that are equal to comp lex numbers are unchanged by uni- 
tary transformations: [A', B ] = [A, B]. 



Comparing Matrix: anti Continuous 
Representations 

Werner Heisenberg developed the matrix-oriented view of quantum physics 
that you've been using so Far in this chapter. It's sometimes called matrix 
mechanics. The matrix representation is fine for many problems, hut some- 
times you have to go past It. as you're about to see. 

One of the central problems of quantum mechanics is to calculate the energy 
levels of a system. The energy operator Is called the HmUUioman, H, and 
finding the energy levels of a system breaks down to finding the eigenvalues 
of the problem: 

H ijf > = E iy ;■ 

Here. E is an eigenvalue of thelf operator. 
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Here's the same equation in matrix terms; 



H„-E 


", 






-0 


H, 




H, 


H ,, 




H, 


IU 


H„-E 


h m 




H,. 




H, 


H.-E ... 





The allowable energy levels of the physical system are the eigenvalues RL 

That’s fine if you have a discrete basis of eigenvectors — if the number of 
energy states is finite- But what if the number Of energy states is infinite? In 
that case, you can no Longer use a discrete basis For yniir operators and hras 
and kels — you use: a amiinuous basis. 



Going continuous With cahutus 

Representing quanto nr mechanics in a continuous basis is an invention of the 
physicist Erwin Schr&dinger. In the continuous basis, summations become 
Integrals. For example, take ihe following relation, where I is the identity 
matrix: 



It becomes the following: 



J tJ". o x o. I - I 



And every ket ! ip» can be expanded In a basis o! other kels. I o„>. like this: 



w > i J jJrr p •;>' > 



Doing the wave 

Take a look at the position operator, P, in a continuous basis. Applying this 
operator gives you r, the position vector: 

R m ? 3 p- y> 



TlClh' 
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In this equation, applying the position operator to a state vector returns the 
locations, r. that a particle may he found at. You can expand any Tcet In the 
position basis like (his: 

I,/ > — 1 d 'r|r x i ■ nr > 

And this becomes 



w> I K^[r)| r r ly > 




Here's a very Important thing to understand; q(Y) = i|) j is, the tmst function 

for the state vector ip - — it's the feet's representation In the position basis- 
Or in common terms, it’s just a function where the quantity I i|r(r}l-(# 5, j' repre- 
sents the probability that the particle will he found In the region rfVat jr. 



The wave function is the foundation of what’s called wave mechanics, as 
opposed lo matrix mechanics. W hat's important to realise is thal when you 
talk about representing physical systems in wave mechanics, you don'l use 
the tasteless bras and kets of matrix mechanics; rather, you usually use the 
wave fund ion — that is. bras and feels in the position basis. 



Therefore, you cjo From talking about I to --rl ip>, which equals y£r), This 
wave function appears a lot in the coming chapters, and it's Just a ket in the 
position basis. So in wave mechanics, H I i*i> = E ip becomes the following; 



< r | H ¥> > = F < r i;/ > 



You can write this, as the following; 
< r|H ip =■ ■ E i, 11 1 r | 



But what is <r\ 1 1 i ipr>? it .$ equal lo lte(r). The 1 lam i I Ionian operator. 1 1, is lihe 
total energy of the system, kinetic (p-!2m) plus potential (V(r)) so you get the 
Following equation: 
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But l he momentum. operator is 

P 1 w -■ * ur i j { +-ih 3_; v / > j + ift ‘if > k 

tor 1 Dy 1 dz 1 



Therefore, substituting the momentum operator lorp gives you this: 

=& 1 .\ 



H = 



'hn 



JT , £_ + 

{ r.Lv dy ■>=- 



4 V 



I/) 



Us ini* the Laplacian operator., you u]et this equation: 



You eon rewrite this equation as The folLniving r called the. Schmrlingpr wjiuifiany. 

H vl/ 1 = 4 v{ r ] + V f r ) y \ r J = j 



So in the wave mechanics view of quantum physics, you’re now working with 
a differential equation instead of multiple matrices of elements. This all came 
from working in the position basis, y(r) - <r\ y.- instead of fust i y>. 

The quantum physics to. the rest ol the book ts largely about solving this dif- 
ferent 1*1 equation for .% variety of potentials, V{r>, That is. your focus Is on 
Finding the wave function that satisfies the Schrodinger equation for various 
physical systems. When you solve the Schriid inger equaritm lor f(r), you can 
Find the allowed energy states for a physical system, as well as l Ilk probabil- 
ity that the system will he in a certain position state. 



Note that, besides wave (unctions in the position basis, you can also give ts 
wave function in the momentum basis, if ip), or in any number of oilier bases. 




The Heisenberg technique of matrix mechanics is one way of working with 
quantum physics, and It’s best used for physical systems with well-defined 
energy states, such as harmonic oscillators, The Schridlnger way of looking at 

thing. 1 ;, wave mec hanks, uses wave Functions, mostly in flic positm-n basis,, Ip 
reduce questions in quantum physics to a differential equation. 




Part II 



Bound and 
Undetermined: 
Handling Particles 
in Bound States 



The 5 th Wave Ely Rich Tennant 




Aion^ with. 'Antimatter/ 1 and 'Dark Hatter/ 
we h VE recently discovered- the existence of 
'l^esn't Hatter/ which appears to have no 
effect on the universe whatsoever.* 






In this part . . . 



M his part is where you get the lowdown on one ol quan- 
w turn physics' Favorite topics: solving the energy Levels 
anil wave kmc lions lor particles trapped in various bound 
states. For example, you rnay have a particle trapped in a 
square wuJI. which is itiuch like having a pea in ahox. Or 
you may have a particle in harmonic oscillation. Quantum 
physics is expert at handling those kinds ol situations. 






Chapter 3 






h This Chapter 
Winders Handing potential wells 
Working with infinite square wells 
IDete running enertfy Bevels 

Trapping particles with potential barriers 
I landling 6ree particles 



MM/ hat's that, Lassie 1 ' Stuck in an energy well? Go get help! Iil this cliap- 
Vr ter, you get to see quantum physics at work, solving problems in 
one dimension. You see particles trapped in potential wells and solve for the 
allowable energy states using quantum physics. Thai goes against the grain 
In classical physics, which doesn't restrict trapped particles to any particular 
energy spectrum. But as you know, when the world gets microscopic, quan- 
tum physics takes over. 



The equation of s he moment is the Srhiddiiigrr equation (derived iri Chapter 2 ). 
which lets you solve tor the wave fund km, <|r(jr)„ and the energy levels, E: 



^^fO^V( J >{r) = K V r(^ 



Looking into a Square WeU 

A Sierra .if uwfl is a pOtcnrial (thal is. a potential energy welt} that forms a 
square shape, as you can sec in Figure 3-i. 
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The potential, or goes tu infinity as x = 0 and x * a (where a is distance). 
Like 1 1 Lis: 



I v* V'U) = where* < 0 
/ V(j) = 0, where 0 < x £ a 
v* V(x) - where x > a 

Using square wells, you cam rrap particles- It you put a particle Into a square 
well with a limited amount of energy, it'll be trapped because it can’t over- 
come the infinite potential at either side ot the square well- Therefore, the 
particle has to move inside the square well. 

So does the particle just sort oE roll around on, the bottom of the square well? 
Not exactly. The particle is in a bound state, and its wave function depends 
on its energy. The wave function isn’t complicated: 




n = l, 2, 3 . 
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So you have the allowed wave functions lor the stales n • 1, 2.3, and so oiv 
The energy of the allowable bound stales are given by the following equation: 

E B iX,r n- 1.2.3... 

2m« 2 

The rest of this chapter shows you how to solve problems like this one. 



Trapping Particles in Potential Wells 

Take a look at the potential in Figure 3-2. No Lice the dip. or u&H. in the poten- 
tial. which means that particles can be trapped to it if they don't have too 
much energy. 

The particle's kinetic energy summed with its potential energy is a constant, 
equal tn its. total energy: 



If iis total energy is Less than V,, the particle will be trapped in the potential 
well, you see in Figure 3-2; to get out ol rhp well. The particles kinetic energy 
would have to become negative lo satisfy the equation, which is impossible. 
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In this section* you take a look at the various possible states that a par' 
tide with energy E can take In the potential given by Figure 3-2. Quantum- 
median Ically speaking, those states are of two kinds — bound and unbound 
This section looks at them in overview. 



Binding particles in potential Wells 

Hawd states ha open when the particle isn't tree to travel to infinity — it’s as 
simple as that. In other words, I he particle is cuulined to the potential well. 

A particle i raveling in the potential well you see in Figure 3-2 is bound if 
Lis energy, E. is less than both V ( and v g . in that case, the particle moves 
between x t and ,r ; . A particle trapped in such a well is represented by a wave 
Function, and you can solve the SehrGdinger equation for the allowed wave 
functions font the allowed energy stales, You need to use two boundary con- 
ditions d he Schrodinfeer equation is a second-order differential equation) to 
solve the problem completely- 




Bound slates bur dtsrmte — that is, they form an energy spectrum of discrete 
energy Levels. The Bchrodinger equation gives you those states. In addition, 
in one-dimensional problems, the energy levels of a bound state are not 
degenerate — that is, no two energy levels are the same in the entire energy 
spectrum. 



Escaping from potential Wells 

If a particle’s energy. E, is- greater than the potential V, In Figure 3-2. the par- 
ticle can escape from the potential well. There are |wo possible eases: 

V, K. < v, and ¥. -- v,. This section loots at them separately. 



Case ?; Enemy between the two potentiate W } < F< 

If v, <E< V 2I the particle in the potential well has enough energy to over- 
come the barrier on the left but not on the right. The particle is thus free to 
move to negative Infinity, so Its allowed jt region Is between — ■= and jr,. 

Here, the allowed energy values are continuous, not discrete, because the par- 
ticle isn't completely hminrl. The energy eigenvalues are not degenerate — 
that is, no two energy eigenvalues are the same (see Chapter 2 for more on 
eigenvalues }. 

The Schrfrdtnger equation is a second-order differential equation, so it has 
two linearly Independent solutions; however, In this ease, only one of those 
solutions Is physical and doesn't diverge. 



T 
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The wave equation In this case turns out to oscillate for x x._, and to decay 
rapidly four > a*. 

CftJP 2: Energy greater than the higher potential (E > V 2 \ 

If E V . the particle isn't liowid at all and is free to travel from negative infin- 
ity to positive infinity. 

The energy spectrum Is continuous and the wave function turns out to l>e a 
sum of a function moving to the right and one moving to the loft, The energy 
levels of the allowed spectrum are therefore doubly degenerate. 

That's all the overview you need — time to start solving tire Schrodluger 
equation For various different potentials, starting with the easiest of silt inf i- 
1 1 i te square walls. 

Trapping Particles in Infinite 
Square Potential Wells 

Infinite square wells, in which the walls go to infinity, are a favorite in physics 
problems. You explore the quantum physics take on these problems in this 
s Rptiorr, 



Finding a u'at/e -function equation 

Take a look at the infinite square well llsal appears back in Figure 3-1. Here's 
what I bat. square well looks like: 

1 ^ V( r) an. where x < Cl 
w V(a) * 0, where 0 < x<a 
v* V(,r) = *, w h ere > a 

The Schrtkimger equation looks like this in three dimensions: 

7 r ' r| - Aur ( r ) + V |V|y \ r ) = t>l/ ) 

Wtltlng out the Schrbdlngesr equation gives you the following: 

^(&' + sy + a^)*' (, ') + ' r kP('-) = Ev{d 
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You're Interested In only one dimension — x (distance) — In this chapter, so 
the Schrddlnger equation looks like 

in, 



Because V(jc) = 0 inside tire well, tire equation becoin.es 



=£ dl 
2m dx- 



p(*) = Ey{x) 



Anri in problems all this sort, the equation is usually written, as 

£ s ^) + *>(*]=<) 

f| n 

where ft- = (k is the wave number). 

So now you have a second-order differential equation to solve for lire wave 
function ol a particle trapped in an infinite square well. 

You get two Independent solutions because this equation Is a second-order 

differential equation; 

V|M = A sin(j kx) 
iji. f.n = 1:5 casfiSt.x) 

A anti B are constants that are yet to be determined. 

The general solution of *Wft V( y) =0 is the sum of P-{*l aild V- (-'■ ) • 

ctx ' * s * 

l|i(X) ■* A *ln(fer) • R cosQbO 



Determining the energy tenets 

The equation Hi) = A $infft,v) - R cos(ft-r) tells you that you have to use the 
boundary conditiDris to iuid the constants A and. B (I be preceding sect ion 
explains huw to derive the equation). What are the boundary cmiilitions? The 
wave Junction in c s t disappear at the boundaries of an ihliuite square well, so 

I ^o(Ci)=C> 

** v{d) - ft 
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The tact that iy(0) = 0 tells you right away that Li must he zero because 
cosfftj = 1 . Ami the iact that w(a) = iJ tells you that iji(e) = A sintfei) = 0. 
Because slue; is. zero whcis its anjuiiLcnt is a multiple oi s, this means that 

fttt - rtft n. 1,2,3... 

b'ole that although n * 0 Is technically a solution, it yields ijr(O) - 0, so It's not 
a physical solution — the physical solutions begin with n - 1. 

This equation can also be written as 



km— rr-1,2.3... 

it 

And because ft 1 - 2j?jF/ia-, you have the following equation, where n ■ 1,2, 
3, - — those are the allowed energy states- These are quantized states* 
corresponding to the quantum numbers 1 , 2.. 3, and SO on: 



2 mV, _ ft x' 

if if 

P = n Vi vr 

2mfia : 

Note that the first physical state corresponds to .n = 1, which gives you (his 
nc\t equation: 

E = JL2lL 

2ma 1 

This is The lowest physical state that Itie particles ran occupy. .lust for kicks, 
put some numbers into this, assuming that you have an electron, mass 9.11 x 
fiCT : " kilograms, confined to an infinite square well of width ol the order ol 
ffohrmJrus (the average radios ot an electron’s orbit in a hydrogen atom), 
about 10 _|: meters. 



F. = 



— - _ niveg you this energy for the ground state; 
2itki 1 



(] 05 x UP 1 ) 1 (3.14)' 
3(9.11 x ]« ')(|i.i 



= 6.00xl0' w Joules 



J py 
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That's a very small amount, about 10 electron volts (eV — ihe amount of 
energy one electron gains falling through l volt). Even so. It’s already on the 
order of the energy of the ground state of an electron In the ground state of 
a hydrogen atom (13.6 eV), so you can say you're certainly in the right quan- 
tum physics ballpark now, 



Normalizing the ttfatfe function 

Okay , ytm have this fur the wave equalluji lor a particle in art infinite square 
well: 



^(rr) = A Sir.' ns * j 

Tlie wave function is a sine wave, going to zero at x - 0 and x - o. You can see 
the first two wave functions plotted in Figure >'V 



Fbguw 3'3: 
W (.m 
funcliDns 
in a square 
well. 
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Normalizing the wave function lets you solve lor the unknown constant A In 
a mrmaiizeef lunctlon, the probability of finding the particle between x and 
dx, ly(*i -dx, adds up to I when you Integrate over the whole square well, 
x - 0 to x - a: 

l = ] V-'U:)| J l£v 

I 



Substituting for u(x) gives you the following: 



L =| A ] Sin | ^ jif 



Here's what the integral in this equation equals: 

j sbr | - 1 '" - jtir - % 

So from Hie previous equation, 1 = A J <l |. Solve lor A: 

A= (f)‘ 

Therefore, here’s the normalized wave equation with Elm value of A plugged in: 
And that’s the normalized wave inaction for a particle in an Infinite square well 



Adding time dependence 
to Wave functions 

Now- how about seeing Jiuw l lie wave Junction lor a particle in an inEliiitu 
square well evolves with time? The Schrfitlinger equation looks like this: 

^^fr)e v l>y{r)=Ety(r} 



Tiqhi 



materia 
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You can also write the SchrSdinger equation this way, where H Is the 
Hermitbara Hamiltonian operator: 

IMU - F.ij-fr) 

That’s actually the tme^ndepenttent SchrodinRer equation, The time-ctepen&nt 
Schrudinger equation looks like this: 



A yrfr\ f )=HtyjV, f) 

Combining the preceding three equations gives you the following, which is 
another lorm of the thne-depe^dent Schrddinger equation: 

&!! { y(r. /)= \v(r. tj+Vfr, t) 



And because you’ re dealing with only one dimension, x, this equation 
becomes 



& V^*. *) = 2 *. £ K-r. 0 + v [ x < s ) 

This is simpler than it looks, however, because the potential doesn't change 
with lime. In fact* because E is constant, you can rewrite the equation as 



iy{x, f)=E^[^.r) 

That cquarion makes tile- a Lot simpler — its easy to solve Hip time-dependent 
Schrtidiiiger equation it you’ re d-ealing with a constant potential. In this case, 
the solution is 



t^-v. /]= i^.rje 



* lt A 




Neat. When the potential doesn’t vary with time, the solution to the lime- 
dependent Schr&dinger equation simply becomes y(jf), the spatial part, nuilil- 
plied by the time-dependent part. 
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So when you add in the time-dependent part to the time-independent wave 
function, you get the time-dependent wave function, which looks like this; 




rr = ] 2. It ... 



The energy of the nth quantum state is 



rtVff 3 



I, 2, 3 



2mr 

Therefore, the result is 
where exp \ x ] = e " . 



jj.ir.v 

a 






Shifting to symmetric square 
&eft potentials 

The standard inEinite square well Errcrks like! this; 

I t* V(j,-) = x. where a- -■ t) 
t* V(x) « 0, where 0 s x s! a 
** V(x) * o:. where x > a 



But what If you want to shift things so that the square well is symmetric around 
the origin instead? That Is, you move the square well so that it extends from 
- J /± In ■'/,? Here’s what She new infinite square well looks like in this ease: 

1 ^ V(jr) = x, vi here x ■- ->■/: 
v* V(jr) = 0, where — k <x±-ji 
f* V(i-} - t „ where x w fz 



You can translate from this new square well to the old one by adding "h Iojc, 
which means that you can write the wave function for the new square well ‘m 
this equation like the following; 
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sm 



'tH)) 



Doing a little trig gives you the following equations: 



¥{x' r 


i ■ 

— COS 

{a } 


nx\ x | 
0 


ji- 1. 3,5... 






J 






L?) ,sln 


flK \ X j 

a 


ft *2.4,6... 



So as you can see, the result is a m5>: of sines and cosines. The bound states 
are these, In Increasing quantum order: 




And so cm. 

Note that the cosines are symmetric around tin- origin: c(.i.) = i+r f— .cj . 'The 
sines are aid i-symiiivirie: -v(jc) = lyf-vl. 



Limited Potential: Taking a Look at 
Particles and Potential Steps 

Truly infinite potentiate {which I discuss in the previous sections) are hard 
to come by. in this section, you look at some real-world examples, where 
t lu* potential is set to some finite V,. not infinity. For example, take a look at 
the situation in Figure 3-4 . "I in- re. a particle is traveling toward a potential 
step. Currently, the particle is in a region where V’ = 0, but it’ll soon be in Lbe 
region V - V M . 
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J 

e - - - - 

v * 


V 

L 






Figure J-t: 




A p ore rii-il 


X 


step, E :=■ V,. 





o 



There are iwo case a to look at here In terms of E, the energy ol the particle: 

* F. > V ( : Classically, when F, > V', r you expect the particle to he able to con- 
tinue on to the region x ■- 0, 

v* 1L< V 4 : Whtn E « l' , you'd expert the particle to bounce hack arid not be 
ab§e to gel to the region x -■ 0 at all. 

Iji this section, you start by taking a look at the case where the particle’s 
energy, E, is greater than the potential V M , as shown in Figure 3-4; then you 
take a look at the case where E < V M „ 



Assuming the particle has 
ptent y of emtqif 

Start with the case where the particle's energy, E, is greater than the poten- 
tial V . FroiJ] a quantum physics point of view, here's what Hie Schrodulgcr 
equation would look like: 

v* Fur the region * < 0; ^ , 1 ( jr ) +■ x | ■ 0 




'r- 
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r/ J, fcy , _ , 

^ For the region x > 0: ** [,? | + A,V, I v | = 0 

In this equation, k,' ■ 

In other words, ft is going lo vary by region, as you see In Figure 3-5, 



N«r« 3'S: 
Tlio value el 
.V by region, 
whera 





Treating the first equation as a second-order differential equation, you can 
see that the most general solution Is the following; 

ip , (y) ~ Ae* ■ » Be"* 1 ', wbere x 0 

And lor the regloEi jr > 0, salving the second equation jjjvKS you I his: 

n+i .r.ii j = ^ - Lfe where x ■- Cl 

Note that r* Y represents plane waves traveling in the ~x direction, and. er m 
represents plane waves traveling in the ~.v direction. 
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What this solution means is that waves can hit the potential step from the 
loft and. be either transmitted! or reflected. Given flint wav of looking at the 
problem, yon may note that the wave can he reflected only going to the right, 
m>r to the left, so n must equal zero, That makes the wave equation become 
the following; 

I v* Where .r <P: ip, (rj = Ae'*' ' + fie "h ' 
v* Where * > y.,(x) = C<?" Vf 

The term Atr" 1 1 represents the incident wave. Be"*""* is. the reflected wave, anti 
Ce^is tlie transmitted wave. 

Calculating the probability ef reffecthit vr traitsmissiffn 

You earl calculate the probability that the particle will he reflected or trans- 
mitted through the potential step by calculating the reflection and frunswi/s- 
Sd'orr coefficient if J, is the reflected current density, I, is the incident cttrrenL 
density, and I, is the transmitted current densiLy, then Ft, the reflection coef- 
ficient is 




T, tho transmission coefficient, is 




You now have to calculate J r . J,. and J r Actually, that's not so hard — start 
with J,. Because the incident part of the wave Is v,(A) ^ Ae' l, : v , the Incident 
current density Is 



J, 



2 m 






~dx~ 



- 'A ( '• ; 



£^!a 0 

itx 



And ibis just equals 




. J r and J, work in. the same way: 
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So you have this for the reflection coefficient: 




T, Uie Irauscnissiun coefficient, is- 



T = 



L 

J., 




finding A, & t and C 

So tiow du yuu figure Out the constants A, B, and C? You do that us you figure 
out ituj coefficients with the infinite square well potential — with boundary 
conditions (see the earlier section “Trapping Particles in Infinite Square Well 
Potentials"). I lowever, here, you can't necessarily say that goes to zero, 
because the potential is no longer infinite. Instead, the boundary conditions 
are that y(jc) and dv(x)?dx are continuous across the potential step's bound- 
ary, In other words, 

I t^ ip | (ff)=Hr E (hl 



You know the following: 

I * Where X < «; if , (.rj - At?'*** + Be*** 
r" Where r>fl: VvOO ■ Ce'* T 

Therefore, plugging these two equations into y ,(()) - y,(0) gives ymi A + H =f!. 

And plugging them into : ( 0 ) ■ [til gives you 

CtV ' flV 

- kf 

Solving for B- in terms of A gives you this result: 
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Solving for C In terms of A gives you 



C- 



2k, 

k + * 



A 



You can then calculate A from the jiormalLzatitiii condition of the wave 
function: 



But you don't actually need A, because it drops out oi the ratios for the 
rei lection and transmission coefficients, Rand T. In particular. 





Therefore.. 



"(m*J 

T= 

i>, t^;r 

That's an interesting result, and it disagrees with classical physios, which 
says that there should be no particle reflect icm at all. .As yon can see, if A, * 
k . . then their will indeed he particle reflection. 

Note that as fr t goes to A*, Rgoes to 0 and T goes to 1. which is what you'd 
expect. 

So already you have a result that differs from the classical — the particle cats 
be reflected at the potential step. That's the wave-like behavior of the par- 
ticle coming into play" again. 
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Assuming the particle doesn't 
haPe enough energg 

Okay, now Iry the case where t ■- V , when there's a potential step, as shown 
in Injure 3-6. In this case, rtic particle doesn't have enough energy to make il 
into l tic region x -■ 0. according to classical physics. See what quantum phys- 
ics has tu say about it. 



1 1 

i 



F 



o 



FtflUW 3 6: 
A |Uin:fi1inl 
1 ■;! |3 , E -i V 



Tackle the region x -■ « first. There, the SchrSdiiiger equation would look 
like this: 



dx 1 






where k 



t _ 'ImEf 
1 " /&* ■ 



You know the solution lo this Irom die previous discussion] on potential steps 
(see 'Limited Potential: Taking a Look at Particles and Potential Steps'): 



! l- 



^(jr) = Ae“ + Be ** 
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Okay, but what about the region x > 0? Thai's a different story. Here's the 
Sc h rodln ger equation: 






(where *>0} 



wl, e r«' = ME-V.,l/ 

But hang on; E- V' M la less than Kero, which would make k imaginary, which is 
impossible physically. So change the sign in the SchrBdiagcr equation from 
plus to minus: 



dx 1 



(*)-*>, [*) = o 



x >D 



And use the tallowing for *. (note that this is positive if E < V„): 
2 m [V r -E] 

•• ” 



Okay, so now ynu have to solve the rlLFFcmiti.il 



^7'!*) *>..(*)- o 

solutions: 



(where jt > 0}. There are twn linearly independent 



1 ^ h Cr*^ 
v* i^-c) = rW 



Anri the general solution to 



■v- 






| .r | — k u/ , i .v j ■ 0 



(where jr > 0) is 



f/,|> ) ■ O M i EV a Ct 

However, wave functions must be Unite everywhere, and the second term 
is dearly not Finite as r goes to infinity, so D must equal zero (note that IF x 
goes In negalive inlinily, the Fust term also diverges, hut because the poten- 
tial step is limited lo x > (l, that isn't a problem 1. Therefore, here's she solu- 
tion For A" > II: 



¥,{x) ?.Ce' kA X -■ 0 
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So your wave functions for the two regions are 

i|> ,(.*•) - * Re* 1 "'' X’-d 

iy„i> j - Ce-^ jr>(l 

Fatting this in terms, of the incident, reflected, and transmitted wave func- 
tions. litt-i"), and hj i i (.v'i. you liavc the Julluwing: 

I V* NI|X) = rV“ A ' V 
V* 4i r i» ■= Be"*' 1 ' 
v* vix) * Ce^ 

finding ttattsmtesivn and tefiectwn coefficients 

Now you can figure out the reflection and transmission coefficients, R and 
T (as you do for the ease £ ■ V . in tho earlier section “Assuming the particle 
has plenty of energy"): 



Actually, this is very easy in this case: take a look at J,: 



f M x ) 












But because ^X\) - Ctr*--', hj£0 Is completely real, which means that In this 
case, the following Is true. 



J. 




dr 






* J 



A nd this equation, of course, Is equal to aero. 

So j r = 0: therefore, T = 0. BE T = ft, then R in List aqiLftl l. That means that you 
have a complete reflection, just as in the classical solution. 



The name rtf saiutwtt; Finding a p&rticfe in x > 0 

Despite the complete reflection, there’s a difference between the mathemati- 
cal aud classical solution: There actually is a nonzero chance of finding the 
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particle In the region x > 0. To see that, take a look at the probability density 
for x > 0, which is 

poo - 1 m ) - 

Plugging in for the wave function W>) gives v<n.i 
P(.r) = I i|r,(.v) I * = IC -e-- 1 *' 

You can use llie continuity conditions to solve for C in terms of A: 






Using the rontinnity conditions gives you the following: 



Ak , 1 A ; e a - L 






This does fall quickly to zero as x gels large, but near x - 0, it has a nonzero 
value. 



You can see what the probability density looks like lor the F « V M case of a 
potential step in Figure 3-7. 
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Okay, you’ve taken care of infinite square wells and potential steps. Now 
what about the case where the potential step doesn’t extend out to infinity 
but Is itself bounded; That brings you to potential frame rs, which S discuss in 
the nevt section. 

Hitting the Watt: Particles 
and Potential Barriers 

What if the particle could work its way through a potential step — that is. the 
step was ol limited extent? Then you’d have a potential barrier,, which is set 

up aomettilrii.r ifc- thist 

I v* V(.t ) = o, where x < n 
v* V(y) = V . where P s x £ a 
v* Vi.v) = 0, wherrs x ■> a 

You can see witai this potential looks like irs Figure 3-8. 



V 
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In solving the Schrddlnger equation lor a potential barrier, von have to con- 
sider two cases, corresponding to whether the particle has more or less 
energy than the potential harrier. In other words, if E Is the energy ol the 
Incident particle, the two cases to consider are F > V„ and F. < V„, This section 
Starts with F. V., 



Getting through potential 
harriers token E>V 0 



In the case where E the particle lias enough energy lo pass through the 
potential barrier and end up in ihex - a region. This is what the Schrddinger 
equation looks like: 



b* For the region x -r. (I: f x ) + V ( x j = b 



where to 2 = 



' tjY ii ■ i 

V Far the region OlrS or — 7 I x | + x) = It 



where fe , 3 = 



_2m(E~V n ) / 



V 



v* Far the nglon * > or ^£7 [ x \ + k, ty , ( x } 0 



where k. 

The solutions fo r ijj L ( jr), v,(x), and ijj,(x)are the following: 



^ Where x < #: y , i» - Ae'*' x * 

v* Where 0 S x £ m iy,{.v) = C.f*?' t IlKr'E" 

v* Where x > ui v (x,i = !•> h-' t I'e-' 

In fact, because there's no leftward traveling wave Ln tliex> fa region, 
]■' = E), so o (xi = Ee ,L i v . 



Copyrighted material 




80 



Part II: Bound and Undeiemiinfid: Handling Pari in Ins in Bound States 



So h.ow do you determine A, B, C. D, and E? You use the continuity condi- 
tions, which work out here to he the following: 






c/^f, 

dsr 








Okay, from these equations, you gel the following: 

P* A ■ R I C l 1} 
v* J*,(A - R ) = J*,(C - T>> 
v* Cfl ,A - rt t I3(r ll i u = Yjs ltr \ a 
v* ikCtr*-- - rt^Def 11 '-" = 

So putting: all of these equations together, you get this for the coefficient E in 
terms Of A: 



E - 4k,k,Ac * [ 4k, it, cos ( ft.tr ) - 2 1 ( it . 1 + k . ) sin f k,u } J ' 
Wow. So what’s tlie transmission coefficient, T? Well. T is 




And this works out lo he 



T 




k A 



sin/ 



i 



(M) 



Whew! Note that ns A. goes to ft,, T goes to l, which Is what you'd expect. 
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So how about R, the reflection coed Sclent? I'll spare you the algebra, here's 
what R equals: 




Von can see what the I- > V, probability density.. I v(x') I J . looks like for the 
potential barrier In Figure 3-0, 



Figure 3-S: 

a patn-nlial 
barrier 

t-v, 




That, completes The potential barrier when E ;■ V . 

Getting through potential harriers, 
eOen When E < (/ 

What happens If the particle doesn't !m>e as much, energy as the potential 
of the barrier? In other words, you're now facing the situation you see In 
Figure 3- 111 



'TIQrT 



srial 
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v 

A 



E 



c- v 



figure 3-ID; 
A potmilial 
barrier 
E -i V .. 



0 



il 



X 



Now the Schrdtlinger equation looks tike this: 
** For the region x <■ 0; rp L (jc) - A^ 1 ' + 



d tif 



v* For the region 0 ^ * < «: 1 x I + & V j(jf) = 0 



where k = 



_Zm(F.-V..)/ 

/v- 



But now E - V .is less than 0, which would make k imaginary. Arid that’s 
Impossible physically. So change the sign in the Schr6dinger equation 
from plus to minus: 



And use IrhLs lur k,\ kj = ^ j/ 

■ ■ ■' 

For the irf^jflon _v ■ ^r(*H>W=o 

where if . 1 = ■- 
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All tills means that Hie solutions lor v/i'}, i|r. 1 (jr)„ anti, ny./.v) am the following: 

v* Where .v < ty , i.vj = + fie f 

e 1 Whftr 0 5 .t£ a; w .(-0 - Ct jl, - V * Dtr l '- r 
Where x>tr. ■ Ee'* |J + Fe _ *-- V 

In fact, there's no leftward (raveling wave In the region * > a\ F ■ (L so 
F ^ lVfT - 

This situation is similar to the rase where t > V . exrep t for rhe region 
I) £ x 'l a. The wave function oscillates in Hie regions where it lias positive 
energy, x < LI and x ■■ a. but Is a decaying exponential in the: region 0 £ x £ a. 

You can sec what the probability density, ijj(.t) c \ looks like in Figure 3-1 1. 



hl'WI 1 




Finding the reflection md transmission coefficients 

How about the reflection and transmission coefficients, R and T? Here's what 
they equal: 





T = 
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As you may expect, you use the continuity conditions lo determine A, R, and E: 

^ uv'i)> » V.IMJ 
dys, 

^ W.Xa) = 









dx 



A lair bit of algebra and trig is involved in solving lur R and T' htircs wlilt R 
and T turn Out "0 be: 



R = 



* L±*1 

ft,*. 



sin/i 



M 



A cos h 



M 



(k.’+kr 






sin/i (A,o| 



T = 4 



■1 cos/i | k .a \ I 



V-^/ v 

" k,k, 



sisi /2 |6.tr j 



Respite the equations complexity, it’s amazing (hat the expression for T eau 
be nonzero. Classically, particles can't enter the forbidden zone it £ x i a 
because E « V. , where V (l is the potential Ln that region; they just don’t have 
enough energy to make it into that area. 

Tunneling th rvugh 

Quantum mechanically, the phenomenon where panicles can get through 
regions that they're classically forbidden to enter is called Htnnttimg. 
Tunneling is possible because in quantum mechanics, particles show wave 
properties. 

J'unnnlmjj is one of tilt! most exciting results oE quantum physics — it means 
that particles can actually gel through classically forbidden regions because of 
the spread in their wave functions, fhis is. of course, a microscopic effect — 
don't try lo walk through any closed doors — but si's a significant one. Among 
other effects, tunneling makes transistors and integrated circuits possible. 

You can calculate the transmission coefficient, which fells yon the prob- 
ability that a particle gets through, given a certain incident intensity, when 
tunneling is involved. Ruing so is relatively easy in the preceding section 
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because the barrier that the particle has to get through is a square barrier. 
But In general, calculating the transmission coefficient Isn't so easy. Read on, 

Getting the transmission ufith the WKB approximation 

The way you generally calculate the tram mission coefficient is to break up 
the potential you're working with into a succession el square barriers anti to 
sum them. That's called the Wenael-KrumernSrilbum (WKB) approximation — 
treating a general potential V(„y}. as a sum of square potential barriers. 

The result of the WKR approximation ig, that the transmission coefficient for 
an arbitrary potential. V{x ), ior a particle oE mass m and energy F. is given by 
this expression (that Is. as long as V{jr) is a smooth, slowly varying function): 



T - exp 




tlx 



So now you can amaze your friends bv calculating the probability that a par- 
ticle wili tunnel through an arbitrary potential. It's the stuff science fiction is 
made of — well, on the microscopic scale, anyway. 



Particles Unbound: Soloing the 
Schrodinger Equation for Free Particles 

What a haul particles outside any square well — that is. free particles? There 
are plenty of particles that act freely in Hie universe, and quantum physics 
has something to say about them. 

Here's the Schrddinger equation: 

What St the particle were a free particle, with V(.y) - 0? In that case, you'd 
have the following equation: 
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And you can rewrite this as 

^_V ^ (jr ) + * V{ -V ) = d 

where the wave number, k h k f ■ \ 

ft 

You can write the genural solution to this SchrfSdinger equation as 
41 (jc) = Ai? 1 * 1 ' - Be -1 * 1 ' 

If you add time-dependence to the equation, yon get this time-dependent 
wave fiinc t ion ; 



v 1 '.] r, I ) = Aevp ■ + R expj -| ik - r - | 



That's a solution to the Schrttdinger equation, but it turns out to be un physi- 
cal. To see this, note that lor either term In the equation, you can’t normalize 
the probability density, 1 [see the earlier section titled "KormalmnR 
the wave function" lor more on normalizing); 

I Y'lVi I - = I A I - nr I R 1 2 

What's goiiii) on here? The probability density for the position ol the particle 
is uniform throughout el! r: 1 In oilier words, you can’t pin down the particle 
at ah. 

This is a result ol the form of the time-dependent wave function, which uses 
an exact value for the wave number, k — and p ■ fik and E - b&.Qm. So what 
that equation says is that you know R and p exactly. And if you know p and R 
exactly, that causes a large uncertainty in x and t — in fact, x and / are com- 
pletely uncertain. That doesn’t r nrres pnn rl to physical reality. 

tor that matter, tire wave function i|j(.r), as it stands, isn't something you 
ran normalize. Trying to normalize the first term, for example, gives you this 
integral; 






.Anri Inr the- first term of ipf.v, 1), this is 



J y( x )y'(x)dx= a| j J Lie- 
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And the (amt is true dE the second terras Ln w(t. t). 

So whal do yon do here Id gel a physical particle? The next section explains. 



Getting a physical particle 
(Pith a ufaPe packet 




if you have a number of solutions to the Sehrodiiriger equation* any linear 
combination of those soJaations La also a solution. So that's the key to ijetling 
a physical particle: You add various wave functions together so lhal you yet a 
tcave packet, which ts a cuJlcctiori oi wave functions ol the iorm c ,, ‘" |:,V '' 11 sucti 
tJiat ihc wave functions interfere constructively at one Location and interfere 
destructively (go to zero) at alt other locations: 



4r.fj X*,-"' ' 1 



This ts usually written as a continuous Integral: 






■'*' IJ ..I 



(*0 : - 



(IP 



What is <y(k. t)'l It's the amplitude of each component wave function, and you 
can find /] from the Fourier transform of the equation; 



P(*. /)= — 1 4r J V'(* »)» 
( 2 ' T )“ " 



Because k = p- ft, you can also write the wave packet equations like this, in 
terms of p, not k: 



') = 




' ' ( 2 *)=- 



Copyrighted material 




Part II: Bound and Unde term in fid: Handling Parti des in Bound States 



Well, you may be asking yourself Just what's going on here. It loots like 
fO. O Is defined. In terms of o(p, f). but ofa t) Is defined In terms ol v(x f)= 
That looks pretty circular. 

The answer Is that the two previous equations aren't definitions of y(r. f) 
or $0?. fl\ they're just equations relating the two. You're free to choose your 
lswu wave [racket shape yourseli — lor example, you lELay specify the shape 

+“■ 

ml 6(p, 0- and lyf x, s ',= — ! — r je A dk would let you fimJ iy(_r, t). 
\ 2 x) 2 - 



Going through a Gaussian example 

Here's art example iu which you get concrete. Selecting ah actual wave packet 
shape, Choose a so-called Gaussian wave packet, which you can see in 
Figure 3- 1 2 — localized in one place, zero in the others. 




The amplitude •«*) you may choose for this wave packet Is 



\v j &| = a exp 



/■} 



You start by normalizing <H*) to determine what A is. Here's how that works: 



=JH*i 



dk 



Substituting in gives you this equation: 
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Doing the Integral (that means looking It up In math tables) gives you the 
following: 



1 A 



2 */ 
L / 






■/J- 



TherdOre. A - 
So here's vDiir wave hjiictiosi: 






[^) 



-o' \k~k-)/ 



e*'dk 



Tli is little gem ol an integral can he evaluated to give you the following: 



*MynA'*4'%Y‘ 

So that’s l h e wave fund ion for this. Gaussian wave packet (AW.- The 
e* p j is ihe Gaussian part that gives the wave packet the distinctive 
shape that yw see in Figure 3-12) — and it‘s already normalized- 

Now you can use this wave? packet Function to determine tins probability that 
the particle will be hi, say-, the- region 0 £ x £■'/.. The probability is 




In this case, the Integral 3s 




And this works out to he 




So the probability that the particle will be In the region 0<jc£ is is '/j. CooLI 



Tiari; 
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Chapter 4 

Back and Forth with 
Harmonic Oscillators 



h This Chapter 

Hamiltonians: Looking at total energy 

Solving for energy stares wiTti creation amt annihilation operators 
Understanding the matrix version nl harmonic oscillator operators 
Writing computer code to solve the Schr&dlnger equation 



M^urtmtruc mcilkuofs are physics setups with periodic motion, sue has 
¥ I things bouncing on springs or tlck-tocklng on pendulums. You're prob- 
ably already familiar with harmonic oscillator problems In the macroscopic 
arena, but: now you're going microscopic. There are many, many physical 
cases that can be approximated by harmonic oscillators, such ns atoms in a 
cits Eat structure. 



In this chapter, you sec both exacl solutions to harmonic oscillator problems 
as well as computational methods for solving them. Knowing how to solve 
tile Sdirodiilger equation, using computers is a usciul skill lor any quantum 
physics expert, 



Grappling with the Harmonic 
Oscillator Hamiltonians 



Okay, lime to star? talking HaiiLiLtuinans (and I'm nut referring to Fairs of the 
U.S. Founding Father Alexander Hamilton). The Hamiltonian will Jet you Mud 
the energy levels of a system.. 




n 
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Going classical with harmonic oscillation 

In classical terms, tire- force on an abject In harmonic oscillation is the follow- 
ing (this is Hooke's law): 



J ; = -kx 



Iri tlsit equation, k is rtic spring constants measured in Newtons/ meter, and or 
Ls displacement. The key paint here is ihat the restoring force ■uji whatever is 
in harmonic iimU ion is prupOrtim)al to its displacement. In other words, the 
farther you St retd) U spring, the harder it’ll poll back. 

Because F *■ ma t where m Is the mass of the particle In harmonic motion and 
a is Its instantaneous acceleration, yon can substitute lor F and write this 

equation as 
mat * fev - f j 

Here a the equal iui) tor instantaneous acceleration, where x is displacement 
and r ii time: 




So substituting for o. you can rewrite the force equation as 

ma r fee = . - + fee a 0 

rfr 

Di viding by the mass of the particle gives you the following: 



(i) + kx = £ 
di 1 m m 

If you take */„ - m (where m is the angular Frequency), this becomes 

" x = U 

tit 

You can solve this equation lor x, where \ and Hare constants; 



-VV 
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X = A s Inosf i B COS <<’4 
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Therefore* the solution Is an oscillating one because It in wolves sines and 
cosines, which represent periodic waveforms. 



Understanding total energg 
in quantum oscillation 

Now look at harmonic uscillaturii in quantum physics terms. The fftnrrflltwid'fln 
(H) is the sum of kinetic and potential energies — tlm total energy of file 
system: 

1 1 - KE * PE 



For a harmonic oscillator, here’s what these energies are equal to; 



o* The kinetic energy at any one moment is the Following, where p is the 
particle's momentum and m is its mess: 



KE = 



iL. 

2rrr 



** The particle’s potential energy Is equal to I he Following, where ft L? the 
spring constant and x is displacement; 

PE l hr- 



{Able; The ft cancels out because or = ) 

Therefore, in quantum physics terms, you can write the Hamiltonian as H = 
KE ■+ PL, or 



H *XL + lirmr'X t 
2m 2 



where P and X are the momentum and position operators. 

You can apply the Hamiltonian operator to various eigenstates (see Chapter 
2 for more oil eigenstates), i w>. of the harmonic oscillator to get the total 
energy, E, ol those eigenstates; 

I I ',v > | if/ > + " mm ; X ' ur > E v > 



opyrig 
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The problem now becomes one of rinding the eigenstates and eigenvalues. 
However., this doesn’t turn out to be an easy task. Unlike the potentials 
V{jr) covered in Chapter I, V(jt] for a harmonic oscillator is more complex, 
rfepen ding as it does on x' . 

ho you have to he clever. The way you solve harmonic oscillator problems in 
[LLianrum physics is with operator algebra — drat is, you introduce a new set 
of operators. And they’re coming up now. 

Creation and Annihilation: Introducing 
the Harmonic Oscillator Operators 

Creation ^itd annihilation may sound like big make-or-break the universe kinds 
□f ideas, but they play a starring role hi the quantum world when you're work- 
ing with harmonic Oscillators. You use the creation and annihilation operators 
to solve harmonic oscillator problems because doing so is a clever way id 
handling the tougher Ham lit on Ian equation (see the preceding section). Here’s 
what these two new operators do: 

t*" Creation operator: The creation operator raises the energy level of an 
eigenstate by one level, so IS the harmonic oscillator Is in the fourth 
energy level, the creation operator raises M to the fifth level. 

& Annihilation operator: The annihilation operator does Itic reverse, low- 
ering eigenstates out level. 

These operators make it easier to solve for the energy spectrum without 
a Lot Of work solving for the actual eigenstates. Ill Ollier words, you Cap 
understand the whole energy spectrum by looking at the energy difference 
between eigenstates. 

Mind your p*$ md c{*s: Getting the 
energy state equations 

Here's how people usually solve for the energy spectrum. First, you intro- 
duce two new operators, p and q, which are dimension less; they relate to the 
P (momentum) operator this way: 
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You use these two new operators., p and q r as the basis of the annihilation 
operator, a. amt the creation operator, o ' : 

^ a ' = ^~ tp) 

Now you can wrim t Jits harmonic oscillator Hamiltonian liku this, in terms of a 
and i?: 



H=J | J w[( f , d + i] 

As for creating new operators here, The quantum physicists wen! crazy, even 
giving a name fa rj a: the .Vac number operator. Sa here's how you can write 
tJio Hamiltonian: 

The JS operator returns the number ai the energy level cl rtie harm runic oscil- 
Lalor. if you denote tlie eigenstates of H as !n>, you get this, where rr is Hid 
number of the nth state: 

N ii (T> = n I n* 

Because H nMF * '/■), and because H l rr - * E, n hi-*, then by comparing the 
previous, two equations, you have 



E. =| n + i'/rfli n=«,1.2... 

Amazingly, that gives you the energy eigenvalues of the nth state of a quan- 
tum mechanical harmonic oscillator. So here are the energy states: 

* The ground .state energy corresponds to n - 0: 

F, = l/ f( y 
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e* Tlie lirst excited State is 

E.-jAnJ 

The second excited state has an energy of 

e,=|A™ 

And so on, That Is, the energy levels are discrete and nowtegenemte (not 
shared by any two states). Thus, the energy spectrum is made up i>f equidis- 
tant bands. 



Finding the Eigenstates 

When you have the eigenstates (see Chapter 2 to find out all about eigen- 
states), you can determine the allowable states of a system and the relative 
probability that the system will he in any o! those states. 

The commutator ot operators A, ES is | A, E5] = AK - HA, so note that the com- 
mutator ol a and a* is the following: 

[«,u'J = --[y + fp.y-«p] 

This is equal to the following: 

[ fl - °' ] = *H q + * p ' q ~ H = “4 ■?- P I 

This equation breaks down to [o, a J = 1. And putting together this equation 
with N+ -Lj, you get |o, H = Aoki ansi | ta 8 , 3-1 ] = -ham' . 
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Finding the energy &( n\n> 

Okay, w i t] i the commutator relations, you're ready to go. The Sirs l Question is 
M the energy of slate in* is E,„ wfial is the energy ol the state til rr>? You can 
write the energy of a I this way: 

H(/j n •> 

= (rtlt - li uul ) I n> 

= (E n -}]to}(n!n:>) 



So am ■ is also an eigenstate of the harmonic oscillator, with energy E, - tiro, 
not E„ That's why a is called the annihilation or lowering operator: ft lowers 
the energy level of a harmonic oscillator eigenstate by one level. 

Finding the energy &{ A V b> 

So what's the energy level of ti !«=■? You can write that can like lliis: 



H|o rt>) 

a ( a *1-1 I fHwa | n > 

= (E n + ft*j)(<r , |n>) 

All this means that d I n> is an eigenstate ut the harmonic oscillator, with 
energy E„ + r?w, not Just — that Is, the a operator raises the energy level of 
an eigenstate of the harmonic oscillator by one level. 

Using a and a f directly 

If you've been following along Irom the preceding section, you know that 
H(ol m J = (fc„ - o>J and H(a T l it>J = + ficoXrr lo^). You can derive The 

foltowing Irom the these equations: 

I V a n ■ = C n- [■' 

V a' rr> = [?l n * 1^ 
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C anti D ate positive constants, but what do they equal? The states I n - 1> 
and in h have to be normalized, which means that *n - 1 1 rr - 1 > ■ 

<tt i l n • l> - ]. So take a look at the quantity using the C operator: 

^WXcrlnv'j-C^rr-llfl-K 

And because I n - 1» is [LuriiLaJizcd, *..n - 1 n - 1* = 3: 

(■=nltf f XLrl rr*) = C ; 

• ri : rr rr I a > - C 

Bui you also know lhaL d'o - N*. the energy level operator, so you get the fol- 
lowing equation: 



• fl I N I ii> - if 

Nln> = n\n>* whore n Is the energy level, so 
n*n l n* - C" 



However. ^rrln> = I. SO 

a - C- 
n W -C 

This finally tells you, from o \n * - Cl n - 1>, that 
a ! n •> = n In- I -■ 

That's mol — now you know flow to use the Lowering; operator, a. an eigert- 

states o! the harmonic oscillator. 

What about the raising operator, o'? Following the same course of reasoning 
you take with the a operator, you can show the following'.. 

aW - u, , |) lfl \n*. i> 

So at this point, you know what the energy eigenvalues are and how the rais- 
ing and Lowering operators aEEer.l the harmonic oscillator eigenstates. You've 
made quite a list of progress, using the a and a* operators instead oi trying to 
solve the Schrodinger equation. 
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Finding the harmonic oscillator 
energy eigenstates 

The charm ol using the operators a and#' Ls that given the ground stare, H)», 
those operators let you fitici all successive energy states. It you want to lisid 
art excited state ol a harmonic oscillator, you can start ivilti tlie ground state, 
! Lt>, and apply the raising operator, a . For example, you can do ttd*: 



|i>au'|0>a|l> 



And so on. In general, you have this relation: 



War ft: tig in position space 

Okay, h J 0 > is fine as tar as it goes — but just whist is i0>? Can't 

v'^l 

you Ret a spatial eigenstate rtf this eigenvector? Something like ipnt-vt, not 
just 10*? Yes, you can. in other words, you want to :bd Ho ymi 

need the representations of a and a in position space. 



The p operator is delta ed as 
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Because P= -Hi ^ , you can write 



. _ -ih d 
P ' ' Tldx 



\ mfit u | 



And writing jit,, = this becomes 



D = — — = — i x 

P - ■ dx dx 



i m/rflf l 



Okay, what about the a operator? You know that 



y£ 



Aj]d that 






Therefore. 



°=i(t +I 4) 



Vuu can also writer ihrs equation as 



u 



xj> 



(*«■■£) 



Okay, so that’s a in the position representation. What's o - That turns out to 
he this; 
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Now’s the El me to be clever. You want to solve for 1 0> In the position space, 
or Here‘s the clever part — when woe use the lowering operator, a, 
on S FI =• , you have to get a because there's no lower state than the ground 
state, so All): = if. And applying the <rl bra gives you ode 10 - = II. 

That's clever because its going to give you a homogeneniui differential equa- 
tion (that is, one that equals, zero). First, you substitute lord: 



■; .v 



■ ■=: x 



o|0 v- € 



0 >= 0 



xj> 

Then you use < .v I) = ^ ,1 x \ : 



■* 1 ^ (*)+-* 



. dy t (rl 



(lx 



m 0 



Multiplying both sidles by x.& gives you the following 



X- W r I a" I — X,* 



^M =n 

dx 



dyjx} _ -.my \x\ 
dx x. ' 

The solution to this compact differential equation Is 

That's a gaussiah function. so the ground state of a quaiatusu mechanical lusr- 
mojiic oscillator is a gaussian curve, as you see in Figure 4-1. 



right© 



*ria 
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Vfl hi 



Figure 4-1; 
Thu ! : iii; u ! id 
state :: J a 
quantum 
mechanical 
bafmasiie 
oscillator. 




Finding the vtafe function of the ground state 

As a gauss ian curve, Lhe ground state of a quantum oscillator is ^ h i>) - A 
exp(-.r‘/2uv). I low can you Hgurt out A? Wave (unctions must he normalized, 
so- the following has to be time: 

' = ]k..(*)> 



Substituting forw IV) gives yoL] this [test equal mu: 



“I- 



dx 



= A j exp| ~ x 



You can. evaluate this integral to be 



l=A J jexpj' ., x . Jctor = A-fl- x t 



Therefore. 
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This means that the wave function for the ground state of a quantum 
mechanical harmonic oscillator Is 



^.!>l = ■ 



*V'i 






Cool. Now you’ve got an esact wave function. 

A tittle excitement; Fit: ditty the first pj pelted state 

Okay, (he preceding section shows you what looks like. What about the 
first excited state, hi, O')- Well, as you know, iihO) - vj l l > and 1 1 = o T iOv, so 

iji 1 1 \ ) = ••-j-ifl' i n • 

And you know Hint s' is the Following: 



ti 




JC 




Therefore, i|i i = Oi‘ o I LN becomes 



jc a I'l > = 




I 

y lV 



<A-(x-jr r z — I (f> 
t. dr Jj 



And because ^..Or) - «x Q\ you gel the following equation: 
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You also know the following: 



Therefore, y, I* l = — Yy r (ji j becomes 






V 1 (^] = “| £ c- Jf “p( ~ X /ix 0 
.1 'x, J '■■ - y ' 



What's look Jikt-j You can see a graph of ijijCjt) in Figure 4-2, where It has 
one rjuiJe (transition through the ri\is). 



Fi yore 4-2: 
The first 
averted 
slata si a 
quantum 
mechanical 
kronuric 
ou cil kilur. 




Finding the second excited state 

All right, how about finding and so on? You can find from this 
equation; 







<-X\o') lf> 
<x[*'f ff> 
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Substituting for a the equation becomes 



lt$w$ termite fwlyitomiete to find mg excited Mate 
You can generalize the differential equation for (ji„CO Mice this: 



V,Uj = 



] 



[ 



# ' 1 ( 2 ’ « ! ) ' " ' J 






To solve this general differential equation, you make use of the fact that 

( *-*•* i) ) 1 ■«’*(■%„= )"■( %■, ) 



is the nth tiermite poiynamiat, whirls is defined this way: 



H , ( x ) = ( -1 ]' fijep ( -v- ) ‘j -esp( —x- } 



rfi 



Holy irmtkRrcl! What do ! hr harm its polynomials look like;? Here's H M (xJ, 
j, anti so oil: 

^ H„u> = l 
^ H,(x) = Zv 
^ H-.(X') - -l..o - 2 
- S.X - 12* 

u* H,(.f) = 1 Hv 1 - * 12 

t+ H-(*> = - Ihllr 1 , ]2ily 



What does this buy you? You can express the wave functions for quantum 
meclianicaJ ttamioniu oscillators like this, using the hermite polynomials H.X.yj: 



where = 



)■ 



w 

I (mi j 



TIQh'et 
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And that’s what the wave function looks like fora quantum mechanical har- 
monic oscillator. Who knew It would've Involved her mite polynomials? 

You can see what 'M *) looks itke in Figure -1-3; note (hat there are two nodes 
here — In general, uT/1 for the harmonic oscillator will have n nodes, 




Putting in seme numbers 

The preceding section gives, you 4r Y (.ir)L and you've already solved fur E*,, so 
you’re on top of harmonic oscillators. Take a look at an example. 

Say that you have a proton undergoing harmonic oscillation with » - -1.58 x 
I |> J 1 sec' 1 , as shown in Figure H. 



Figure 4-4: 
A proton 
undergoing 
kiunuric 
:: 'i: llation, 
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WJhat are the energies of the various energy levels ol the proton? You know 
that in general, 

t, =||n + -| Jflitj rr = 0,1.2... 



So here are the energies of the proton, In megaeiectron volts (MeV)j 
tS F„ = ^ = l,!iO MeV 
t* F, = 3^ = i.5n MeV 
* F. I =^- = 1S<\ MeV 
^ E, = ^ = in, 50 MeV 



And so on, 

Now what about t lit 1 wave Jimnticms? The general form of is 



V ..[*)- 




-xV 



/lx 



! 

P 



where a 



■(; 



y I, 



n 



Sox., = 3-71 x](f m, 



Convert ail length measurements into Jem to meters (I fm * I x nr lf jii). giving 
you x„ *3.71 fm. Here’s %,(>), where x Is measured In femtometers: 



^ i>) = 



1)1 K * 



CX]J 



.. /27.S 



Here are a couple more wave functions: 

»••(*) " iraT* *(^ 7, H ~ s /n.a 

- J,-. . i x/ 3~t)' 
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Looking at Harmonic Oscillator 
Operators as Matrices 

Because the harmonic oscillator lias regularly spaced energy levels, people 
often view it in terms oi matrices, which can make things simpler. For exam- 
ple. ihe tallowing may be the ground state eigenvector (note that it's an inli- 
nile vector): 

1,5 > = 

i 

0 

0 

0 

0 

0 



Anri this may lie the lirst estatecl state: 



l l > = 
0 
1 
0 
0 
0 
0 



TICirT 
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And so on, The IM operator, which Just returns the* energy level, would then 
look like this: 







rs = 




0 


0 


0 


0 


0 


] 


0 


0 


0 


ft 


2 


0 


0 


1} 


0 


3 


(1 


<1 


tj 


u 


0 


h 


0 


0 



So N 2 -‘ gives you 



N|2-> = 



0 0 0 Cl 


0 


0 1 0 0 ... 


0 


0 tl 2 c 


1 


0 0 0 3 ... 


0 


Ct it fl 4) .... 


0 


0 0 0 0,,. 


0 


1 " ' r 


■ 
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This Is equal to 



N |2 > = 



0 




0 


1) 


0 


0 ... 


0 


LI 




Cl 


1 


u 


u ... 


i) 


2 




0 


L> 


2 


0 ... 


l 


0 




0 


0 


0 


3 ... 


0 


(1 


= 


Cl 


0 


□ 


11 ... 


0 


0 




0 


■ : 


0 


0 ... 


0 


0 






. 


_ 


. 




0 




, 




. 


■ ■ ■ 


, 


0 




. 


- 


- 




i. 



In other words. N!2> = 21 2>. 

f low about the a (lowering) operator? Thai looks like this; 



(I 


v'T 


a = 

0 


n ... 


(1 


(i 


J2 


o ... 


(1 


ci 


i) 


- 


<y 


tt 


0 


o ... 


a 


a 


0 


a ... 


(i 


a 


a 


fi ... 


- 


■ 


- 


■ ■- 


■ 


- 


- 


- 




. 


. 


. 
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In this representation, what Is a l !>? In general, a l rr> ■ rr 1 '- ft - 1>. so a 1> 
s h ould equal 1 0> , T ake a look: 



a|l> = 



it 


v'T 


a 


li ... 


Ci 


Ci 


0 




u ... 


1 


Q 


li 


(i 


s ... 


0 


0 


0 


0 


0 ... 


0 


0 


0 


0 


0 ... 


Cl 


it 


it 


!) 


a ... 


Cl 


■ 


, 


■ 


■ ■■ 


■ 




, 


r 


■ ■■■ 





This matrix multi plication equals- the folio win 13: 



ai>= 



m'T 




0 xi'l n 0 ... 


0 


0 




0 0 g‘2 0 ... 


t 


(t 




0 0 0 v‘3 ... 


0 


0 




0 tt 0 Ci .... 


0 


0 


= 


0 0 0 0 ... 


0 


(t 




0 0 0 0 ... 


0 






. 


- 


■ 







■ 
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In other words, a I ]> *IO>, )ust as expected. 

So how about the o (raising) operator? Here's how M works In general: o' 

■ (n i I ) J hi * 1 , In matrix terms, a looks like this: 

0 0 0 fl ...j 

m'T o o o .. 

0 Jz 9 0 J 

o n $ c J 

a o i> Ja J 

o o o o J 



For example, you expect that ti 1 1 1 = Jl i 1 Does it? The matrix multiplica- 
t ion is 



a'|l>- 



0 0 0 0 J 


0 


Ji o o o J 


1 


0 vf 0 0 J 


0 


0 0 -J' 0 J 


0 


0 0 Cl v'4 


o 


0 0 0 0... 


0 


i a j “ 





Copy righted material 




Chapter 4 : lack a ad Forth with Harmonic Oscillator# 113 

This equals liw lollowlnj;: 



u' |l 1- = 



ft 




0 0 0 0 


ii 


5 




v'T o o o ... 


1 


J2 




0 y'2 0 0 ... 


0 


it 




G 0 0 


Cl 


Cl 


= 


Cl 0 0 J4 ... 


n 


0 




Cl 0 0 0 ... 


ii 


’ 




. 


■ 


■ 






- 



So u l ' - 'U 2>, as it should, 

How about taking a look at the Hamiltonian, which returns the energy of an 
eigenstate. Hln» = K„ n--'? hi matrix form, the Hamiltonian looks like this,: 



M 

2 



EU 

o o n 

3 0 0 

0 5 0 

0 0 7 

0 0 0 0 

i) o (t a 



ited material 



So it yon prefer the matrix way of looking at things, that's how it works for 
the harmonic oscillator. 
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A Jolt of Java: Using Code to Solve the 
Schrodinger Equation Numerically 

Here's the OtiC-diiauuisiOuaj Schrudingtr equation: 



And ior harmonic os dilators, you can write the equation like this, -where 



k ' ~ 'im 



K - V| x I ) 

/a, i 



t Lx' 



I 3(1 



In general, as the potential V(r) gels more and more complex, using a com- 
puter to solve the Schr Winger equation begins to look more and more 
attractive. In this section. I show you how to do just that tor the harmonic 
oscillator SchrWlnger equation. 



Making your approximations 

In computer terms, you can Approximate v(jr) as a caller t ion of points, <jr ( . rp-., 
ip ,, ijr .. i|f 3r anti so on. as you sec in Figure 4-3. 



qIje) 
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ILaeh paint along vU'j — Vj, V:, 4* „ qr.,. Vs- anti so on — is separated from its 
neighbor by a distance, fr u , along the xaxis. And because ttufiiix is the slope of 
iji l.'ct, you cau make the approximatum that 

dy iy,„ -yt, 

d-< A, 



Iji other words, the slope, tkt itx is approximately equal to AyiAx, which ts 
equal to iii n . i - (= Ay) divided by fr„ (= Ax). 



Yon can rearrange the equation to (his: 






That's a crude approximation lor iy n , ,, given qq. 5a for example, if yon know 
qq, you can find Ihe approximate value of %„ if you know <t\yfdx in the region 
of S',. 

You can, of Course, find better approximations for qr fl . |. Iei particular, physi- 
cists often use the jVjjmtfPiM 1 aigurithm when solving the Schrfldingtr equa- 
tion. and that algorithm gives you , \ in terms of qi n and v„_ f . Here's what 
the Numerov algorithm says: 



^„ =■ 



V V 



f Sh;K[x\' / \ 
n 



V'- I ' 12 

V 






1 i h /ie... , (jn- )■ 



12 



Ur-,-™-* 

I n this equation, k, [ jc } = ^ 

and the boundary conditions are v(-*> = v(*0 = h. Wow, Imagine having to 
calculate this by hand. Why not leave it up to the computer? 

For a proton undergoing harmonic oscillation with a = 4.5ft x lb" sec - ', the 
exact ground state energy is 



/ fr 





E„ = M = i.i 
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You solve this problem computationally earlier in rh.kg chapter. The Following 
sections have you try to get this same result using the Nnmerov algorithm 
and 3 computer. 



Building the actual cade 

To calculate the ground stale energy ol the harmonic oscillator using tire 
Nnmerov algorithm, this section uses the Java programming language, which 
you can g£-it lor I ret from java . son . com- 

Here's how you use the program: You choose a trial value of the energy for 
the ground state, E h , and then calculate u(y) at x, which should he Kero — 
and If It’s not, you can adjust your guess for E and try again. You. keep going 
until %tfx) = I) {or it not actually (5, a very small number in computer terms) — 
and when it does, you know you’ve guessed the correct energy. 

Apprwmtttwg ) 

How do you calculate u(x)? After all, infinity is a pretty hlg Dumber, and the 
computer is going to have trouble with that, In practical! terms, you have to use 
a rum her that approximates infinity, In this case- you can use (he classical turn- 
ing point* of Flic proton — the points where all flic proton's energy is potential 
energy and it has stopped moving in preparation for reversing its direction. 

1 j 

At the turning points. .y„ E n = 1 Filial is, all the energy hi in potential 

energy), so 



jr, ■ 



2F / 



I MW 



I 



\ y i 



And this is on the order ol 15 Ecintomctcrs (Em)- so you assume that yfjr) 
should surely bo zero al, say, IIS fin. Here's the interval over which you cal- 
culate tf(x): 

I P*' C-i i, = -IS fill 

*" x-,. lV - 15 fm 

Divide this 50 f«i Interval into 201} segments, making the width of each sctr 
mem, h,, equal to (x riri , -a ll>lf ,),2(i0 = h : = Cl, 15 fm, 

Okay, you’re making progress. You'll start by assuming that ^(jr, ) = II. guess 

a value ol E,„ and then calculate (because there are 2U0 seg- 

ment, at x = x„ M . y,. - which should equal zero when you gel E : .. 
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Here's what the results tel! you: 

Corral: if is zero, or in practical terms, less than, say* vour 

maximum allowed value of v - 1 x t€“\ then you're done — the F„ you 
guessed was correct. 

t* Too high: If ahsCv^u) is larger than your maximum allowed i|», 

(= t x Ur \ and v 2)K is positive, tin; energy you c hus c lor E,, was too 
high. You have to subtract a small amount ul energy, AE — say I x 10" 
MeV — from your guess for Hit* energy: then calculate abs(uj,,r.) again 
and see whether it's still higher than your maximum allowed v, v : , „ x - If 
so. you have to repeal the process again. 

** Too law: If absCv»iu) is larger than your maximum allowed v. f- 1 x 
ifT’), and i|r»;i is negative, the energy you chose for F„ was too low. You 
have to add a small amount o* energy, AF., to your guess for the energy; 
then calculate absfqr-.,.,) ajjaiu and see whether It'.? still higher than your 
maximum allowed tji, i{r„ ai . If so, you have to repeat the process. 

So how do vuu calculate u au .? Given two Starting values, ht,, and u h use I I'll: 
Numtrtv algorithm: 



™. v MaT 



12 



V. ~ 



( ] + [Mr-a (■*}'] 1 



12 



V. 






M 






12 



Keep calculating successive points along yix): ^ i^ ;< , iy t> and so on. The last 
point 3s 

Okay, you’re on our way, You're going to start the code with the assumption 
that r|f n - 0 and 41 : Is a very small number (you can choose any small number 
you like). Because you know that the exact ground level energy is actually 
1.511 MeV, start the code with the guess that K tj - |.4£Kk)fl(X) MeV and let the 
computer calculate the actual valise using increments of iE = t x UT MeV. 

Note also fins equation depends an fc n {.v)\ k n _ ,(jri\ and k„ , Here's how 
you can find these values, where t,^ rrml is the current guess tor Ihc ground 
state energy (substitute ft, ft - i.amJn + 1 lory): 



ft ‘ | v )= 2 



F -m : x 1 / 



>py ric 
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And you know that K - 4.58 x I O' 1 sec J , so 

I V = ’Im 1 

^IjL i = M6 IVfeV tm ! 



Therefore. tt,'( v,) = ti.tiht,. —5.1)3 x 10 x.", where .v. for a particular se-t}- 
menljisjr. -jk t -x, 

Writing the cede 

Okay, now I’m going to pur together alt the info Irnm the preceding sec- 
tion into some Java code. You start with a Java class, up (for Srhrodinger 
Lqualinn), in a file you name sejnca: 

public elaaa ae 



} 

Then you sel up the variables and constants you'J] need, including an array 
For the values you calculate lor tji (because to find t|i„ . you'll have had to 
store? the already-calculated values of (jr„ mid J: 

public claaE se 
t 



double pal [ ] i 
doubl e ECu i : tnt j 
dc-uble Emin * 1.43 Or 
double xMin * -15.; 
doubl e Kite?; = 15, ; 

doub ! c h 7,c z o ; 

double EDsitu = 0 . OOmidOl j 
doubl e in^y Pal = 0.00000 00 L r 
int cumberDi vi e ion e - 200* 



1 

The se class's constructor gels run first, so you iniLialize values there, includ- 
ing ui,i vfx ) ■ 0) and y, (any small number you want) to get the calcula- 

tion going. In file main method, called after the cortsfruchv. you create an 
object of the se class and call, it calculate method to get things started: 
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nublic clrtSS n: CT- 
{ 

douii'l e p 5 i [ ] 
double ECurrenfc; 
doubl e Emin * 1 . 4 $ 0 j 
d'“ub lc xM in = - 15 - ; 
dou bl e jtWsx = 15 - ; 
double h 2 4X0 ! 

doubl e EDe 1 ta a 0 . D 0 0 0 d O L .? 
doubl e maxes i ■ 0.00000 001 ; 
int number Divi s ion e «■ 200; 

p u b lie se ij } 

[ 

ECurrent ■ Eminr 

psi * new double [numberPivi si or.s +■ 1 j 
psi[0] - 0? 

psi p ] = -a.COOOCQQQl ; 

pa i T nUmtei D i v i is ie fljs ] = L . 0 r 

hZeto = ■; j-LHa.it - jtHi/Lj / LiuinbiitDlvi alone- ; 

public static void rnai n t 5 tr ing [] argvi 
( 

3e di = net* 1 (1 r 
de . calculate | J ; 
i 



The nial work Takes place it: the calculate methml, where you use Hie current 
yiucss tor l he ciierjjy and calculate 

II abs(nr-.r,ij is Jess than vuur maximum allowed value of ip, ^ !ll4X , you've 
found the answer — * your current guess for the energy is right on. 

* It abs(yM,„) is greater than y.„ ir and Is positive, you have to subtract 
it from your current guess for the energy and try again. 

** if absOf is greater than y-, m and is negative, you have to add af 
to your current guess lor the energy and then try again.. 
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Here's what all this looks Like in code: 

pu b 1 ic v<3 id co l cu la to [ i 
t 

wh i 1 e ( Ma t Li . tiba ( p a i [ number Bi v i » i one | | > maxPel ( 
for (int i ^ lj i -cnumi-er Di vi s ion a ; i++) { 
psi[i * 1] ■ calculateNextPsi ( i | 

) 

if (pni [ numbc rB i v i nlbn si > 0.0) f 
ECuiczent = ECui i crit H : ■ i ■ 1 1 1_ : j 

} 

else ( 

ECurrent ■ ECurrent 4- EDelta: 

) 

Sys ben . out . pr i n 1 1 n | i Ps i 2 □ 0 ; i + p# i \ nizmbe :!! ; yi I ' 
+ 1 Bl 1 + pound | ECUc i trLt) ) j 

} 

system, out, pr in tin (l\nThe qr-our.d state eneroy is i + 
round | ECurrent ) i i HeV.S); 

1 

Note that the next value of f (that is, f„ . i) is calculated with a method 
named raiculateNextPai. Here's where you use the Humerov algorithm — 
yiven T|) n , nr„ . ,, you ran oak nlale 



public double calcula^eNex tPai (but n} 

C 

double KS-nUMi ousOne s ca lcul at eXSaua red ( a - 1J ; 
double K5qH =■ oalculatBKSqtiared (n) j 
double KSqNFlunOno = calCulatcKSqoorci I n + II ; 
doiib le nexc Psl = 2.0 * ( l . D - | S . 5 * hzero * hs er o * 
KSqJi / 12. D I! * psi [n? i 
nextFsi • nextPsi - |1.0 +■ ^hZero * liZero * 
KSqHKirausOns / 12,0)1 * psi[n - 1] ; 
nsJctPsi = nextPsi / 11.0 + (hEero * hEoro * KSqNFlusQne 
t 12 - 0 )).; 
return next Ps i r 

! 

Finally, note Hut to calculate . ,, you oeetl k„. k, and k„ . h which you 
find with a method named mtetflaieKtyuamt which uses the numeric values 
you’ve already figured out lor this problem: 

public riaubl b cn 1 cuTatsKSqUia.r«d{ iflt O.J 
[ 

double x = IbSero + n) + xHin; 

return ( {<0,05) * ECurrent! - |!x * x) * &.6ie-3i)r 

1 
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Whew. Here’s the whole program, se-jovo: 

public clflpn ! 5 C 
{ 

double pai [ ] t 
doubl e ECurrent ; 
double Emin ■ 1.430; 
double 3cWi.n = -IS.: 
doub 1 e :■■..'■! = 15.: 

doilbl e !hZe£o; 

dc-ubl e EDe 1 1 a ■ 0.0000001 = 
double maxes I ■ 0.00000001: 
int n.umbsrDi vis ions = 200; 

pu bl 1 o ge / [ 

{ 

ECurren t m Em in r 

psi » new double [numbers! visions i I); 
psi [0] = 0; 

psi rn = -a . zoooDDoai.: 

psi [number Divisions ] = 1 . 0 ? 

hZero = (vHsk - xHiii > / nuriber Divisions; 



public static void mni n ( 5 tr ing [] argvj 
[ 

a e de = new a.* t J ^ 
de . calculate [ } j 



publ i c vo id ca. 1 cu la te I ) 

[ 

wb.l la (I4a th . abs (asi[ mmberDi vi siona J | > maxPa 1 ! ( 
for (Int 1 = lj i ^number Divisions ? i++H 
ps i [ i + 1 ] - cal cu La teNext Es i ( i ) ; 

) 

if (psi [number Division si > CKO) { 

ECurrent = ECurrent - EDeltd; 

] 

elee ( 

ECurrent * E Currant + EDelta; 

j 

Sys teta - out- . pr in t In 1 1 Ps i 2 0 0 ; i + psi; nemhe r Div i s i on s 1 
f L E: 1 + round {ECurr ent ) |i ■ 

1 

System. cut .printing i \nThe ground state energy is 1 +■ 
round | ECurrent ; * i HsV, ij; 
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pub I i c doubl e ca lew la “ eKijq'ja rediint n I 
t 

double y, = Ibzero * d) + xMinj 

return < ( -;rj .d-5) * ECurrent! - ( (x + x) + E. . 63e-3 ;i ;i .r 
I 

publ i c doubl e ca lew la - t F 5 i ! i r.fc jji 

[ 

double KSqHMinueOne = aalfmlateSSqtiared{fl - I); 
double KSqll - calculateKSquared in) ; 
double KSqMFlusOne « calculateKSquer ed | n + l!r 
double jioxifeP*i - 2,0 * (1,0 - {5.0 * htero * bZwzo * 
K5qK / 13,01} * psilnji 
nKKh.P.^i = n.KxtP.r:i - { } .0 1 - {Hlfiro * h7,rtro + 
KSqNKinusOne / 1S.0) } * pal [ji - l]r 
nextFsi • nextPai /{ 1 .0 +■ thzero * hZero * KSqlvPlusOne 
? 12-0)); 
return n ext Ps i : 

:■ 

public double round ( doub le val) 

( 

double divider ■ 100000.? 
val = vaI * divider, 1 

double temp = Math . round (vftlj ; 

1 ■ '! 'jin ^doublel Lemp i divider; 

:■ 

} 

Okay, now you can compile Ike code with javac. (he Java compiler (it javac 
Isn't in your computer’s path, fee sure to add the correct path to your com- 
mand-1 ine com man d, S«C h as (1 i >l; : \ j a va \ h an \ j r sr-inirn). 

C : > jiftvftc >s«s e j a vh 

Thi& creams se.rfaffs from .<x.jninj. and ynu can run .ce.rtaff.t with Java ilselF 
(adding the eorretl path il needed): 



C->java se 
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Running the cede 

When you run the java code for the harmonic oscillator Schrodinger equa- 
tion, It displays the successive values ol v ; „: as It adjusts the current guess 
for the energy as it narrow®, In on the right answer — which It displays at the 
End ol: the run. Hercf s- what you see: 

C : > j ava se 

PSI200: -1.O3OM44097337778E-4 E: 1.4® 

P5:200: - 1 . 050354423295303E-4 E: 1.4® 

PS:2O0; -1. 0503444368533 10SE -4 E; 1.4? 

PSTHOQ ! -1-05033445043604555-4 Ej 1 .49 



PSX2QQ; *6 . 12820872 814324E- B E; 1.50056 
POISOtl : -6- 0311,2752 1J56655E-E E; 1 . 300 6 6 
PSISaa? -5. BS40463493D7554E-I5 1!: 1.50066 

ESX200i -S,S36&S5iso6Daai5E-S Bi 1,50066 
FSXSDOt -5.73988397946 177BE-8 Hi 1.50066 
P5:200: -5.6428029151212084E-8 E-. 1.50066 
PS:200; -5.54572 132 5289922 4E 8 Ei 1.50066 
P3TS00: -5.44064 OSO 6 6 5 ] 9 9 £ 6E - B E : 1.50066 

PST30P* -5.35l55'37022Dl63$E-S Bi 1,50066 
ps:200i -5,25447S733S76335S-8 Bi 1,50066 
FSXSnOt -5.1573 977143 26237 E-8 Hi 1.50066 
FS:200; - 5 . 060316SG1012202E- 8 E; 1.50066 
P5I2O0 ; -4 . 963235841725704E-8 B; 1,50066 
PSISOO: -4.BB6154915227413B-S E: 1.50066 
PSI200 : -4.769074O41B271214E-0 Bi 1,50066 
PEX200i —4 . 67LB , 932032691944E-8 Hi 1.50066 
P5I200 : -4.574912 36 8974434 E-8 E: 1.50066 
PS = 2 00 ; -4 , 477 83 15322 5 87505B-8 E-. 1.50066 
PSI200 s - 4 . 3 8 0 7 5 07? 04 7 6 514 E - 8 Ei 1.50066 
ps:2aa* -4.2B3670O57B3992E-S E = 1.50066 

P6l200i -4.1665e934521757SE-8 Bi 1,50066 
PSI2aa: -4 . 0B9bQ£ 5S731B4QG4.E-8 Hi 1.50066 
F5:20Os “3 . 99242703522620 lE -8 E; 1.50066 
PSI200 ; -1,8953472 07 3066213E -8 E; 1,50066 
RSI200 s -3.79826665057731E-S E= 1.50056 
PSI2aa: -3 .70113603 B5D2B2GE-B E- 1 . 500 6 6 
PEX20Di -3. 604105453620 2 66E-S Hi 1.50066 
PSI20a ! -3 .50702 494 9809914E-B E: 1.50066 
PS 1 2 00 ; -3 ,409944421787517415-8 Ei 1 . 50066 
PSI200 s -3 , 31286 3911389194E-B E; 1,50066 
PS 1 2 00 : -3 ,215 7 i3 471996181 5E -8 E: 1.50066 
PSI20G : -3 . 1 1 B 7 03 00 5 9 9 0 2 $5 6E - B E 1 1.50066 

PSI200: -3 . 021622619b94b36E-B E: 1.50066 
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PS : 2 OO ; -a.92454219B5136167B- : 8 B: 1,50066 
PSI200 j -2 r 6274618172-37 53 S58”B Is 1,50066 
ps:2aa : -3..73Q391 53 443697 03 E-B E: 1 . $0066 

PEX2&&I -2,€333011960G9577E-€ Bi 1+50066 

PS 1 2 □ 0 r -1.336220SS3B3iaa65E-B E: 1 . 50066 
P5:20O: -2. 4:391406320658141-8 E: 1.50066 
P5I200S - 1 ,3420604246230753-6 Es 1,50066 
P5I20&S ”2,24498Q221960756E-e Es 1,50066 
PSI2aa ! -2.l479aC0O52-47249E-S Ej 1.50067 
PEEJODi -2.050619828502253215-8 Bi 1+50067 
FS=200i -1.953739Y616B23192E-B E: 1.S0067 
P31200 ; -1,856653 0-0 0 2 0 6 6 1 0 5E • 8 Is 1.50067 
PSI200 : -1.75957S5286272332:&»8 E; 1,50067 
P512Q-0: -1, 662499470377955 5K-8 Es 1.50067 
PSI2Q0 : -1,56541946169206 2 E - B E: 1.50067 
PS 12 00 : -1.&6B339473D036424E-B E: 1 . §0067 
ps:200; -1. 371259459203416 5E-B E: 1,50067 
PS=200s -1 ,27417351596385978-0 Es 1,50067 
PS1200; -1,177099 62 296634 BE -B Es 1,50067 
Ps:200: ”1-0000197 14 27336S3E-P Es 1.50067 
psxaoai - 9 . B 2 93 9 7 98 5 2 9 6 3 2 E - 9 B: 1.50-347 

The ground ata^e energy is 1.50067 MeV. 

And there you have si — the program approx I mates the ground state energy 
as 1,50067 MeV, pretty darn close to the value you calculated theoretically 
before: 1 , MeV, 

Very cool, 
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Part III 

Turning to Angular 
Momentum 
and Spin 




The 5 th Wave 



By Hich Tennant 



“Great— di££erentiai equations brought tj& 
KevHoYi h s Lavi ot \5niversst Gravitation, 
MaxvtelTs iieki equations, and new SHuatVs 
Hate oi Hair Loss.’ 1 







fa this part . . . 



T hings that spin and relate — that's the topic -of this 
pari. Quantum physics has all kinds of things to say 
about how angular momentum and spin are quantized, 
and you see it all in this part. 
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Chapters 



Working with Angular Momentum 
on the Quantum Level 



h This Chapter 
Angular momentum 

Angular momon-tun* airac! Hits Hamiltonian 
Matrix representation qE angular momentum 
Eigenfunctions of angular momentum 



» n classical mechanics, you may measure angular momentum by attaching 
m- a golf ball to a string and whirling It over your head. In quantum mechan- 
ics, think In terms of a single molecule made up of two bound atoms rotating 
around each other. That's the level al which quantum mechanical effects 
become noiiceabie. Anri at that level, it lianas out that angular momentum 
is quantized. Arul since that has tangible results in many cases, surii as the 
spectrum oE excited atoms, it’s an important topic. 

Besides having kinetic and potential energy, particles can also have rvttilional 
energy. Here's what flit Hamiltonian (total energy — sec Ctiiaptcr 4) looks 
like' 



II- 



I- 

21 



Here, L Is the angular momentum operator and [ Is the rotation moment of 
Inertia. What are the eigenstates of angular mo men tool; if L is the angular 
momentum operator, .and t is an eigenvalue of L, then you could write the 
following; 

i-* Incomplete! 



But that turns OUt to he Incomplete because angular momentum is a vector 
in three-dimensional space — and It can he pointing any direction. Angular 
momentum is typically given hy a magnitude and a component in one 
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direction, which Is usually the Z direction. So In addition to the magnitude i, 
you also specify the component of L In the Z direction. L- (the choice of Z Is 
arbitrary — you can just as easily use the X or Y direction). 

If the quantum number of the Z component of the angular momentum Is des- 
ignated by m, then the complete eigenstate is given by If, m>, so the equation 
becomes the luJlowirig: 

H | J , m :>= j=y l,m > 



That's tin: kind ul discus a ion about eigenstates that t cover hi this chapter 
and I begin will] a discussion of angular uiOiiLuituifl. 

Ringing the Operators: Round and 
Round with Angular Momentum 

Take a look at Figure 5 - 1 , which depicts a disk rotating in 31 ] space, because 
yoii’r? wording in 'M'i, you have So go with vectors to represent both magni- 
tude and direction. 



i 



Fi gore 5-1: 
A rotating 
disk will 
angular 
momentum 
vector L. 




■i' 
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As you can see* the disk's angular momentum vector. L, points perpendicular 
to the plane of rotation, Here, you can apply the right-hand rule: II you wrap 
your right hand In the direction something is relating, your thumb points in 
the direction of the I. vector. 



Having the L vector point out of the plane of rotation has some advantages, 
For example, it something is- rolaling al a constant angular speed, the L 
vector will be constant in magnitude and direction — which makes more 
Sense than having tile L vector rotating in the plane of the disk's rotation and 
constantly changing direction. 

Because L is a 3 D vector, it can point in any direction, which means that it 
has*. >■, and * components. U. and L f which aren't vectors, just magni- 
tudes). You can see l.. in Figure n-i. 

L is the vector product of R (position) and P (L = Rx P). You can also write 
L,, and U at any given moment in terms of operators like this, where P„ 

P„ and P. are the momentum oppynU'ii^ (whir h return the moments:™ Ln the x, 
y, and 2 directions) and X, V, and Z are the posifiun operators (which return 
Hie position in the x. v. and 2 directions): 

1 ^ L, = YJv-ZP, 

^ L. - ZP, -XP T 
^ I,. = XF,.-YF, 

You can write The momentum operators F v . F,, and F. as 




P. 





Therefore* substituting these operators in the L„ ],, and U. equations, you 
nan write the equations as 

*'l,-h«(x£-y£) 
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Finding Commutators of L# L v , and L z 

First examine L,, L,. and L. lay takin g a look at I low they commute; IE they 
commute (lor example, if [L,, L] = (I), then you can measure- any two of them. 
(L, and L r for example) exactly. SE not, t lien they're subject to tlie uncertainly 
relation, and you can’t measure them s i mu llaneous ly exactly. 

Okay, so what's the commutator of L and L,? Using U - YP ; - ZP.. and Ly - 
7.P...-XP., you cam write the following; 

1 1.,, f.,1 = IVP.-ZP,, ZP T -XPJ 

You can. write This equation as 

[L, L/j = |YP., ZPJ- [Yf\ XPJ - [ZP jr ZP,] 4 [ZP„ XPJ 
-Y[P,ZPJP. *X[Z,P. J ]P J 
-ift(XP,-YP,) 

But XP V - YP V - L : , so [U, L,| - ml.. So L^. and L, don't commute, which means 
that you can’t measure them both simultaneously with complete precision, 
You can also Show that 1 1^, L.| = and 1 1__ I.J = /fil... 

Because none nl the components oF angular momentum commute with each 
other, you can’t measure any two simuitaneousty wilii complete precision. 
Hats. 

That also means that the L„. L,. and L- operators can't share the same eigen- 
states. So what can you do? How can you Find an operator that shares eigen- 
states with the various components of L so that you can write the eigenstates 
as 1 1, m>7 

The usual trick here is that the square of the angular momentum, L' , is a 
scalar, nol a vector, so It’ll commute with The L„ L,. and J_ operators, no 
prohiem: 

1 ^ IL-, LJ = Eh 

^ [L“, L, | - € 

[L a , L : ] - 0 

Okay, cool, you're making progress. Because U. L,, and L : don’t commute, 
you can't create art eigenstate that lists quantum numbers for any two of 
them, But because i. J commutes with them, you can construct eigenstates 
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than Tiav-e eigenvalues for L" and any one ol L, . L, and L : . By convention, the 
direct Ion that's usually chosen Is 1... 



Creating the Angular Momentum 
Eigenstates 

Now’s the lime to create the actual eigenstates, if, m > , of angular mo men turn 
states in quantum mechanics. When you have the eigenstates, you 'll also 
have the eigenvalues, and when you have the eigenvalues, you can solve 
the Hamiltonian and get the allowed energy levels of an object with angular 
momentum. 




ikm'l make the assumption that the eigenstates are If, m>; rather, say they’re 
( a, Jj =-, where the eigenvalue- of Lr in L ; let. fU = h'ncl a, ji>, So Uie eigenvalue of 
L is O'nc, where you have yel to solve lor Similarly, the eigenvalue of L.. is 
L - 1 cl, p> s ftp a, p». 



To proceed further, you have to introduce rttfs/rtgand lowering operators (as 
you do with the harmonic oscillator jo Chapter 4). That way, you can solve 
for the ground state by. for example, applying the lowering operator to the 
ground stale and setting the result equal to ?,ero — and then solving for the 
ground slate itself, 



In this case. Hie raising operator is L. and the lowering operator is L__ These 
upcralors raise and lower the L. quantum number. In a way analogous to- the 
raising and lowering operators iti Chapter 4, you can define the raising and 
lowering Operators this wav: 



I y* Raising: L, ■ L, + r'L, 

W Lowering: L - U - tl. 

These two equations mean that 



L, = .;IL. + (. ) 






You can aLsu see tliat 



d- 



dlUl Id 



L L = L. ♦ L, 5 - nL. = L a - 4 s - fiL. 
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That means the lolloping are all equal to L 
v 1/ = I...U t L J -ftL_ 

n* L- = L.L. + L J + ijL_ 
v* L- = 1 f .L. L * L L ) * I J 

I 

You can also see that these equations are true: 

[L ', L] - 0 
^ [L., LI :m,. 
t" | L., L.J = ±fiL, 

Okay, now you can put all tt*is to work. You 're yetting to the gonti stull. 

Take a look at the operation of L. nil In, fk: 

L. I if. JS> ■ 7 

To see what L. In, [5> 6s, start by applying the L. operator on It like this: 
l_l. I a, |i> -7 

from [L.. L. j = you can see that L. L - L, L. = tiL.. so 
L. L, a. fk = L L. In, |’k + JrL. I a, jk 
And because L. I tt, p* = p. you have the following: 

I. L. ,f. p> . ^i|t + J >L. hf. J5=- 

This equation means that the eigenstate h l u L |5> Is also an eigenstate of ihe 
L. operator, with an eigenvalue of f|l ■ 1). Or In a more comprehensible way; 

I.. I ra, JJ> = -cl fit, fl t tv 

where r is a tuns taut you find later in "finding the Eigenvalues ol the Raising 
and Lowering Operators ~ 

So the L, operator lias the effect of rasing the (i quantum number by I. 
Similarly, the lowering operator does this; 

LlB,ft>-ClH, [l-l* 
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Now take a look at what 1/1. i n. 15 > equals: 

!_“]*, a,|iv = ? 

Because 1." is a scalar, it rommufss with everything. L" L. - L. L." = 0, so this is 
true: 

L.ajiv. p> L. L- 1 rv. p> 

And because L n. [5> - r&, |K you have the following equation: 

Ijl.ifi, [5> i- ufi" L, I ft. |5> 

Similarly, tha lowering operator, L„ givR&ynu this: 

I." I. Jr, |i > = ft.tr L_ I nc. [i-- 



So the results uf these equations mean that the L. operators don't change the 
n eigenvalue ui I n:. 0 - at all. 

Okay , so Just what are « and ft? Read on. 



Finding the Angular Momentum 
Eigenvalues 

Tilt eigenvalues of tin: angular mumeiil'.iin are the possible values I .he angu- 
lar lUOmenlurl] Cftii take. 60 they're worth 1 hiding. Lefk take a look at lluw lu 
do jufit th.il. 



Deriving eigenstate equations 
With and 0^ 

4 | n n *1*1 

Note that |/ - [ = ]. v " + which Is a positive number, so [.” - I./ £• it. That 
means that 



«x.$ L J -Lj« f p^>0 
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And substituting in L «, fi> ■ and l : ' I n, p> » lift gives you tills: 

<«. $ ]. J - iJ Ut, |l> = *■<« - ii-j i n 

Therefore, ct a |i" So there's a maximum possible value ol p, which you can 
cal] fi,,,,...- 

You can be clever now, because there has to be a siale f a, p mx .s- such tliat 
you can't raise |) any more. Thus, Lf you apply the raising operator, you get 
zero: 



- ie ■5|ra,5£ > - 0 

Applying the lowering operator to This also gives you zero: 

L_L.la,|W=0 

And because L.L. = L -L - tlL_, Thai: means the following is true: 

a-_L/-r ? L : )l rt, p M ,„>-0 

Putting in L a, * aft and L- : l it* |Jv> - |S Mll ^i gives you this: 

« = LA.,vT 11 Hi 

Cool, now you know what n is. At this point, its usual to rename |i„. lT as l and 
p as rrr. so ( a, p» becomes 1, m> and 

I t* L S U. nr* = i(j * Lift 3 ! ^ m> 

** L : l l m> ■ fipi 

You can say even more, In addition to a there must also be a |5 Mlhl such 
that when you apply the lowering operator, yon get zero, because you 
can't go any lower than |J r : 

LI/. p m h,> = 0 

And you can apply L on this as well: 

L,Ljlt|W^:i 
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From LL. *L - L. - bL : . you know that. 

which gives you the following: 

O - P, ' - Pm„, * 0 

a “■ Piiiln — P^iO 
rl * CP.-Hl, " 

And comparing this equation to a = |i r71 . L , j[i iniJ « I) = 0 gives you 

ft™ 



Note that because you reach lot, p by n successive applications ol L_ on 

1 «■■ P„„ L ..=-. you get the following: 

Pfv-lv - P.-|ii, + ft 



Coupling these two equations gives you 




!!•, „ - % 

Therefore. [i 1M „ ran be @3 1 her an integer or half an hunger (tie pending on 
whether 12 is wen or odd}. 

Because I = m = p. ajld. n is a positive number, you can find that -i £ m m l L 
So ik. iv.' you have it: 

The eigenstates are l f, m>. 



v 1 The quantum number of the total angular momentum isf. 



^ The quantum number of the angular momentum along the z axis is nt. 

v* Ldf, m> = KJ + l)lf. m*. where t - if. 1 . Vi. ... 

>* L.I t me = fll I (, rtt>, where ar = -i r -(l - L } I - L, L 



u* -I < m £ i 



For each /„ there are 21 * i values of rw. For example, if / = 2. then ni can equal 
-2,-1, 0, ) . or If l ■ '/i, tliea ?rr can equal -'.h, 'h., ■!, and % 
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You can see a representative L and L. in Figure 5'2, L Is the- total angular 
momentum and L. is the projection of that total angular momentum on the 
zaxls. 




Getting rotational enema 
of a diatomic motecuie 

Here's an example that involves, finding Hie rotational energy spectrum o( a 
diatomic molecule. Figure 5-3 shows lire setup: A rotating diatomic molecule 
is composed of two atoms with masses wr, and mj. The first atom rotates at 
r h y... and the second atom rotates, at r - i-„ What's the molecule's rotational 
energy? 
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Figure 14 : 
A rotating 
diatomic 
mnlacule. 




The Kamil Ionian (as you can -see at the chapter intro) is 




1 Js the relational moment of Inertia, which is 

I = mi* - ;' * fJV- J = ,0 f 

wtiure r = I r - r, and 1 j - — ' V ' -. 

r rrt , + ft? . 



Because!- = !•»>. E, = Therefore, the Hamiltonian becomes 
L/__ LL 

n 2 fir 1 
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So applying the Hamiltonian to the eigenstate*, 1 1 m>, gives you the following: 



: L -l T {ft 
l).\r 



And as you know. L’ if ju> - i{i » ])fi’if m>, so this equation becomes 



S — |Oh* 



2 pr 



2 fir ■ 



And because H 1 1. m> = E. !, m»„ you can see that 

2 fjr 

And that’s the energy as a function oJ /. the angular momentum quantum 
number. 



Finding the Eigenvalues of the Raising 
and Lowering Operators 

This section tanks at finding the eigenvalues of the raising and towering angu- 
lar irinmenlum operators, which raise and lower a state’s z component rd 
angular momentum. 

Start by taking a look at L.. and plan to solve tor c: 

L. H m ■ - c t in + 1> 

So L /. m> gives you a new state, and multiplying that new state by Its irons* 
pose should give you <r: 

| L. 1 L rn») ! L. S f. m > = r" 

'Jo see tills equation., note rtiar {L. I, m») L I L, m> = c 1 <1, m +1 I. m ~ ] > = tr. 
On the other liEmd, also note that(L. 1 1, mv-l’L, 1/ m> = « /, m I L„ L_ I (, m-. so 
you have 

■ I mil... I if, m> * c " 
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What do you do about L L? WelL, you see earlier In the chapter. In "Creating 
the Anguiar Momentum EsRenstates that this is true: L,I- = I- - 1= * fit-.. So 
your equation becomes the FoLlnwing: 

■= I, m I L - — L_“ + ti L. [i l, m > = c 

Great! That means thm c is. equal to 

e - (4 ml L" - L/ * #iL-i l. m*)'" 

So what is (/■■!. m I L'' - L : + M : 1 1, m> ) "? Applying the L : and L : operators 
gives you this value for c: 

l: = h [ )(} ♦ ])-m {m t 1)| 

And that's the eigenvalue? oE L . which means you have this relation: 

L. 1 1. m=- = * 1) - m{ m * 1}] ' 1 1. m * 1> 

Similarly, you can show that L gives you tire following: 

L s I, m> - n\J(J + I)- - J)] l/l l f, m - |> 

Interpreting Angular Momentum 
With Matrices 

Chapter 4 covert a matrix interpretation uf harmonic Oscillator States and 
operators, and you can handle angular momentum the same way (which 
often makes understanding w hat's going on with angular momentum easier). 
You get to take a look at the matrix representation of angular momentum on 
a quantum level now, 

Consider a system with angular momentum, with the total angular momen- 
tum quantum number 1= 1. That means that m can lake the values -ML and 
1. So you can represent the three possible angular momentum states like tins: 

|L,- l >= 

0 

0 

1 



J py 
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iUJ>= 

0 

i 

0 

I 

(I 

a 

Okay, so what are tlie operators you've seen in this chapter in matrix repre- 
sentation? For exam pie. what is L J V V'ou can writ e L" tJiis way in matrix (omi: 



V 





<u r. : i,o > 


<1.1 f. : 1-1 > 


],o t; i.] > 


h, 0 |L' 1 1.0 > 


<1,0 J. |]--l > 


< 1-1 v\u> 


<] -l|L l |],0 ^ 


< ].— i|l = i.i > 



Okay. < l , ] ! V 1, 1> - i(S * ])«' ! - 2t\: < I, 1 i L J 1 1, 0> - 0; < L 0i 1/11, 0> - 2ft 5 '; 
and so on; Therefore, the preceding matrix becomes the following: 



].' = 



■>tr 


ft 


0 


0 


2/r 1 


0 


0 


0 


2 fl- 



And you can also write this as 



L : = 



2 



] n 

i) i 



!) 

0 



It D E 



So in matrix lurm, I he equation! L" ! L, 1> = 2ii ? l ], 3 » becomes 





1 


o n 


i 


2 it' 


0 


i n 


0 




(t 


ft i 


Q 



py righted materi 
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This equals llae following: 





1 




I ft 0 


t 


'Iti- 


ft 


= llr 


0 1 ft 


0 




ft 




ft ft 1 


ft 



How about the L, operator? As you probably know (from the preceding sec- 
tion), r., i nj> - 1 l}-/n(ru i l)| ,: l|, ijj * l > r tn this example,/ - l and m 

- 1, ft, and -l, So you have the following: 

I I* I..I l, lv = I! 

^ L. 1 1,-1 >- g5fi|l.0> 



So the I-, operator looks like this in matrix form; 



L.= 



0 1 0 



J2F< 



ft ft 1 
0 0 0 



Thcrdorc. L. ! 1, 1 1 woliIe! be 
L. |L,lV> = 





P 1 ft 


0 


■M>; M u 1 


1 




0 ft ft 


0 



Anil this equals 



L.iUt> = 





i 




ft 1 ft 


ft 


Jzn 


ft 


= J2h 


ft n i 


1 




0 




0 ft ft 


ft 



In other words, ^ h 1 ■ 1 L ] C . 
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Okay, what about L ? Yon know that l i i m - • * I) - m{m - I)] 1 -- 1 f„ m 

l». In this example, l - 1 and m *1,1 and -I. So that means the following: 

^ L_ 1,1 >= J2h L, 0 

^ L l.n >= - l > 

V L_ | .1.-1* =0 

So the L. operator looks like this in matrix form: 

L_ > 

o o o 

■Mh 1 (I It 

o i & 

That means that f- II, I > would he 



L_ 


tl’= 




0 U U 1 


■J2t j 


1 0 11 0 




& ] U 0 



This equals 

L_ |U>= 

U |1 0 0 3 

J2h 1 = .,' 2 /r I 0 00 

& ti L 0 € 

'►Vh id i lulls you that 
v'2A |1.0;i-L 1 1. 1 > 

Just as you'd expect. 
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Okay, you’ve found L . L„ and L. Finding the matrix representation of L a is 
simple because 

I t^ ,'M 1. 1> - ]. : l|, L> 

() = L f l Ml? 
t ^ -S'. ],]=■ = L.I ], -i> 

So you have that 

■-,= 

L 0 0 
a y o a 
oo-i 

Ti ms IL 1,-1 -- equals 



L 


0 


0 


n 


0 


0 


0 


0 


0 


0 


-a 


-i 



And this equals 

Ml.- '>= 

0 10 0 0 

-h ft = k IJ 0 0 0 

i U 0 =i -] 

So L = i.-K 

Now what about finding the L, and L, operators? That's not as hard as you 
may think, because 

L, = jj^lL. + L-) 

and 

L> = -L,) 
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Take a look at L, first. L r equals 
L. - 

n i (! 

VZfiO 0 1 

0 0 0 

And l._ cqi la Is 
L = 

0 0 0 
M 1 0 0 

0 L 0 




Cool. This, is going prelly well — how about calcuLalijig [L,. L, )? To do that, 
you need to calculate [I*. LJ ■ L.U - L.L,. First find L.L.r 



L , L 



ft 


] 


v\ 


i> 


-i 


0 


i 


ii 


l 


f 


0 


- 1 


Vi 


i 


u 


n 


i 


i> 
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This equals 



L..L,- 



h 


L 


0 


0 -> 


0 


l 


0 -i 


i 


0 


l 


r 0 


J B ! 

-i\ M Y 


1 


0 & 


n 


i 


a 


n i 


0 


i 


i> -i 



AtkI similarly, L,U equals 
L..L, ■ 



0 -l 


0 


0 


] 


0 


i & 


■ t 


1 


0 


1 


0 r 


0 


0 


1 


0 



And this -equals 
L.L.- 



n -r 


Cl 


0 


1 


a 




-d 


i> 


-? 


r 0 




I 




1 






0 


& 


rj r 


u 


0 


1 


ti 




■ 


0 


F 



'll U i) 

tr- 

y 0 0 0 

1 0 0 -2i 

And this equals 

| L..L. ]= L 1. - L 1. = 
i o a 

e® o o 
a o -i 
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But because 



1 0 0 
13 0 13 
0 0 -I 



yo 5 i can rewrite | L r .L„ J — L,L r L , L , 
1 1 0 0 



h 



£i £i D 
Cl 0 -3 



like lkia:[L Ll L 1= L,L* - L.L J = rfaL, 



Cool, sn [ L LJ = rtiU. 



Rounding It Out: Switching to the 
Spherical Coordinate System 

So far. this chapter has been dealing with angular momentum liras and kels: 

v'2Jij ],«>=: L |],l> 

The charm ol bras and kels is ihsi Nicy don't limit you to any specific system 
ol representation (see Chapter 2). So you have the general eigenstates, but 
what are the actual eis&nfurtcttout of L and I/'. 1 That is, you're going to try to 
find the actual functions that you can use with the angular momentum opera- 
tors like l ' and 

To find the actual eigenfunctions (pot just the eigenstates), you turn from 
rectangular coordinates, *. y, and to spherical coordinates because it 'll 
make the math much simpler (after ail, angular momentum is about things 
going around in circles). Figure 5-4 shows the spherical coordinate system. 



J py 
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In tile rectangular tCarmsianJ eonr-dlnam system, you use x, y, and z In orient 
yogrseLF. in the spherical coordinate system, you also use three quantities: r 
9, arid cj*. as Figure 5-4 shows. You can translate between the Spherical coor- 
dinate system and I lie reclangular one this way: The r vector is the vector to 
tile particle that has angular momentum, 9 is the angle of r from the z axis, 
and q is the angle of r from the x axis. 

I k" x -rslni) c -0341 

/ i'a r s ifiil 5. in ci 
z = r ros9 

Consider t he equations lor angular momentum: 

L = YF - 2 P = -fllf T ; -' - Zl 

I, dz dy J 

L.-a-.-XK, =-*(g-gj 



L , *(» •) 

When you take the angular mu men turn equations wit)] the Spherical-COOrdL- 
nalu-syslem conversion equations, you can derive the following: 




>pyrig 
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+ jl. 




j'cosP h l I 
sinfi J 



I. =L -fl- 



=- w fl 



mill i 

kmi# ctyl 



Okay, Ihese eq ual loos look pretty Involved. Ei.it I here's one thing to notice; 
They depend only on ft anti 4>, which means their eigenstates depend only on 
ft and not on r. So the eigen Functions Ql the operators in the preceding list 
can be denoted like this; 

< 6 . 0 i i 1 , m -■ 

Traditionally, you give the name Y,„.(0. <$<) to the eigen fund ions ol angular 
momentum in spherical coordinates, so you have Hie following: 

Y 1IS («, o) - <0, 4' S. m> 

Ail right, time to work on linding the actual lorm ol Y. m (f+, p}, You know that 
when you use the 1/ and [..operators on angular momentum eigenstates, you 
get this: 

I.' I r'. m • = j'(,' ♦ l)*ri 
U I i rri = /liti i J 1 , fii--- 

So the following must he true; 

1 ^ LY.UmO -i(i> l)oY,„fft,^) 

^ LY^.fa. 40 = nmam-) 

In iart, you ran go further. Note that L. depends only on 0, which suggests 
that you ran split ¥!.,(#, 0) up into a part that depends on 0 and a part that 
depends an o. Splitting Y |W (6.4>J up Into parts looks like this: 

Y,^.+)-e ln C*X,C*) 

That's what makes working with spherical coordinates so helpful — you can 
split the eigen Functions up into two parts, otse that depends only on u and 
one part that depends only on 41. 
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The eigenfunctions of L z in 

spherical coordinates 

Start by limiting the cigcnEunctiuns ul L. in spherical coordinates. Ln spherical 
coordinates, the L- operator looks like Lb is.: 




So L^/0, <) - Is 

}* iff 



which Is the lollowing: 

= -^ i, 

And because t.-Y^u^) - this equation can be written In this version; 

Cancelling out terms from the two sides of this equation gives you this differ- 
ential equation; 

This looks easy to solve, and the solution is just 
*,M = C*™ 

where C is a constant of integration. 

You can determine C by insisting iliat c„,(» be normalized — that is. that the 
following bold true; 

2r 

J<b ,l>)^(^)^=< 

1 




150 
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which gives you 





So <^ lip (c) is equal to this: 

)= 

Yuij'rc making progress — -yuu've been able to dctcrullue the lurm oi 
soY,,,(e.<tO = 0 a„(&) which equals 

(2*0 ' 

That's, great — you're hallway there, but you atm have to determine the form 
of ft the eigenfunction of I.'. That's coming up next. 




The eigenfunctions of L 2 in 
spherical coordinates 

Now you're going to tackle the eigenfunction of I. . ft flI .(0). You already know 
that in spherical coordinates, the L operator looks like this: 

L ” h sto Ij t “ & do } + Sf ] 



That's quite an operator. And you know that 



Y ta ( fl . 0 )* e ,.(0 j^-r 
I,--) 



So applying the L' operator to j'. 11 (h,d> gives you the following: 






2* 



1 f 



sill t 1 nit 



m I Hi:,e DS ) ' sin W 



0 f 6t]i 



Tiqht 



m atari 
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And because L"Y lV jH,g.) = j(f +. I lfl"Y jT1 (9, tb) = * IXj^e/mtOJdVtoj, this equa- 

tion* he cum rs 



-h z 

(2JT)' 



i a 

sinfl 






^[j+i)A : eue)-f 



(2 a 1 



Wow, what have you gotten in to? Cancelling terms and subtracting the right- 
hand side from i h e Left finally gives you this differentia] equation; 



L 



sin 9 






Combining terms and dividing by e""* gives you the Following: 



Holy cow! tgn't there someone who's tried to solve this kind of differential 
equation before? Yes, there is. This equation Is a t.egpitdre d Wewttftti equa- 
tion, and the solution# are well-known. (Whew!} fn general, the solutions take 
this form: 



0i»(e) = CA n P 1 - ltf (cos0j 



where P A „(cose) is the Leipndrt function. 

So wliat are the Legendre functions? You can start by separating out the m 
dependence, which works this way with the Legendre functions: 

where P,<c) is called a Legendre polynomial and is given by the Rodrigues 
formula: 
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You can use I his equation to derive the first few Legendre polynomials like this: 

^ = 1 
* r,U) : X 

V* Fax) = Vs( 3 -t 3 - IJ 
^ P-i(jt) - V* (5jf ' - 3x) 
u" Pifx) - Yb (35.it 1 - 3UY * 3) 
t-" P,C t) - V» (63/ - TOf 1 i L&v) 

and so nn. That's what the first few P,(jr) polynomials look Sifts. So what do the 
associated Legendre functions, EY,,(jr) look like? You can also calculate them. You 
can start, off with EYiU'}, wfiern m = 0. Those are easy because F„(a) = lAl’-v), so 

1 ^ Pi,(.r) = jt 

t" l) 

^ (5t ' - 3x) 

Also, you can (tod that 

^Pi,{;0 = fl-r) :rj 
P, s {x) = 3x(l - .r} 1 -'- 
t* P a ;{jc 3 =a(l -.V .I 1 ' 3 
V* f n {■* I = | { 5 a” J - I ) 1 1 - -Y ‘ ' 

^ P :ii (;0 * JTvf] -jf) 

* P Ax) = lox(T,-r ) J ' 

ThiU equations yivc you art overview of wtiat the P; n , functions look like, 
which means yuu 're almost doin' . As you may recall, ©,*,(0^ the 0 part ol 
Y hl ,C0, $). is related to the P,- rt Junctions like this: 






And now you know whal the P Al , functions look like, hut what do C...„, the con- 
stants, look I ike? As soon as you have those, you’ll have the complete angular 
momentum eigen Functions. 1^,(0, *}. because Y.. n (0, o) = 
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You can go about calculating the constants C*, the way you always calculate 
such constants of Integration in quantum physics — you normalize the eigen- 
functions to L. For Yftju.o) «- H.„(o>i*, l (i|>), that looks like this: 

J <te] Y„ ' ( 0.4 ] Y„ { 0.4 ) sin 0 cl 0 = 1 

j b 

Substitute the following three quantities in this equation: 

^ Y,. ; (n o ) ~ ti\M) 

* *!«(*) ■ Cft,.P,„{cosu) 

You Sot the Eollowlng: 

|C f Z* S J 

LJd.f^JIVfcosfl) S\n9d8=\ 

Z.7T Q p 

The integral over i> gives 2s, so this becomes 

K , J|P Al 'fcGe[t)| < sinS<ltf=i 




You can, evaluate the integral to this: 

Ip J 2 (^ + #»)! 

1 “ 2f + l(f-^)! 



$0 in other w orris: 

(2f + l](f-m)! 









where m 2 (i 



Which means that 






ufj+mji 



l'= 



P T i cos 9 ) where raid 
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$0 Y,- ni (0, fi j = 0^,(0 which 55 thft Angular momentum. eigenfunction in 

spherical coordinates, is 



M a.p)=(-ir 



[2l4l)[l-mj! 
4jr( 1 + m } ! 



P ln ( cos 9 ]e ' lf where .7,' 2 (1 



The functions jjjven by This equation arc called I he nom\nUi€d Kpheriail hay- 
mantes. Here are what the first few normalized spherical harmonics look like:: 

* Y|,i(b, ■*} - (V-n) 1 '" COS« 

^y, 4 M=t} | 

* V B 0 , <i> - f 3 en«'-'u - n 

*- = c "sin* 

^ 1>) - C‘!-nf : e' 1 ' slirfl 



In lacr, ytnicauuse these relations to convert I he spherical harmonics to 
rcctai Lijular ctror Jinatcs: 

= — 

y 

v* $in(?sin$ = — 

cos 9 = i 

r 



Substituting these eq nations into 

fai4 i)f]-wj! 1 ■' 



= « 



An [l-tnij! 



p» ( * ) '■ s ” where rn > (1 g)ves voy , h e 



spherical harmonics In rectangular coordinates; 
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* *r, {**«)•*(%.) '«**'>’>/ 



\ -rl 15 / )'‘{* : -y'± r ^) 
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Chapter 6 




In This Chapter 

Discovering sp-Ln with she Stern-Gerlach experiment 
p Looking at eigenstates and spin notation 
p Understanding fermions anti boson;, 

Comparing the spin operators with angular mom etumn operators 
p Working wills spin L /i and PauLi matrices 



#«r | lyricists liavc suggested that orbital angular momentum is not lire 
f only kind of angular momentum present in an atom — electrons could 
also have intnttsic built-in angular momentum. This kind of built-in angular 
momentum ts called south Whether or not electrons actually split will never 
be known — they’re as close to point-like particles as you can come, without 
any apparent Internal structure. Yet the fact remains that they have Intrinsic 
angular momentum. And that's what this chapter Is about — the Intrinsic. 
huLLt-in quantum mechanical spin of wbalnmic particles. 



The Stern ~Ger(ach Experiment and 
the Case of the Missing Spot 

The Stern-Gerlach experiment unexpectedly revealed the existence of spin 
back Ip 1922. Physicists Otto Stern and Wall her Gerlach septa beam ol silver 
atoms through the poles of a magnet — whose magnetic field was in the 2 
direction — as you can see hi Figure £>1, 
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FtflHp* 5 - 1 : 
Trie: Siem- 
Oc? In r:h 
pjtpenniem. 



Si vs r s i n ms 



Spin up 
Spin dawn 



Magnet Screen 



Because 46 ui silver's 47 electrons art: arranged in a syimUt'lrical cloud, they 
cmribuie nothing to the urbital angular momentum of the atom. Hit 471 h 
electron can be In 

I k" The 5s state, in which case Us angular momentum Is / - 0 and the z com- 
ponent ol that angular momentum is i} 

y* The op stale, In which case Its angular mo men turn Ls f = 1 , which means 
that lhn s component of its angular momentum can bn— 1, (1, or 1 

That means that Stern ami Gerlach expected la sec one Or three spots tin the 
screen yOu sue at right hi Figure 6-1, corresponding to the different status of 
the i component of angular aioiiiniitjet, 

Bui famously, ihey taw only two spots. This puzzled the physics community 
for about three years. Then, in 19-25-. physicists Samuel A. Coudsmlt and 
George F„ Uhlenheck suggested that electrons contained intrinsic angular 
momentum — and (hat Intrinsic angular momentum is what gave them a mag- 
netic moment that Interacted with the magnetic field. After all, it was appar- 
ent that some angular momentum other than orbital angular momentum was- 
at work here. Ariel that built-in angular momentum came to be called spin. 

Tin* beam of silver atoms divides in two. depending -on the spin of the 47th 
electron in the atom, so there are two possible states of spin, which came to 
be known as up aud down, 



jpy 



nan too 



material 
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Spin Is a purely quantum mechanical effect, and there’s no real classical 
analog. The closest you can come is to liken spin to the spin of the Earth as 
it goes around the sun — that is. the Fourth has both spin (because it's rotat- 
ing on its axis) and orbital angular momentum (because It's revolving around 
the sun). Rut even this picture doesn't wholly explain spin in classical terms, 
because it's conceivable that you could stop the Earth from spinning. Rut you 
can's stop electrons Irurn possessing spin, and that also goes lor other sub- 
atomic particles that poasuss spin, such as prutous. 

Spin doesn't depend on spatial degrees of freedom: even if you were to have 
an electron at rest (which violates the uncertainty principle), it would still pos- 
sess spin. 

Getting Down and Dirty with 
Spin and Eigenstates 

Spin throws a bit of a curve at you. When dealing with orbital angular 
momentum (see Chapter 5). you can build angular momentum operators 
because orbital angular momentum is the product of momentum and radius. 
Bui spin is built in: there’s no momentum operator involved. So here's the 
crux: You cannot describe spin with a differential operator, as you can lor 
orbital angular momentum. 

In Chapter 5, 8 show how orbital angular momentum can he reduced to these 
rfiFFeracU iat operators: 




^L, = Yl^-ZP, = -Ur 



y-eL-z-eL 1 
l\-!z J 



u" L . - ZF - XP, = -it j 




I. 



.2 



= xp - vp = -it, 




And you ran Fin d eigenfi met ions for angular momentum, such as V-,-: 
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But because you can't express spin using (inferential operators, you can’t 
find eigenfunctions for spin as you do lor angular momentum. So that means 
that you're left with the bra and ket way of looking at things (bras and kets 
aren’t tied to any specific representation in spatial terms). 

In Chapter 5, you also take a look at things in angular momentum terms. 
uitrmLuciriglhc eigenstates ul orbital angular momentum like this: If, m> 
(where 1 is the angular mamma turn. quaulum [lumber and m is the quail turn 
number of the z component ul angular iiLOinentum). 

You can use the same notation lor spin eigenstates. As with orbital angu- 
lar momentum, you can use a total spin quantum number and a quantum 
number that Indicates the spin along the z axis {Nnie: There’s no true z axis 
built In when It comes to spin— you introduce a z axis when you apply a 
magnetic lie Id; by convention, the z axis Is taken to be in the direction of the 
applied magnetic Field). 

The- letters given sn the total spin quantum number and She z-axis component 
id tlm Spin arc s and m (you fiOmethrms see them written as s and fld h ). In Other 
words, the eigenstates oi spin are written as is, 

So wliaf possible values can s and m take? Thai’s coming up nexL 

Halves and Integers: Saying Hello 
to Fermions and Bosons 

In analogy with orbital angular momentum, you can assume that m (the z-aris 
component of flit spiuj can take the values -a. -s * I, .... s- l. and s, where s 
Is the total spin quantum number. For electrons, Stem and Gerlach observed 
two spots, so you have 2* -* 1 » 2, which means that s - (6. And therefore, m 
can he *'h or So here are the possible eigenstates, for electrons in terms 
of spin: 




I 7 -h 



I ->f :V 

So do ail subatomic particles have: s = Nope. Here? arc their options: 

I v* Fermions: hi physics, particles with hail-integer spin arc called fermions. 
They include electrons, protons, neutrons, and so on, even quarks. For 
example, electrons, protons, and neutrons have spin s - /., and delta 
particles have s ■ V 3 , 
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And tine eigenvalues of the S : operator are 
5. |s.m > = mfr.s.m > - ^ $,m := 

You can represent these two equations graphically as shown In Figure 6-2, 

-where the two spin stales have different projections along the z axis. 




Spin V2 matrices 

l iant: to take a look at the spin cigens tales, and operators lor particles of spin 
'f: in icrais of matrices. There arc only two possible stales, spin up and spin 
down, so this is easy. First, you can represent the dgcnstalc I V\ '/»> like this: 



\Yx 



1 

2 r /2 




And the eigenstate I VW -'/«> looks like this: 





Ci 

I 



Now what about spin operators like S ? The S“ operator looks like this in 
malrix terms; 
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s j 

< 



<Vr ] An4 l A> 

Vl~ 1 Li 1 2 '} 2 > 



<- 1 / l/k J 1 . I s 

^ /2'/2P /V / 2 ■ 

< 1/ _ ]/ S : ! 1/ - 1/ s 
- •>' /I 1 / 2 ’ / 2 ■ 



And this work* out to be the following: 



S : 




I ii 
(i ] 



Similarly, you can represent tlieS, operator this way: 
s. = 



1 > 



'k/'a'/a s H\K :> < Vi') 2 s hr /2 

<K-K SJ IK-K> <K'-KI SJ IK--K> 



This works out to 







(f 

-I 



Using tlae matrix version of S., for example, you can find ihs 2 component 
of Hie spin of, say. die eigenstate ->}:>. Finding l lie 2 component looks 
like 1 1 tie : 



5= 1 -/.> 

Putting this in matrix terms gives you this matrix product: 



ft 


1 0 


;i 


2 


Cl J 


] 



Here's what you get by performing the matrix multiplication: 



ft 


1 0 


n 


-h 


1: 


2 


h -1, 


i 


~Y 


1 
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And putting this hack into ket notation, you get the lolloping: 




How about the raising and lowering operators S. and S_? TticS. operator 
lu-uks like this: 

S,= 
ft ll ] 

1 V 



And the lowering operator looks like this: 

S = 
ft Jti h 
]|l 0 



$0, For example, you can figure out whal S,l '/a, -*}:> is. Here it is in matrix terms: 



li 

L 



11 1 ii 

0 0 1 



Performing the .multiplication gives you this: 



h 


o 


1 


to 


0 

z- 


I 


i 


t! 


0 


i 


] | 


0 



Or In ket form, it's S. i -*/s> = ft I Cool. 



Pauli matrices 

Sometimes, you see the operators S* S.. and S : written In terms of Pauli 
malrices. <s A . a,, and or.. Here's what the Pauli matrices look Like: 
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Dang kids tit ike Thysics 
Departiment, messin h -with 
old Hitch's flashlight again! 



■S 



The 5 th Wave 

iS * _ — — 



fly TJL^h TgRUdml 
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Chapter 7 

Rectangular Coordinates: Solving 
Problems in Three Dimensions 



h This Chapter 

imploring the Sell r&di nger equation in thex y. and * dimensions 
^ Working with Free particles in 
► Getting into rectangular potentials 
^ Seeing harmonic oscillators In 31) space 



■ jfiie-dimensiotial problems are all very well and good, but the real world 
lias three dimensions, This chapter is all about leaving one-dimensional 
potentials behind and starting to take a look at spinDess quantum mechanical 
particles in three dimensions, 



Here, yon work with three dimensions in rectangular coord .Inaras-, starting 
with a look at the SchriHluiger equation in glorious, rcai-lde ill) . You then 
dcivc into fnr particles, bos potentials, and harmonic oscillators. i\ot€: t!y 
the way, tJie nest chapter uses spherical coordinates because some problems 
are belter in One system than the Other. Problems with spherical symmetry 
are best bandied lit spherical coordinates, for example. 



The Schrddinqer Equation: Not# in 3D! 

In one dimension, the time-dependent Sdi rodinger equation, (td the type u i 
Chapters 3 and A that let you liml tlife wave function) looks like this: 

2m ib 7 * ■' xjl} + v 1 r f ] ^ xJ m &t it ^ x ' f ■ 
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And you can generalize that Into three dimensions like this-; 

«( W + W + .J - ) v 1 x ' y ' s ' 1 f j ‘ V ( r ) v (*->'■ ZJ ) * ifi |r ^ ( *->■ zl ) 



-fr 

2 m 



Using the LaplaciEm operator, you can recast this kilo a more compact form. 
Keres what El iu Liplnriafi looks like: 



( .il + J?- + JL 1 

1 rl.V 1 rJy- dZ : J 



= V 



And here's the 3D Schrorlinger equation rising the la placian: 

V : y { x, y, t,t l + V [ Jf,y , z , f ) ip ( x..y, z J ) = if/ Jj- y(x,y , zJ ) 



■him, 

2 m 



To solve this equation, break nut the time-dependent part ol the wave 
function: 



' ku¥ 



Here, y, z) is the solution of the ilmefodependent Sehrodinffer equation, 
and I-: is the energy; 

f *■ y 1 J ) +v ( *■ y* * ) H * ■>', * ) = ¥.v \ x. y ..z ) 

So far, so good. Hut now you've nm info a wall — the expression V ) 

Its in general very hard to deal with. $ o the current equation Is In general very 
hard to solve. 

So what should ynu do? Well, you ran torus on the case in which the equa- 
tion is septirnbfs — that is, whem you ran separate out the jc, y. and z depen- 
[Lence and Find the soJutinn in Each dimension separately. In oilier words, in 
separable rases, the potential, Vf.r. y, s}, is actually the Kuril ol the jc, y, and s 
potentials: 

VCr, V, Z) - vj» 4 V, ( v) 4 V.(;) 

Now you can break the Ham Ikon Ian I n v V ( v s J + v f z) v (-*. y. z ) 

= Ei^r,y,z} inlo three Haniililonizris, H p H,, and Hz 
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(H r 4. H, * JHJytx V. z} = Evt.v. y. z) 
w h ere 



H,=^ + V l.f | 

r 2 m civ* ' ■ 1 

U m ~L L * v -W 



2 in By 

-fr_^ 
1 2m Hz 






When you divide UP Hie Hamiltonian as in fli. * H, * v, z) = Fyi>', y, z), 

you can also divide up the wave function that solves. that equation. In par- 
ticular. you can break tlte wave function into three parls, one tor a*, y. and z: 

Vfa ¥. £) = X(.y)YCv)Z(z) 

That's going to make life considerably easier, because now you can break the 
Hamiltonian up Into three separate operators added together: 



2 m Bx* 



4Vj;x) 





{i 



-/l : 



m Bz : 



+ V. 



= li 




The total energy, F, Is now the sum of the x component's energy plus the y 
component’s energy plus the z component's energy;: 

F = [•: - K, t F. 

So you now have three independent Schmrtintjer equations for the three 
dimensions: 

^&F 2 < ; > +v i>> z (d=MU 
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- ft 1 r> 



u* 1m y/= 



*-Z(z)~E £{z) 



If you rewrite these equations In terms oi the wave number, k, where 
* = 1 retentions become lhe following: 



„£j*w-*.**w 

In ihb sec l ion. you take a look at the solutions to these equations, find the 
total energy, and add time dei>endenee. 

The x, tj f and z equations 

Trike n look ;s( ihe x er|u:bliim for 1 lie free particle, ~yX(je ) = -k, J X(jt) . YciLi 
4- ii3i write ils general Kolulia n ns 

Xi.vi - At- ,L " 

This la a plane wave, and normalising ll [as I discuss m Chapter 3) gives you 
th Is. 



xuh 



Tii« y mid 2 G{uu]itmttnLs EoLlnw 1 tit; smut; Jtirm: 
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Note that k* * ft/ + ft/ is the square of the magnitude of k — that is, ft-. 
Therefore, you can write the equation for the total energy as 








Note that because F, Is a constant, no matter where the particle Is pointed, 
all the eigenfunctions of ^ X(j ) = -A/X(jr}, ^rY(y) = -ft, ' Y( > j „ and 

I 

h/Z{z) are infinitely degenerate as you vary and ft.. 



Adding, time dependence and 
getting a physical solution 

You can add time dependence to the solution lor v(*, y. i). giv ing you qi{>. > J , 
z, !}, if you remember that x,y,z,l) = npr( x,y,z le "* . That equation gives 
you tills lorrn lOf ‘•f(r y. Z, S'): 

¥ {x,y.zj) = —Lre' ir *$ 

■ [2*) ! 

p 

Because &= , J , I he equation turns info 

V U_w)=^ U"'“) 



In lad, now that the right side of the equation is in terms of the radius vector 
r, you can make the left side match: 



0>.' j = 7rV"‘" 






That's the solution to the .Schrodingpr equation, hut it's unphys iral {as I 
discuss for the une-d i mens icmaL Scliradinger equation lor a Irae particle in 
Chapter 3). Why? ‘frying to nnnaaJize litis equation in three dimensions, for 
example, gives you the following, where C is a constant: 



-W 



nqnti 
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Inside the bos, say that V(*, y, z) - fi L and outside the box, say that Vfx, y, z) 
So you have the following; 



vi>-.y.-! = 



0, where 0 c x <L h . Oc y < L.. 0 ■" z < L 4 
m otherwise 



Dividing V(x. y, *) Into V.fjf), V, i>), and V.(jj) gives you 



^ V.[jf] = 

^ V :i [y| = 
l> V.fz) = 



0. where 0 < x < L„ 
■=•=■ otherwise 

0, where < y < L , 
— otherwise 

f l , where 0 <. z < I . . 
w titherwise 
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Okay, because she potential goes to infinity at the walls of the box, the wave 
function, vOf, }•, +), must go to zero at the walls, so that’s your constraint. 1 ji 
3D, the Schrodlnger equation looks like this In three dimensions: 

iy [ V, y ,x )+ V( jr .Jr*, r) V ( x, j. £ ) = Ev( Jr. t ) 



Writing this out gives you the following: 



57 , : I ^*| ? ^)r(xy,.-) 4 V(,. J ,z) v (,,y,z;i=Er(^.z) 



Take this dimension by dimension. Because the potential is separable, you 
can write si(.e. r r a) as yix. y, z) = Xi>-)Yi;yjZ(£j. Inside the box, the potential 
equals Kero, so Ihe Schr&dinger equation looks like this Jor x. y. and &■ 



^^•|,XP)-E.XP) 









The ncjit step is to rewrite these equations in terns of the wave number, k. 



Bemuse 




, ynu cats write the $rhrmlmi$er eqnatirms lor jf r y, and z as 



the following equations 



*£rX(x).-t,‘X(x) 

„ §nr)-»My) 
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Start by taking a look at the equation for x. Now you have something to work 

with — a second order differential equation. a“F x ! x ■ v ■' Here are 

the two independent solutions to this equation, where A and B are vet to be 
determ iited: 

I * X |Oi) - A slnCfer) 

^ x,(a-) - B cos(fcy) 

®p2 ^ ^ ^ 

So the general solution of af X|-V 1 ~ k> *'• " •' is the sum of the last two 
equations: 

X(.v} = A S in(A.r) t H c05<*-v) 



Great Now take a took it determining the energy levels. 



determining the energy levels 

To be able to determine Lite energy levels of a particle in a bo* potential, you 
need an exact value for X(.v) — not just one of the terms of the constants A 
and 13. You have it? use the boundary conditions to find A and 6. What are the 
bonn dare conditions'? The wave function must disappear at the boundaries 
of the box, so 

I u* X(d} = Q 

^X(L,)=b 

So the fact that vfp) = 0 tells you right away that B must be 0. because cosip) 
- I. And the fact that X(L.) - 6 tells you that X(L,) - A sin(A,L,) - 0. Because 
the sine is 0 when its argument is a multiple of n, this means that 

&U=n.?r n, "1,2,3... 




And because h ■' 



Ir 1 



it means that 




n, -1,2.3... 



£. frx* 
2mL 
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That's the energy in the x component of the wave function, corresponding to 
the quantum numbers l, 2,3. and so on. The total energy of a particle of mass 

n 'fiV 

m inside she box potential Is E = K * E, + Following R, = r , , you 



2rtiL 



fravu t his lor E v and E ; 



rr.W 

" 2mL n. = 1,2,3. 



E. = 



rt.Vi V n, =1,2.3... 



2 ml./ 

So the total enerpy of the particle Is E - E* * E - F., which equals this; 



E ■ 

n, i.- f:' 
2fnT r 

' 2 ,-nI..“ 

4 - ■ 

2ml ' 



rr, = 1,2.3... 
rr v -- 1.2.3... 
rr. = 1.2,3.,. 



And there you have the total energy- of a parlide lift the box potential. 



Normalizing the uWe function 

Now how about normalizing the wave Junction v(x y. ?)? In the a dimension, 
you have this for the wave equation; 



X(*}= A sin 



n jt x '| 

ttJ 



So the wave function is a sine wave, going to zero at x - 0 and x - L- You can 
also insist that the wave function be normalised, like this: 



•vl 

> 
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By normalizing the wave hincttom, you can solve for the unknown constant A. 
Substituting for Xi>j in the equation gives you (tie following: 





(t f?X 

K 




J™ li— “"-t 



Therdurc. i A|' Jsin^^— ] <t<c becomes 
CM solve lor A: 




which IILCEU1S VUL1 




firs At, new yon. haw the- constant A., s q you can get X(_r) : 




fi, L,2.d... 



Now get y(x. y z), You can divide the wave Function into three parts: 

iF&r- ,v. z) = xi; »:iYCtozi>) 

By analogy with Xfjr), you ran find Y(y) and Z(z): 




So y, t) equals the following: 



's[:>ry,z 




A 




n r m 1,2.3. .. 
n, = 1.2,3... 
rt, = 1,2,3... 
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Figure 7 - 3 : 
ft tarmwilc 
gsciHaSor 




The rFstnring force has the Form F'j, = -k, r x in one dimension, where ft, is 
the constant of proportionality between tin: force cm ihu particle anil the 
Location of the particle. The- potential energy cl the particle as a JmiciL-Q.fi 

rtf location x is, V/ e i = . This is aJ$o sometimes written as 

V(*)- Jmfti.V 

whpra 

In this, section, you lake a Look at the harmonic oscillator in three dimen- 
sions. In three dimensions, the potential looks like this: 



V [ X. y,z) - I mta , x + I tmajy " + 1 mot/* ‘ 



»-y 

r / rn 
h / 
m 



/ m 



eu 



m 



*_v 

f / n 



Now that you haver a iorm tor the! potential, yon can start talking In terms of 
SeTiroeiingeV.s equation : 




ai 4 ■>- + £- 

d.v- ,ty : i)z J 






x,y, z ) ->- V ( jt. y. z ) y ( x. y,z ] = E|y (x y, z } 
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Chapters 

Solving Problems in Three 
Dimensions: Spherical 
Coordinates 



In This Chapter 

Problems in spherical coordinates 
iTFer particles in spherical coordinates 
Square well potentials 
'sol topic harmonic oscillators 



» n your other II le as a sea captain-slash* pilot, you're probably pretty 
4 familiar with latitude and longitude — coordinates that basically name a 
couple of angles as measured Irom the center of the Earl h. Put together the 
angle east or west, the angle north or south, and the all-important distance 
from the center oi the Earth, and yon have a vector that gives a jjnnrl rlescri|j- 
tion oE location in Ihret dimensions. That vector is part oF a xphericaf coordi- 
nate system. 

Navigators talk more about the pair *ul angles than die distance (“Earth’s sur- 
face” is generally specific enough lor them), but quantum physicists Find both 
angles and radius length Important. Some 3D quantum physics problems even 
allow you to break down a wave function Into two parts: an angular part and a 
radial part, 

In this chapter, 3 discuss three-dimensional problems that are best handled 
using spherical coordinates. (For 30 problems that work better In rectangu- 
lar cuordmalc systems, sec Chapter 7.) 
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A New Angle: Choosing Spherical 
Coordinates Instead of Rectangular 

Say you have a 3D bos potential, and suppose that the potential well that the 
particle is trapped in looks like this, which is suited to working with rectangu- 
lar coord ijiates: 



vU.j.r) = 



fl, where Q < x < L* , 
w otherwise 



0<y <l, r , 0<z< I., 



Because you can easily break this potential down in the x. y, and z directions, 
you can break the wave function down that way. too-, as you see here; 

,Kx y z) = X(.v)YCv)ZC^ 

Solving lor the wave function gives you the following normalized result in 
rsc lingular coordinates; 



^{x,y.z) = 




H-]l L 



\n,zx) JKXy: 



stn 



sin 



I i'i ."z j 



n, = 1,2,3. : 

fj, 

n x = U.3... 



Th« eoergy Bevels also break down into separate contributions from all three 
rectangular axes: 



L = 14, * E, - t. 

And solving tor E gives you this equation (from Chapter 7): 



Em 

n, hn 

2mL „ 1 
n W 

-lr 

2mL 1 

rj. '!• 7T 

+ ~2mL/ 



n,= 1.2.3... 
n r = 1,2,3... 
n ( ■ 1,2,3... 
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But what i! ihe potential well a particle Is trapped in has spherical symmetry, 
not rectangular? For example, what If tire potential well were to look Me this, 
where r Is the radios of the panicle’s local Ion with respect to the origin and 
where o is a constant? 



VI ') = 



(?. where t) < r a 
otherwise 




Clearly, Irving lo stuif this kind of problem into a rectangular-coordinates 
kind of solution Is only asking for trouble, because although you can do it, It 
Involves lots, of sines and cosines and results In a pretty complex solution. A 
much better tactic k» to solve this kind of a problem In the natural coordinate 
system In which the potential Is expressed: spherical coordinates, 

Figure $-1 shows the spherical coordinate system along with the correspond- 
ing rectangular coordinates, ,r. y. and z. In the spherical coordinate system, 
you locale points with a radius vector named r, which has three components: 

I v* An r component (line length ol the radius vector) 
t* 9 (the angle from t axis to the the r vector) 
v* ^ {the angle from the > axis to the the r vector) 




r 



jpy 
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Taking a Goad Look at Central 
Potentials in 30 

This chapter focuses oil problems that involve central putentfah — that is, 
spherically symmetrical potentials, of the kind where V(r) - V(r). In ether 
words, the potential is Independent of the vector nature of the radius vector: 
the potential depends on only the magnitude of vector r (which is r). not on 
t he angle of r- 

When yoo wort on problems that have a. central potential, you’re able to 
separate the wave function snip a. radial part (which depends on the form of 
the potential} and an angular part, which H a spherical harmonic. Head on. 



Breaking down the Schrodinger equation 

The Schrbdinger equation looks like this in three dimensions, where a is the 
Laplacian operator (see Chapter 2 (or more on operators): 

.^frJ S V(r)^r). Ep[r) 



And the Laplacian operator looks like this in rectangular coordinates: 

: '. v J dy : ' 3z J 

In spherical coordinates, it's a little messy, hut you can simplify later. Check 
out Urn K-plierina I Laplacian operator: 



■' r rJr r Pr L 



Here, 1. is the square of the orbital angular momentum: 



L m-ft 



f 1 

sintfM du) 



i 1 
sirs 8 do 



So In spherical coordinates, the Schrddlnger equation for a central potential 
looks like this when you substitute in the terms: 



2m r f f r 



vi r )+ 



2^TL>(r) + V(,)^(^ = Fvr(r) 
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Take a look at the preceding equation. The f brst re™ actually corresponds to 
the radial kinetic energy — that Ls, the kinetic energy of the particle moving 
in the radial direction, The second term corresponds to the mtatianal kinetic 
energy. And the third term comes ptmds to The pirtentiui energy. 

So what tract you say about the solutions lu this version o! site Schrodinger 
equation? You -can note that the Birst term depends only on r, as does the 
third, and that the second term depends only on angles. So you can break 
the wave function, 'i'(r) - ijj(r. ij, ij>), into two parts: 

I t*" \ radial part 

k' A part that depends- cm The angles 
This is a Special property ul problems with central potentials. 



The angular part of \y (r, 9, <f> J 1 

When you have a central potential, what can yOu say about ihc angular part 
hi v(r r 0 . $)? The annular pari niusi be an eigenfunction of L .and as I show in 
Chapter 5, the eigenfunctions of L' are the spherical harmonics. V Al ,fU. t') (where 
I ts the total angular momentum quantum number aitd m Is the z component of 
the angular momentum's quantum number). The spherical harmonica equal 

I 

Y t , ( t hti' ) = { - 1 ) " ' !, l, | '^ y 1 ' ; P A ( cos $) o m * where m > 0 

Here are the llrst several normalised spherical harmonics; 

[4jr.ll* 

r* Y'.{^) = (^ I' — 

V r IJl [6),^)»±(^ I e siilS 

^ Y >.i s -P) = ( 

** Y„(y, i>) = ±( gl'^sihS 
^ «>} = ( ;; v ' T -)" <? ’'" S,ini0 
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That's what the angular part of the wave function Is going to be; a spherical 
harmonic. 



The radiat part of \|/ (r f 0p (f) J 

Vo l.i can give the radial pan of the wave function she name where a 
Is a quantum mum her corresponding to the Quantum state of the radial part 
of the wave function, and l is the total angular momentum quantum number. 
Tilt radial par! is symmetric with respect to angles, so it can't depend on m, 
the quantum number ui the s Cuihpuneiit ot flic angular MujuCjI turn. In other 
words, the wave function ior particles in central potentials looks like the fol- 
lowing equation in spherical coordinates: 



v( r - T o) - R, ! i(r)Y ( u.( 0 I *) 



The next step Is to solve for R„,(r) in general. Substituting yr>. et, $} 
from the preceding equal Ion into the Schrodinger equal Ion, 

= y0lS 



I . Y . | r\ : J ) 



-ft* — l— \ y R , ■; r 1 1 + Imr \ V in- tl + 

K,{r|dr-l J YjS.*) 



= (} 



W .i 




Okay, what can you make of this? First, note {from Chapter 5) that ihe spheri- 
cal harmonics are eigenfunctions of L (that’s the whole reason for using 
them), with eigenvalue !(i * IjV': 

L ! Y t ,[fto) = E[f-tl)ftX{M 

So the last term In this equation is s Imply i(i - Ijfl . That means l hat 
-A ! - — ( -T — r -^-r r R f ( r ) 1 - 'Imr ’ f V [ r } - E + ' — - 1 - U- takes the form 

-Sr r , ' r | r R rf (r | J - 2 mr ! [v(r)-F. |+/(f + i )h 2 = 0 , which equals 

■ of \ J 



.-hid: 

2m dr 



T rR »(')/ 



v (h- 



■»! / + l).r. 
2mr : 



[rR,;{ri]=E[i‘R iJ (iO 



The preceding equation is She one you use to determine Ihe radial part of the 

wavs Emctinn, H Ml (r). It 's railed the- yatfirtl equation for a nfliitraJ potential. 



T 
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The wav you usually handle this equation is to substitute n tor kr, where k ■ 
(ZmEJ ft, so that R,A>) becomes R,(4r) - RAp). This substitution means that 

r R ^(^)] + /l '^/ r [r R b( (k)]=e[ r R J> )| becomes the following 



£M£j_ £ dR.(p) 

dp 2 p dp 



]- 



ItJ + ti 



Ri(p)=o 



In this section, you see how the spherical fl easel and Neumann functions 
come to the rescue when you're dealing with free particles. 



The spherical Bessel and 
Neumann functions 

The radial part o! the equation* rf K| (/*) _ A riK| + 

tip p dp 



p' 



K,(p} = d> 



Looks tough, but the solutions lorn out to be well-known — this equation is 
called the spherical Bessel epmi/on and the solution is a combination of the 
spherical Bessel functions Lh(n>] and the spherical Neumann functions |n/Cp)3: 




Ht.fP ) = A./XP1 * MPl 

So what are the s pherical Bessel Fund ions anrl the spherical Neumann func- 
tions? The spherical Bessel lunctinns are ijiven by 




Here's what the First few iterations oE /H’p } look like: 



/ l.'M = 

/, {p} = 

i.-lp) 



an p 



an;; cosp 



3sinp Joosp ship 



How about the spherical Neumann functions? The spherical Neumann func- 
t ions are Crven lry 



n <[p)=-(-pt 



]_d_ cusp 

, P dp p 
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Here are the llrst few Iterations of m(p): 



v n, {p)-~ 

i> ",(>) = - 



ousp 



rrjsp siup 



p f P 



The timits for small and large p 

According to l lie spherical Bessel equation, the radial part of the wave func- 
tion tor a free particle looks like this: 

R'Cp) " A'j.<p) * BifrXp) 



Take a look at the spherical Bessel functions and Neumann intact ions for 
small and large pr 



n" Small p: The Bessel Functions reduce to 



i.-(v) 



Z : l Ip' 
(2'+l)! 



The Neuman it (unctions reduce to 



n '[fi )“ 



H:2 f-l) 1 p M 
2'I\ 



Larffe p: The Bessel Functions reduce to 
The Neumann fund Ions reduce to n - (p 1 






Note thaf the Neumann functions, diverge tor smaitl p. ThereEore- any wave Func- 
tion that includes the Neumann Fimctions also diverges, which is unphysicat. 

.So the Neumann functions aren't acceptable functions in the wave function. 

That means the wave Function oi r. 0, $), which eciuals K,,ii'r> Y,,jn. c- j, equals 
the following: 

■I'Cr. ft,o) Yj^fu.o) 

where ft - {2fftEJ ! "Vft, Note that because ft can take any value, (he energy 
levels are continuous. 
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Handling the Spherical 
Square Well Potential 

Take a look at a spherical square well potential ul the kind you Can See in 
Figure S-2 (J introduce square wells in Chapter 3). This potential traps particles 
inside it. Mathematically, you can express the square well, potential like this: 



V(r) = 



- V... where 0 < r < a 

I), where r > a 



Mr 



Figure 8-2: 
The 
spherical 

:s:|.i:srr': well 

potc Aiial. 




Nott that This potential is- spherically symmetric: anti varies only iri r. riot in 0 
or 4>- You're dealing with a central potential, so you can break the wave func- 
tion into an angular part and a radial part (see the earlier .section “Taking a 
Good Look at Central Potentials in 3D"). 

This section has you lake a look at (he radial equation, handling the two 
cases of 0 < t -■ a and r> a separately. 



'TIQrT 



mi 



>ri 
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Inside the square utetf: 0 <r < a 



For a spherical square welt potential, here's what live radial equation looks 
like forthe region r^o: 



dil* 
2m hr 



y [rRJ>-}]- 



V( r ) — 






| r R -( r )]=E.rR^)| 



lii tills, region, V(r) = -V L ., so you have 






2m 3r 



-V> 



2m r ' 



L' M').! 



Taking the V„ term over to the right gives you the following: 
“* 3 "[^, R ., I/ )] = ( F.+ V, ).r R,(rj 



2 m itr 



And here's wltal dividing by r gives you: 



Then, multiplying by -Zmift . you get 

f Jrfr «. ! = - !™ ( r -* V, )MO 



Now make the change of variable p - kr, where Hr - (2f«(E +V h )) ‘/ft. so 
that TC/r] becomes R,(fr) = R.f p >. UsIhr this substitution means that 

|-V fll whereO<rea , . 

takes t he lullowing form: 



VI 0 = 



Q, where r > lj 



dR : (p\ 2 dK f (p) 
dp ’ 1 p dp 



p : 



k,{p)= o 



This is the sphericaL bessel ^quallon (just as you see* For the Free particle in 
"Handling l , 'ree Particles in AD with Spherical. Coordinates"). This time, k = 
[2rtt(E +V,i)] , '"/ft, not (2rnE p '/ft. Thai makes sense, because now the particle 
is trapped in the square well, so its total energy is E + V f „ not just E. 



-W 



rianiec 



mi 



srial 
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£ f ; K,(p) + 2 dK^^ 
dp •' p dp 




''"‘‘I) 

P z 



K,(p]. = U 



The solution is a combination raf spherical Mess-cl luntdions and spherlraJ 
Neumann functions, where B, is a constant: 

R^ci} = S,|j'iCp) * ri,t,vi] 

So the radial solution outside the square well looks like this, where p - 

riZmE ) W A: 

( rA * ] = B,| j r [ p„ m ) ■+ ji, ( p,.„ t ) ' Y h {e, 0 ) 



Prum l III preceding seel ion, you know that the wave function inside the 
s nuare well is 

- A,f, (p.^ ] 



So how do you Sissd the constants A ; and B,.? You find those constants through 
conlimiity constraints: At the Inside/ outside boundary, where r = a. the wave 
function and its first derivative must be continuous. So to determine A, and 
B f , you have to solve these two equations: 



^ Y^( a, £>.*] = 



Uf 



[r-GM 




Getting the Goods on Isotropic 
Harmonic Oscillators 

This section takes a took at spherically symmetric harmonic oscillators In 
three dimensions, in one dimension, you write the harmonic oscillator poten- 
tial like this: 

V|V| = :J nua ! x- 

i if , , 

where or = — (here, k is the spring constant ; that is, the restoring force <ni the 
harmonic oscillator is F - -kx% You can tom these two equations into three- 

dimensional versions qf |he harmonic potential by replacing A' with ?: 



J py 
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VlVl' iffluV 

where <u J *^ . Because this potential is spherically symmetric, the wave 
FurinrLon is going to he of tlifis Following form: 

V r{r.^]=^.(r)Y 1 ,{M) 



where you have yet to solve for the radial function R„ii>) and where Y hl {6, *) 
describes the spherical harmonics, 



The Sobriktinger tf[[Liarlnn looks likethl$ in Utree dimensions: 



=£H 

2m Sr 



r R j‘ r )] T 



V(r] + 



f(j+L)ift J 

Ymr 



|cR„(r)]=E[riUr)] 



Substituting For V(r) Irom v[r] = | mturr gives you the fallowing: 



Tj 

2m hr 



y[r RJ r)]- 






mi' M r I 

2 2mr 



[r «,(r)] = E[fH.(r)| 



Well, the solution to this (.filiation is pretty difficult to obtain, artel you're not 
going to gain anything by going through the math (pages and pages ui it), so 
here's the solution: 



R v (r) = C.^r J exp £ -m w j . | 1 .... „ 1 ‘ ' J ( m» j ] 
where exp(aj = e‘ ; and 



' , l B l 1 


1 

2 


i . J fl^O L 




^ ' A / 


rn^!r*±iii 


: J “'1 


h ih j 


[(rr+f+l)!]' 



And the L y V) functions are the generalized. Laguerre polynomials: 



L/(r)« 



r V d' i 
ul dr'' 




J py 
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Wow, Aren't you glad you didn't .slog through Ihe math? Here ate the first few 
generalised Laguerre polynomials: 



V A 2. •'•(> ! = - i- + f) * ! 






(fe+ 2 ')l>+]) 



^ L :( O-tt 



'.b + 'A'ir' i6+2|l6+3|r (ft+l)f 0+2 )( f>+3) 
~T “ 1 — 1 -i. 



All right, you have ihu form for R, :l (r). To I ind list; complete wave lunCtiOu. 
t\ o), you multiply by Ihu spherical harmonics, Y,*(S, a): 



Now take a look at the first few wave I muttons for the isotropic harmonic 
oscillator to spherical coordinates: 



r,„. ( r, ] = -3|-( 'f - ] ' exp| -m Y„. ( S.d ) 



[ « f 

S' , {r,e. p) - I- - ( n I'"' ) 'rcsp( -mat ^ | V.,. ( 8,*) 



i f lG| = 

^ )V *«p( jY:J>^) 



As you can see. when you have a potential that depends on r , as with har- 
monic oscillators, the wave function gets pretty complex pretty fast, 

The energy of an Isotropic 3D harmonic oscillator is quantized, and you car 
derive the Following relation for the energy levels: 

E a ■ | rj + :> J/nu fr = 1,2.3. .. 



So the energy levels start at 3fito/2 and then go to Sflru/2, 7fioj/2, and so on. 
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Chapter 9 

Understanding Hydrogen Atoms 



In This Chapter 

The Sdir&dlnger equation for hydrogen 
The radial wave functions 
► Fnergv degeneracy 
t Location of Ihe electron 



ALt only \s hydrogen She most common element in the universe, but it's 
/ V also the simplest, And one thing quantum physics is good at is predict- 
ing everything about simple atoms. This chapter is all about the hydrogen 
atom and solving the Schrodlnger equation to find the energy levels of the 



hydrogen atom. Kor such a small little guy, the hydrogen atom nan whip up a 
Lot of math — and I solve that math In tins chapter. 



Using the Schrodlnger equation tcJls yon just about all yon need to know 
about the hydrogen atom, and its all based on a single assumption: that tire 
wave function must go to zero as rgues to infinity, which is what makes- solv- 
ing the Schr&dLnger equation possible. I start by Introducing the Schr&dlnger 
equation for the hydrogen atom and lake you through calculating energy 
degeneracy and figuring out how far the electron Is from the proton. 



Coming to Terms : The Schrodinger 
Equation (or the Hgdrogen Atom 

Hydrogen atoms- lire composed of a Single proton, around which rotates a 
single electron. You can see how that looks in Figure 9-1. 

Note that the proton Isn't at the exact center of the atom — the center of mass 
is at the exact center, in fact, the proton is at a radius of jr p from the enact 
center, and the electron is at a radius of r r 
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Figure £1-1: 
The hydro- 
gen atom. 




So what does theSchnddinyer equation, which will give you the wave equations 
you need, look Like? Well, it Includes terms Lor the kinetic and potential energy* of 
the proton and the electron. Here’s the term for the proton's kinetic energy: 




where V^ 1 = ^~T + . Here, x f , is the proton 1 * .v position. y,, Is the 

proton's y position. and i fI Is Sis z position. 

The Schrodlnger equation also includes a term for the electron’s kinetic energy: 

y 3 



where V ■' = 






+ ■ 



tz 2 

electron’s y position, and z, is Its z position. 



. Here. x.. Is the electron’s .v position, >■.. Is the 
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Simplifying and Splitting the 
Schrodinger Equation for Hydrogen 

Hurt's tlit usual quantum ilittliEulical SdlrociiilgCr equation lor lilt hydrogen 
atom: 



-v 

2m, 



^V(^p) 



2m . 



‘t r r | 



Tlit problem is that you're taking into account tlit distance Hit proton is 
from the center ol mass of the atom, so the math is messy, if you were to 
assume that Hie proton is stationary and that r t . = U, this equation would 
break down to the following, which is much easier to solve: 



Unfortunately, that equation isn't exact because it ignores the movement of 
the proton, so you see the more-complete version of the equation in quantum 
mechanics Lexis. 

To simplify the usual Scht6dinger equation, you switch to center-of-mass coor- 
dinates. The center of mass of the proton/electron system is at this location: 

R = !flr - * |, ' ,V - r - 



And the vector between the electron and proton Is 



’<• - r-. 



Using vectors Rand t instead rsf r„ and r. makes the Schrodinger equati™ 

easier lo soh'e, The Laplacian for R Is V„ : = +■ ^7. And the Laplacian 

For r is V 1 i r . 

* fix' iiy OZ' 



How can you relate and V/ to the usual equation’s v, J and V,’ ? After 
the algebra settles, you get 



_L T? 2 + _L v 7 
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Well, well. well. This. equal ion lias, terms tliat depend on. either tff(ti) or f(r) 
but not both. That means you can separate this equation into two equations, 
Like tLiis (where 8 Lie total energy, £, equals E* -» EJc 

" m*j v ■>(*!- 1 E - 

Mulli|)!yina gjij j,-V,V(»)*E. by y(lt) 0nx you 



And multiplying 



2miy[r) 



l 

YW> )-pr= E, by fir) gives you 






Mow you have two Schrodlnger equations. The nevt two sections show you 
how to solve th i? i st Independently. 



Solving for \j f(H) 



In V „’VI R) = E b ^(R), how do you solve loir f (R). which is the wave 

Function oi the centered mass of the etectmu/pruton system)' This is a 
s traightforward differential equation, Euldthe solution is 

H>(R) Ce * r 

Elere, C is a constant and k is the wave vector, where k ■ — * . You can 
FLnrtC by insisting that f(R) be normalized, which means that 

i = ]vr[R)p~{R)d’R 



This equation tells you that t: = - -Therefore, 

i'iv'y 

c^r- 



Wric 




Chapter 9' Understanding Hydrogen Atoms 



f 3 h 



'Wy 



In practice, K„ m sns mail that people almost always lust ignore y(R) — that is, 
they assume it to he S. In other words, the real action is in not in ip(R); ijr{K) 
is the wave hunt ion for The center of mass nl the hydrogen atom, and i^Vi is the 
wave Junction for a (fictitious) particle of mass m. 



Solving for\\f(r) 

The Schrodlnger equation lor 41(0 la the wave function Fora made-up particle of 
mass m (in practice, m -m, and y(r) is pretty close m ijn> F ), so the energy, E„ Ls 
pretty close to the electron's energy). Here's the Schrhdinger equation for qrfrj : 



You can break the solution, w(rj. into a radial ]jart and an angular part [see 
Ciiapt cr &}: 

¥*•} = Rnl(^V A „C0, 10 

The angular part of vO) Is made up of spherical harmonics. Y r „.(£J. ■!'), so that 
part’s okay, Mow you have to solve for the radial part, , Here's what the 
Schrddkiger equation becomes lor the radial part; 



z±Ld 

2m (It 






■ r ' _k R I y \-£ y R [ri — F. r R . | j- '■ 
Zmr- • r " > r f 



where r - I rl . To solve this equation, you take a look at two oases — where 
ris very small and where * Is very large. Putting them together gives you the 
rough form of the solution. 



Solving the radial Schrddinqer equation 
for small r 



For small r, the radial wave Junction must vanish, and you have 



And multiplying by 2 rnfir, you get 
-J 2 



dr 



r[ r M r l| + - ' r : -^-R„lr J-C' 



J py 
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The solution to this equation is proportional to 

K UJ - A/ h Fir 

Note, however, that R lt ,{n must vankgh as ' sloes to zero — hut the r " ' 1 Term 
goes to infinity. And that means that Ft must he zero, so you have this solution 
Fur small r. 



ITTr.l - r 



That tnket care of small r, The next section takes a look at very large r 



Solving the radial Schrodinger equation 
(or large r 

j 2 1 

For very large r, V r : if/ ( r )- f \ yr | r | = E r iff [ r i h-eco me* 

|i.[,Kj r ]] + 2ME r H rf (,) = 0 



Because the electron is in a bonnet state in the hydrogen atom. E -= (l; thus, 
[lie solution to the preceding equation is proportional to 



R .,<0 - Ae'*' * Bt> ,r 



where X = 



{-ZmF.y 



Note that R^Vl - A# ' •* lit' ' diverges as r goes to iiilinity because of Hlr Be 1 " 
term, so B must be equal to zero, that means thMR,„4y) - e~ v . In the .next sec- 
tion. you put the solutions for small rand large r together. 



you got the poWer: Putting together 
the solution for the radial equation 

Putting together She solutions lor small rand Large r {see the preceding sec- 
tions), the Schrodinger equal con gives you a solution to the radial Sc FirncLmger 
equation ot K,..(r) = j f{r)e ', where E{r) is some as-yel-ucicFeternunecI funcTinm 
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This equation gives the recurrence relation of the Infinite series, 

jr fei 0 ,r M: + sf ttj?- - A ( A + f + ] ) Vr ** = Q, That is , If you have 

one eociJicLenl. you can get the neat one using this equation. What does that 
buy you? Well, take a look at the ratio ol 




fl*_, t(* + 2 f + ]) 



Here's what this ratio approaches as k goes to ■< 




n 

k 



This resembles the expansion lore 1 *, ivhich is 




As for et. the ratio of successive terms is 

l)E 

* {a*r 

2 

And in the Limit A -* «, the expansion for e' approaches 




That's the rase Inr <r\ Knr E(j), ynu have 

lini — — * ~ 
i -a,., k 

Comparing those two equations, it's apparent that 
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For this series to terminate, n^,. a^, and so on must all bo zero. The 
recurrence relation for the coefficients a* is 

k(k + 2I + l)a t = 2 + tr^, 

For o Nl ] to be zero,, the factor multiplying o 4 ., must be zero for k - N * 1 , 
which means that 

A|ft + f)- -': v «0 

Substituting In * = ISi i I gives you 2 A(N +* + t )- n ^~ = I) , And dividing fey 
2 t'ivrjE you jl(N + f + 1) ™ ■ 0 . Making the Substitution N - / 4 ] — >■ n. where 
n is called the principal quantum number, gives you 

■mt- ™' =d n =1,2.3 ,. 

This is the quant bat ion condition that must he met if the series For F(r) is to 
be Finite, which It most be, physically: 

E {r) = i°y 



, (-2/nFA 1 > 

Because — . the equation nX-* jp =0 puts constraints 

on the allowable values of ihe energy. 



Finding the allowed energies 
of the hydrogen atom 

The quantisation condition for ifi{r} to remain finite as r goes in imfimry is 

n X~m^=U n=\.2M... 
tr 

where } ^ . Substituting A in!u the quajitiiatiocr-condillu]] equalum 

gives you the following: 
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Getting the form of the mdial solution 
of the Schroetinger equation 

In this- seel ton, you complete the calculation of the wave Eductions. Go to the 
calculation oE K,„(r) (sue the earlier section titled “You got the power: Putting 
together the solution tor the radial equation”}. So far, you know tlul 

R, s (r) = r "f (r)e * , where lfn = ^o,r", Therefore, 

‘*=•3 

&J.r) = r'e 



lu lact, this isn't quite enough; the preceding equation conics From solving 
the radial Schradinger equal ion: 



t,« £■ r R - Id 1 l )£p r S - 1>>- 7" r MO- E.r 



The solution is only good in a multiplicative constant, so you add such a con- 
stant, A,, r [which turns out to depend on die principal quantum number n and 
the angular momentum quantum number 1), HKe this: 

M r )= 



You find A*,, by normalising Tt„4'r), 

Now try to solve for R,Xr) by just flat-out doing the math, For example, try to 
find R „,(/). in this case,, n - 1 and / = i). Then, because N + l = n, you have 
IV = n - f- l. Sd N = 0 here. That makes R,„[r) look like this: 

Rjr^Ajre '% aj * 

* -i 



Anri the summation in this equation is equal to 



so 



H,.fr) = A^-'e-^a, 

And because i = ft. * f = l . so r im (j-) - A,. ; e r tt„, where 
you can also write R |,i(d * A- li f}' >t a,,as 




. Therefore, 



R„4 r ) = A»W 



-r 

ttr.. 
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where n. is the Bohr radius. To find A lM acid a H . you normalize v.mfr * o) To 
!, which means integrating i|r K . ( ,(i; 6, 4)} I ~d rover all space and setting the 
rusuJJ to 1. 



Now dr - rslno dr da eft, and Integrating the spherical harm on Its, such as 
Y„_, over a complete sphere, j| Y^| sin Bd&dip, gives you 1. Therefore, you ‘re 
left with the radial part to nonnaiiz#: 

l = ?r R lt (r)|> 

ti 

Pluming K (j ,(r | = A, s ckp — a s into 1 = J r ft.„l/)| dr gives you 

V rif '. ‘ 

f = V<vj r ^p( ^ J-fr 



You can. solve this kind ol integral with the following relation: 
f Jr'iap ] ax )dx ~ 



With this relation, the equation ] - A^a f r ; exp 

ri 

f _ \ % 

l = A I( Y?, [ r espl — ^ dr = h^tL. ^7- 

1 1 J 4 



-lr 



dr becomes 



Therefore. 




A 



c.',- 




A m o„ 




This is a lairlv simple result. Because A ,, is just there to normalize the result. 

r 1 

you can set Am to I (this wouldn't be the case if = I involved 

multiple terms). Therefore, *r„ « . That's tine, and ft makes K Jtt (r), which is 



»py rights 

r J O ■ 
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Vcu know that vf n% Jr, 0, o) = ^.(r) Y.Jp. nj. 



And so ijri,,//, S, 4') becomes 



V, 



,M-p) 






Whew. In general, here's what the wave function i|r nAlI (, j: 0- looks Like Eor 
hydro^m: 







w h ere l,,-,-i' f (2rfn Ni) Ls a. generalised I aRuerre po lynom lal. Here are the f I rat 
Few generalized laguerre polynomials: 



* UV) = i 



v 

is* 



L| i/j -?• + i.i • I 
L/{r) = ^-l>+2)r- 



l> + 2)(*+l) 
2 







l6+;i)r : {b+2)(o+3}r ( ii| ^ + ji ji fr+ :i ; 

2 4 " 2 - 4 " 



Some hydrogen Wwfe (motions 

So what do the hydrogen wave functions took like? in the preceding section, 
you find that tjfnjn(r e.+) looks like this: 

V 



Here an? same other hydrogen wave h motions: 

^ *■*} = -rr e 1 ' ■ ■ v - f ) 

Vr-S 

* Vn. i r -&’ P) " -i- « ;r X* ( *■*) 
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V,Jr,6,$) = — L - (fl.fr) 

^ JJ - 



Note that i|i R!jn (r; 6, if.) behaves Hike r 1 fnr small rami therefore y.oes to zero. 
And lor large r, y nl „{r; 3. <p) decays exponentially In zero. Sn you've solved the 
problem you had earlier of tine wave Emotion diverging as r becomes large — 
arid all because dF the quantization condition, which cut the expression for 
f(r) from an exponent to a polynomial ol limited order. Not bad. 



You can see the radial wave function Ri.iCy) In Figure 9-2. R : «(r) appears in 
Figure 9-3. And you can see Rv(/) In Figure M. 



friliire 9 7; 

Tiie radial 
wave 

function — 

Fyn. 
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Figure 9 1: 

FUri 




<yn 




Calculating the Energy degeneracy 
of the Hydrogen Atom 

Fjch quantum 4 tale of the hydrogen atom 5s specified with three quantum 
numbers: ji (ihe principal quantum number), t (the angular momentum 
quantum number of the electron), and m (the z component of the electron's 



'riontec 
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Quantum states; Adding a (title spin 

You may be asking yourself — what about the spin of the electron? Right you 
ate! The spin of the electron does provide additional quantum states, lip to 
now in this section, you've been treating the wave Function ol the hydrogen 
iiltDLTi as a product of radial and angular parts: 



VniJr e. t J = K„Er)V, 11I (0. t) 



N'ow yOu can add a, spiia part, LOrTuSpOuding to the spin of the elu£ truii, 
where i is tile spin of the electron and m< is the t component of the spin: 

|vtt.) 



Tiie spi]] part ol the! equation can take the JuLluwing values: 

I r* i 'h. 7»> 
t* I "/s. 

Hence, n, 40 now becomes «. ♦): 

v i r M ) = R *■ ( r ) Y„ ( 0.6 )|s, m, } 



And this wave function can take two different forms, depending on m v 
Like this: 



f » m ,(r J .A») = R„(r]V„(* 1 ,)|l-l) 



hi fact, you can use the spin notatioji (which you use in Cliapterfi), where 




im 



For example, for l V->. you can write line wave function as 

W Jn j(M l p)=R ( .(r}Y te (fl.^) | l | ; 

_ 

0 

And tor I you can write the wave function as 

u 

What does this do to the energy degeneracy? ]| you include the spin of the 
electron, there are two spin Matei for every state In. I, n?>, so the degeneracy 
becomes 

6-a 

Degeneracy = - 1 J = 2n’ 

r-p 

So If you include the electron's spin, the energy degeneracy ol the hydrogen 
atom is hr . 

In lant, you ran even arid the spin ol the proton to the wave function 
{although people don't usually do that, because tbs proton's spin interacts 
only weakly wilh Kiayjietic iields applied to the hydrogen atom), [n that case, 
you have a wave function that Looks like the following: 

f *'•«'< *) ■ R - ( ? } Y *. & f j| -v, • rrt„ }| s,, m „ ) 

where s.. is the spin of the electron, m w as the z component at the electrons spin, 
s L , is the spin of the proton, and nr ,, is the z component ot the proton's spin. 
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If you Include the proton's &pkn, the wave function can now take four different 
forma, depending on rcu, like this: 






rij 



1 _ i ■■ 

r 2! 



i _ i ■ 

r i: 



L 1 \ 

2 ' 2 



^ a -l ('A*) = R„ [r ) V b {#.$ ) i J 



\ _ j'- 
V 1 1 



The degeneracy must. now include the proton’s spin, so that's a factor of four 
l or each I rr r i iti ■> : 



Degeneracy = ^4(2f+ l) 

II 

= 4 /] 8 



On the tines: Getting the orbitals 

When you study heated hydrogen in spectroscopy, you get a spectrum con- 
sisting of various lines, named the s (for sharp}, p (for prnteipaf). d (for dif- 
fuse), and f (for fundamental) tines. And other, unnamed tines are present as 
well — theg. h, and so on, 

The s, p, fi, f, and the rest of the tines turn out to correspond to different 
angular momentum states of the cJeclron, railed orbitals. The s stale corre- 
sponds to r = 0 ; the p state, to l = I; the d state, to i = 2 ; the /slain, to t = 3 ; and 
so on. Each of these angular momentum stales has a tliFEerenlty shaped elec- 
tron cloud around the proton — that is, a different orbital. 

Three quantum numbers — rr, i. and rn — determine orbitals. For example, 
the electron cloud for the i L, 0, 0 > stale (]& with m - 0) appears lit Figure $Wk 



-VV 
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Fifliirt % 4 \ 
The- |l r 0, D> 
state. 




TJlt 1 4 . 3 , 2 * State ( 4 f, wltJi r/i = 2 ) appears ir'L Figure 5 - 6 . 



Figure 9-E: 
ThE |4 r ^ 2 > 
stale. 
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And because 4 1 }, this equation becomes the 

following: 

r |i - T 

I K , [ r | Y . sin B dB [dd r dr 

4 4 

The p ret (Jiog equation is equal Id 

j" 2 j 

K -4 r 1 r j | Vs. [ sin e ^ J 

□ 3 

q r R „ ( ' )| #■ '-dr\ Y„ ' [ 0,*)Y m { Q.6 } sin 0 jf dfi 

i» j 

Spherical harmonica are normalised, so Shis just becomes 
I R„ilV)l V tfr 

Okay, that's The probability that the electron Ls inside the spherical shell 
From rfu r- dr. So tire expectation value of r, which is - r--, is 

(r) = ljr|R,(r)| dr 

i 

which L? 

(f) = h ! |R.,(ri] , f ! * 

si 

This is where things gel more complex, because r) involves the lagtieri? 
p-oL>-iioiikiaJs_ Bui after a lot of math, litre's what you get: 

( r ) = jr 1 ^(f)r^ = [ 3 ^- 4 ^l)]y 

t.£ . . 

where r is the Bohr radius: r, = — - . The Bohr radius is about 5.2& x 10 

me 

meters, so the expectation value of the electron's distance from the proton is 
*r> = t 1) 1(2.65 x 10 ") meters 
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Chapter 10 

Handling Many Identical Particles 



hi This Chap ter 

Looking al wave functions and I lamillonLans in many-parLicle systems 
Working with identical and distinguishable partl-cfes 
Identifying and creating symmetric and antisymmetric wave functions 
* Explaining electron shells and the periodic table 



^'Jydrogen atoms (see Chapter 9) involve only a proton and an electron, 

V I but all other atoms Involve more electrons than that. So how do you 
deal with nniltl pi e-electron atoms? For that matter, how do you deal with 
mum-particle systems, such as even a simple gas? 

In general, ymi c.nn t deal with problems like EhLs — exactly, anyway, Imagine 
the complexity til atst Two electrons moving in a helium atom — you'd have 
to take into account the interaction at the electrons not only with the nucleus 
ol the atom bur also with each other — and that depends on tiielr relative 
positions. Sonot only does the HajuiltUiaa have a term in Jor the poten- 
tial energy of the first electron and 7* tor the second electron, but It also has 

a term In “jjr 1 ^Tj for the potential energy that comes from the Interaction of 

the two electrons, And that makes an exact wave function just about impos- 
sible to Find. 

However, even without finding exact wave functions, yon can still tlo a sur- 
prising amount with multi-particle systems, such as deriving the PtiuH exclu- 
sion principle — which says, among other tilings, that no two electrons can 
be in the exact same quantum state. Set fact. you‘ll probably he surprised at 
how much you can actually say about multi-particle systems using quantum 
mechanics. This chapter starts with an introduction lo many-particle systems 
and goes on to discuss Identical particles, symmetry (and anti-symmetry), 
attd electron shells, 
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Many- Particle Systems, 

Generally Speaking 

You can ici: a multi-particle system in Figure ILf-1, where a number Of par- 
ticles arc idea titled by their position (iijnore spin ior the moment). This sec- 
tion explains how lo describe that system its quantum physics terms. 



Figure 10-1; 
A rmulti- 
particle 
system. 



Particle t 




Considering Wove functions 
and Hamiltonians 

Rf's|in bv working with the wave Function. The State oi a system with many par- 
tides, as shown in Figure 10-1. Is Riven by i|<»v T', «*,. And here's the proba- 
bility that particle 8 is in tft u particle 2 is in riV-. particle $ is in rfV : „ and sc on: 

wfr.r.r-.) d'r,d 'r,d V,... 

The normalisation ol w(r-. r., ...) demands that 

J y/j t\ X-. . r, , , )| rfVjdf-VjtfV,. ., = 1 



TlClh; 
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Figure 1 0-2: 
Amulti- 
clcclran 
alum. 



Beclntn r 




11 yens have L electrons, the wave function looks like ui{ , r., ,,,, r?. K). And the 
kinetic energy ol the electrons anti the nucleus Looks like this: 



And the potential energy ol I he system looks like HiLh: 

PE = -1 - -r, . , R ) ■ + 1 r^-M r„r f .. . -r, . It ) 

■ i k - " i-- 1 1 r 



So adding the two preceding equations, here's what you get for the total 
energy (t = &E * PE) of a multi-particle atom: 

E *(ii,r*,~ -r f , «) = I f^V, s p (r M r, t . ..r rh JT)- ^ V ,>f r, , r, . . . . r, . R ) 

~ I ~r^r\ v( r> -r„. . . r, , R) + £7-^— v(r , . r (( . . .r , , R ) 



r -R 



r - r 



Okay, now that looks like a proper mess. Want to win the Nobel prUe in phys- 
ics? Just come up with the general solution to the preceding equation. As is 
always the case when you have a multi-particle system in which the particles 
Interact with each other, you can't split this equation Into a system of N inde- 
pendent equations. 



>pyrit 
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In it a sp.fi where the N particles of a multi-particle system don't interact with 
each other, where you can disconnect the ftchrndingcr equalinn into a set ot 
N independent erpialinnv solutions may be possible. But when the particles 
interact and the Schrndinger equation depends on, those interactions, yon 
can't solve that equation for any significant number ot particles. 

However, that doesn't mean all is lost by any means. You can still say plenty 
about equations like this one if you’re clever — and It all starts with an exam- 
ination of the symmetry of the situation, which I discuss next. 

A Super- Pou/erfut Toot: 

Interchange Symmetry 

Even though finding general solutions lor equations like the one lor the total 
energy ol a mult i-part idc atom (in the preceding sect ion) is impossible, you 
can still see what happens when yen exchange particles with each other — 
and the results are very revealing. This section covers the Idea of Interchange 
symm etry. 




Order matters: Swapping particles 
With the exchange operator 

Von Cart determine what happens to I he wave hulcliuii when you swap two 
particles. Whether the wave function is symmetric under such operations 
gives you insight into whether two particles can occupy the same quantum 
slate. This section discusses swapping particles and looking at symmetric 
and antisymmetric Inactions. 

Take a look at the general wave function for N particles: 

r,, .... r, .... r ( „ Ad 

Mofff In this chapter. 1 talk about symmetry in terms of the location coor- 
dinate. r. to keep things simple, but you can also consider other quantities, 
such as spin, velocity, and so on. That wouldn’t make this discussion any 
different, because you can wrap all ol a particle's quantum measurements — 
location, velocity, speed, and so on — into a single quantum state, which you 
can call Doing so would make the general wave function for N particles into 
this: 

|s "• B ‘l ^|fj "* B I '""I ). Fhn as i said, this section just considers the 

wave ftmet Ion W(r,. a .... r„ a) to keep things simple. 
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Now imagine that you have an exchange operator. P v . that exchanges, pap 
tides i and / In other words. 

P ; ,' l K |r h r :- r. — V ...... r x) - k'( r n r, r„ ..„ r, .... r % ) 

And P. v = P,„ so 

fVfOi. d P r, r-i r, r, t\) 

= PM r i- r .... r s ) 

Also. noLe that applying the exchange operator twice just puls tire two 
exchanged particles back where they were originally, so P„- - l. Here's what 
that looks like: 



Py I it Wt r i p ^"71 .-pi Cj. , 1 - P.ptrtr - !, r-j., ■,,, c J, 1 1 - . 

- r^, ■ - . r. p.pt e. 



However, in general. P u and F J t ,., (where ij* tai) dti not commute. Thai is. 

P v P.. ni * F ia P, (ij a tm). Therefore, [3 :, , ;i . P>,„] * 0 (ij * ini). For example, say you 
have lour particles whose wave (unc Lion is 

‘ F « 



Apply the exchange operators PA. and P M to, see whether P L! P, , equals 
? u P L? p Here’s P M Tp(r lt r ? . r la e,): 

P M ^V|.r,r : .r 1 )= r ; J e’= 

-I 

And here's wltat P| ? P H i|j(r L , r s , r-.. r,) looks tike: 

P,P 

*1 



Okay. Now take a look at P 1( P s: v(r h r : . »■.. r^. I lore’s P 1S shyv r->. r-i, rO: 
P.i v( F i . r . r..r, J « ' 



And here's what P„ Py, ofr,, r-. r iq) looks like: 
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As you can sue by comparing P- ? P :J *'n J - ^ And this East eniration. 

r . 

P, a P H W(V|. r--. r, rj * l J Ui ] J |., ty, r ;l . r.,)- In other words, the order in 
which you apply exchange operators mailers. 



Classifying symmetric and antisymmetric 

ufat/e functions 

P V J - ] (see the preceding section), so note that II a wave function is an 
eigenfunction of P.. then the possible eigenvectors are 1 and -I. That Is, 
for ^( 1 % tv, „„ ry ..., r, .... eigen function ol F„ Looks like 

I ” , mOi . ri, ..„ r, .... r r r\) = iH> t , r*, ...r. .... r r .... j\) 

or — ipO, , r.'u r, -. r. #*) 

That aneail£ there are two kiclds of CLyunlunctiOiiE of tilt- exchange Operator: 

I *" Symmetric eigenfunctions: P^|i,(r h r,, r „ ..., r, r,"i ~ 

V,( r i. o r. .... r„ r N > 

Antisymmetric elgenfiinoMoiu: P.ytr, r., r „ .... »■ r, 4 ) ■ 

-YJ>|. ■■■. r .. r.., J-j. r,,> 

Now take a look at some symmetric atul some antisymmetric eigenfunctions. 
How about this one — Is it symmetric or antisymmetric? 

't'lO'i- rj = (ri - r.y 

You can apply the exchange operator P 
p i*Vi(r„ r j ... r>:, - r , ) J 

Note that because (r, - r.] J - { r, - r,)\ i|j r.} is a symmetric wave function; 
that’s because F,, w : {r 5 , #•..) _ r.}. 



How about this wave function? 




Again, apply the exchange operator, P !; .; 
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r. J -Kr- ! 



i/^rj 



r. Y r ' r +r 1 



Okay, but because n- ■ ■ rr . you know that r\ j o ; (r|. rj) - Vj(f|, r ; ), 

- r ' i I, c - r s= J 

fj) is symmetric. 



Here's another one: 



^ 3 {r,-r.)= 






New apply F |-.: 



p , ! ^( JF r fJ ! ) = 



%-h) 

[WI 



How dees that equation compare to the original one? Well, 

Sfr.-r,) 5[r : - r, ] 

_ r . y [r-r)' •• soP ' - w ! ^ ! ’ r -’* “ _v ’( r i 1 r :>’ 1 before, , r.) is 

antisymmetric. 



What a bout (his otic ' 



<M*^ = — • + r, : 4 ry 
' (M) 



To find out, apply P, t : 



All right — how'i this compare with the original equation? 



r.r, t , rr, ■ , , 

■.i T p i i , ■ # z TJ i 



Ir,-r, I 



t r - r '}' 



Okay — ^[V , r L> ) is symmetric. 



You may think have this process down pretty well, but what about this next 
wave function? 
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.■ -i rjr ± 2 i 

V*1 I" ■■ ,+r;-r{ 

I'-'V) 

Start by applying P| S : 

P . I ^[' , ..'* I ) = 7 i ± V + V -'-, 1 



.So how da these- two equations romp-are? 



r/. 






■ + r,- r : 



r. r 






7-i -r.‘-r: 



That is, iprvOi- r.) Is neither symmetric nor antisymmetric-, in other words, 
if-Xir,, r.y is not an eigenfunction of thePj 3 exchange operator. 

Floating Cars: Tackling Systems of 
Many distinguishable Particles 

All right, I If you' ve- been reading this chapter from the start, you pretty much have 
the Idea of swapping particles down. Now you look at systems of particles that you 
can distinguish — that is. systems of identifiable different particles. As you see in 
this section, you can decouple such systems Into linearly Independent equations. 

Suppose you have a system of many different types of cars floating around in 
space. You can dislinguLih all those cars because they're ah d ill era it — I hey 
have different masses, for one thing. 

Now say Hurt each car interacts with its own potential — that is. Hie potential 
that any one car sees doesn't depend on any other car. That means that the 
potential lor all cars is just the sunt of ihe individual potentials each car sees, 
which looks like this-, assuming you have N cars: 

?E=V(r |1 ^,,r,)=iv(r i ) 



Being able lo cut the potential energy up into a sum ol independent terms 
Uke this makes life a let easier. Here's what the Hamiltonian looks tike: 
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Figure 10-3: 

An electron 
coHding 
With anelhcr 
Sleclrofi. 




But now look -at the scenario in Figure lb-4 — §tiE> electrons cnuLd've bounced 
Like that, nor Jike the hounr.e shown in Figum 1(1-3. And you'd never know it. 



Figure KM: 
electron 
COllilmsi 
with another 
sleclnon. 




So which electron is which? From the experimenter's point ot view, you can't 
tell. You can place delectors lo calcli the electrons, but you can't determine 
which of the incoming electrons coded up in which detector, because of I he 
two possible scenarios in Figures ]i>3 and 10-4 




Quantum mechanically, identical particles don’t retain their Individuality in 
terms of arty measurable, observable quantity. You lose the Individuality of 
Identical particles as soon as you mix them with similar particles. This Idea 
holds true For any N-particle system, As soon as you let N Identical particles 
Interact, yon can't say which exact one is at r, or r. or r :i or n and so on. 
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Symmetry md ami symmetry 

In practical terms, the loss of Individuality among identical particles means 
that the probability density remains unchanged when you exchange particles- 
For example., if you were to exchange electron 80.281 with electron &9.830, 
you’d still have the same probability that an electron would occupy .. r r .. ... 
ami d*r riy 

Hkires what ihix idea Looks like laatbciEiaticalJy fraud s art tin; location and 
spins o i ihu part it tes): 

I ip(r,f, r.S-, ^ r,S t . rJ;,. .... r>,%) I" - u(r|S h r-4 : „ .... rfi f r,S L . .... r^JI" 

The preceding equation means that 

'P(fjSi, r >•■■■; r*. rfo .... -rvSrs'j - + 'l f ( r ,Sh r A; •#, r#,, .... 

So tlie wave function ul a system oJ N identical particles must be either sym- 
metric or antisymmetric when you exchange two particles. Spin turns out to 
be the- deciding (actor: 

u* Antisymmetric wove function's II the particles hive lialF-odd-ultegral 
spin ( /.. V-. and so on), then Hits Is how the wave function looks under 
exchange- of pari teles; 

'HCP]£i. r^i-1 -F r rS|. r ,-' ! .|. ■■■' r \S*) - r >T;" rfy rs , .... r s S % ) 

v* Symmetric wave functions IF the- particles have integral spin (ft. 1, and 
?o on), this is how the wave fund Ion looks under exchange of particles: 

^{riT- ^-Tj. .r,, r,«. ..., *-%H) = V(>|.r„ r, n <\*s) 

Having symmetric nr antisymmetric wave functions leads to some diiierent 
physical behavior, depending on whether the wave function is symmetric or 
antisymmetric. 

In particular, particles with Integral spin, such as photons or pi mesons, are 
<lf O called bosons And particles with haIF-odd-intcgr.il spin, such as electrons, pro- 
)| IJjl tons, end neutrons, ere called fermions. The behavior oi systems qI Fermions is 

tm very different From the behavior oF systems oF bosons.. 

Exchange degeneracy: The steady 
Hamiltonian 

The liamiLtonlan, which you can represent like tills 
H{r,S|, r.j-, rfi, .... r v 5\) 
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tint what il 'lie particle you' re :>Si Lrtying is, a composite particle? What if, f-n-r 
example, you have an alpha panicle, which is marie up «f two protons anri 
two neutrons? is that a fermion nr a hosnn? 



hi fact, protons aEid neutrons themselves art* made upul three quarks. and pi 
mesons arc made up ol two — and quarks have spin '/a. 



Composites can be either fermions or bosons — it all depends on whether 
the spin ol the composite particle ends up being half-odd-lntegral or integral. 
If the composite particle's spin is ! h, J /-.. and so on. then the composite 

particle Is a fermion. If the composite panicle's spin is 0, 1, 2, and so on. then 
the composite particle is a boson. 




In general, if the composite particle Ls made up of an odd number of fermions, 
then it's a Fermion. Otherwise, it s a boson. So for example, because quarks 
are fermions and because nucleons such as protons and neutrons are made 
up of three quarks, those nucleons end up being Jcniuons. But because pi 
mesons are made up of two quarks, they end up being bosons. The alpha 
particle, which consists of two protons and two neutrons, Is a boson. You can 
even consider whole atoms to he composite particles. For example, consider 
the hydrogen atom: That atom Is made up of a proton (a fermion) and an elec- 
tron (another fermion), so that’s two fermions. And that makes the hydrogen 
atom a boson. 



Building Symmetric and Antisymmetric 
Wat/e Functions 



Many of the wave functions that are solutions to physical setups like the 
square well aren’t Inherently symmetric or antisymmetric: they're simply 
asymmetric, ftn other words, they have no definite symmetry, So how do you 
end op with symmetric or antisymmetric wave functions? 



The answer is that you have to create them yourself, and you do that by 
adding: together asymmetric wave Junctions. For example, say that you have 
an asymmetric wave (unction of two particles, es.). 




To create a symmetric wave f miction, add together r-i ; ) and the version 

where the two particles are swapped, yfos*. r,s,). Assuming that y(V|S|, ryfj) 
and qi(rjS £ , r,s,) are normalized, you can create a symmetric wave function 
using these two wave functions this way— just by adding the wave functions: 

(r,s t , rA ] = 4r[ M r,s s , r d s, )+ v[ rji t . r,r, ) J 
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And the total energy Is just the sum of the energies oi the individual particles: 



b=£e 



But now look at the wave function for the system. Earlier in the chapter fsee 
-Floating Cars: Tackling Systems of Many Distinguishable Particles*'), you 
consider the wave function of a system of N distinguishable particles and 
come up with the product of ah the individual wave functions; 

^w. O =11^ < r .) 

1-1 



However, that equation doesn't work with identical particles because yon can't 
say that particle 1 is in state Vi( r iX particle S is iri state w ; (ri), and soon — 
they're identical particles here, not distinguishable particles as bduro. 

The oilier reason this equation doesn't work here is that it has no inherent 
symmetry — and systems of N identical particles must have a definite sym- 
metry . So instead of simply multiplying the wave functions, you have to be a 
llttEe more careful. 



Wa0e functions of tWo -particle systems 

How do you create symmetric and antisymmetric wave (unctions for a two- 
particle: system'? Start with the single-particle wave functions (see the earlier 
section '“Building Symmetric and Antisymmetric Wave liiin:i:ms" j: 

1 ^ V. r>,S , « -s ) = v ( V, ■ 'A )+v( r s , ',*1 )] 

^ ¥. ( r,t, , r s 5, ) » -|L[ rs , .r.s . ) - ( r,S, , r.s }] 



By analogy, here's the symmetric wave Function, this Time made up of two 
Hin^le-particJe wave functions: 

V, (Ml • ty. )=4^[ % ( r.-s, ] ip,, [ ns 4 ) + ly, (r, j , ) y c { r,s, )] 
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And here's the antisymmetric wave function, made up of the two single-parti- 
cle wave functions: 

V. ir,s t ,r 3 s s ) i-Uv.. |r I i l )^|^ i ) - YA )*,, l] 
v* 



where ft, stands for all the quantum numbers of the ft h particle. 

Note in particular that qr.Xns,. rs.'} - ft when ft, = ft,: in other words, the anti- 
symmetric wave Function vanishes when the two particles have the same set 
ul quantum numbers — that is, when they're i.11 the same quantum slate. That 
idea lias important physical ramifications. 

You can also write v,(r,$,. r like this, witere P is Lhe permutation operator, 
which takes the permutation of its argument: 

v. (Vr^.A \ = -J-&V. l/,h ) IP, {r,s, ] 

-gi! f 



And also note that you can write qij iqSi, rss.) Dike this: 

V, I r i s i .'A, )=^rl(- 1 ) Pl ^ ( *>*■ } V,, f r.A ) 



where the term {-!)' is f for even permutations (where you exchange both r,St 
and and also ft, and n , ) and -1 For odd permutations (where you exchange 
i‘ s Si and r f s-j but not n, and n,: or you exchange n, and ru Irut not c,s, and css-,). 



In fact, people sometimes write ) in determinant Eorm tike this: 



ty..fr 1 s,.rs.) = -i 7 dct 



YU*-,!,) 



( r > : ') 
¥\ (v.) 



Note that this determinant is zero If ftj *• »<,. 



Wa0e functions of three -or -more -particle 
systems 



Now you get to put together the wave Function of a system of three particles 
From singLe-partide wove functions. 
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The Pauli axcfttsian prinripte suites that nn two electrons can occupy the 
same quantum state inside a single atom. Anri that result is important for the 
structure ol atoms. Jnsteari of jusl piling on ivilly-niJ!y, electrons have to fill 
quantum states that aren't already taken. The same isn't true lor bosons — Eor 
example, ii you have a heap ol alpha particles (bosons), they can all be in the 
same quantum state. Not so tor lermions. 



There are various quantum numbers that electrons can take in an atom — n 
(the energy), ((the angular momentum), m (the z component of the angular 
momen Hum), and m, (the z component of spin). And using that Information, 
you can construct the electron structure of atoms. 



Figuring out the Periodic Table 




One of the biggest successes of the SchrSdiiiger equation, together with the 
Pauli exclusion principle (see the preceding section), is explaining the elec- 
tron structure of atoms. 

The electrons in an atom have a shell structure, and (hey fill that structure 
based on the S’ mill exclusion principle, Which maintains that no two electrons 
nail have the same slate: 

The major shells arc specified by the principal quantum number, n. cor- 
responding to the distance ol the electron from the nucleus. 

^ Shells, in turn, have subshells based on the orbital angular momentum 
quantum number, ( 

t*" Hu turn, each subshell has subshells — called wbitais — which are based 
qn the z component of the angular momentum, nr 



So each shell rt has n - 1 subs hells, corresponding ] = (}, I, 2, n - 1 , .And In 

turn., each suhshell has tl + 1 orbitals, corresponding to m = -1, I, T. t. 



Much as witli (lie hydrogen alum, the various subshells (f - 1), 1, 2„ 3. 4, and 
so on) arc called the s, p. d. f, g. h. and so on states. So, for example. for a 
given /i. an s state lias one orbital (m = 0), a p stale has three orbitals {m = -1, 
0. and !), a d slate has five orbitals (m *-2. -1, 0. 1, and 2), and so on. 




In addition, due to the z component of the spin, m„, each orbital can contain 
two electrons — one with spin up, and one with spin down. 
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Perturbation theory Lets you handle situations like this — as Long as the pertuh 
batten isn't too strong. Ln other words, if you apply a weak magnetic field to your 
known system, the energy levels will be mostly unchanged but with a correc- 
tion. (Nate That’s why if* called perturbation theory and not dmsliNiil&fetmce 
theory ) The change you make to the setup is slight enough so that you can cal- 
culate the resulting energy levels and wave functions as <onerl/ons to the funda- 
mental energy levels arid wave lunclions ol Uie unperturbed system. 

So what does it mean lo talk oi perturbations in physics terms? Say that you 
have this Hamiltonian; 

11= LI, + AW | A -.4 1 

Here, H M is a known Hamiltonian, with known eigenfunctions and eigenvalues, 
and AW is the so-called perturbation Hamiltonian, where A- 4 indicates that 
the perturbation Hamiltonian is small. 

Finding the eigenstates of the Hamilton tan in this equation is what solving 
problems hke this is all about ■ — In other words, here's the problem you want 
to solve: 

H| v .)=(H 1 *AW]|^J=E n |^J (A-.1) 

The way you solve this equation depends on whether the enact, known solu- 
tions of H u are de^memie (that is, several states have the same energy) or 
nondegenerate. The next section solves the nondegenerate case, 

Working With Perturbations to 
Nondegenerate Hamiltonians 

Start with the case In which the unperturbed Hamiltonian. ir u . has na»deg/sn- 
emte solutions. That Is. for every state !$ T >, there’s exactly one energy, E. . 
that isn’t the same as the energy for any other slate: H |(S„) ■ E 4*1.) (Just as a 
one-to-one function lias only one x value lor any y). You refer to these nonde- 
generate energy levels of the unperturbed Hamiltonian as fr" 1 ,, to distinguish 
them from the corrections that the perturbation introduces, so the equation 
becomes 

From hern on. L reler to itm energy levels of the perturbed system as t . 

Tin: idea be hind perturbation theory is that you can perform expansions 
based on the parameter a (which is much, much less than 1) to find tlm wave 
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So your Task is to cal ruble E :r ' n and E'-\. as well as i|f l: n and So how dn 
yon do that in general? Time lo start s Liaising somi: main. Yuli start with three 
perturbed equal ions: 

I ** Hrtlnil Ionian: H |^} = ( H „ h AW )j ^ „) a E y,) ( A--* l ) 
y* Energy levels; E„ = E l '' l „ + AE ‘ 1 . + }.' E + | X, « I ) 

v* Wove function*: ) =|^„ ) + At| y ' 1 . ) + .A v | ly : . Ji + . . . ( X « 1 ] 

Combine these three equations to gel this Jumbo equation: 

( H, + AW }( |tf, > + -l 1 iy :L -' r } + A •' | tff ) 

= (e\ - At " t + x ! e:\ + .. .){ + X'\v'\) + ■■) (A « 1) 

Matching the coefficients 
of X and simplifying 

You tan. handle the jumbo equation in the preceding section by setting the 
coefficients of a on either side of the equal sign equal to each other. 

equating the revoth order terms its A on either side of this equation, here’s 
what you get: 

Ji.,|p,)=e'MpJ 



Now for the first-order terms in l: equating them on either side ol the Jumbo 
equation gives you 



H.|^}*Wfc}-E«|^.) + E«>,) 



Now e-epante the coefficients of }J in the jumbo equation, giving ynu 
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Okay, that’s the equation you derive from the second order In X. Now you 
have to solve lor E cll n . E ,;| n . and so oil using the zeroih-order, litst-order, and 
second-order equations. 
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Sfari by col i 1 1 l| that the unperturbed wavs Eduction, Id. > Isn't going to he 
very different From the perturbed wave Junction, I \f >. because the perturba- 
tioii is small. Tliat means that • :v ^ ] . in Fact, you can normalize I ijr^ so 
that is exactly equal In It 

(A,K}= l 

Given that y,,} = 0, J+ A ,J-A' '¥'*,)+• the equation becomes 
And because the cnuEEidents of A must both vanish, you get the following: 
This equation is useful lor simplifying the math, 



Finding the first-order corrections 

After matching the coefficients of *- and, simplifying {see the preceding sec- 
tion), yow want to find the first-order corrections to the energy levels and 
the wave functions. Find Uic first-order correction to the energy, E 1 '^ hy 
multiplying H,.||r"-,} + Wp r ) = E w f " ,|*„) hy ^,.1: 

W H ( ^ ' ,) + .;>,! w |^} = (*,| F.W, \ V ", ) + (#„| E ' , ]*, ) 

You can use H y_.) -(H,. * aW )| yr . ) 11, v ,,) ( is - * I ) to s i mp lily this to 

E m n =WwK) 

Swell, that’s the expression you use for the first-order correction, E clJ flr 

Now look Mo finding the first-order correction to the wave function. I v r '^ 
You can multiply the wive-funclioji equation by this nest expression, which 
is equal to ]j 

Xk>-I 

Ll 
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Finding the second-order corrections 



Now find the second-order corrections to the energy levels and the wave 
functions (the preceding section covers first-order corrections). To find F; : - : n . 
multiply both aides of H,. o ' h Ji+W y = E"". 
by I : 






3 I) 



0.,) 



(p..|ll» \y ( \ )-K^„|W \^\ ) • /*.| E m . | rtf 



j 1 * • \+U 
■ • + • ft. 



e«K\v' 






ir. 



* : 



This looks Like a lough equation until you realize tlial -= ij' „ I V i:, ,> ls oqual to 
z Liro. so you get 




And because -d',. ,» is also equal to wo. you get 

WwK*).{p.|if'.|p.) 



E' : - ; n is just a number, so you have 



And of course, because I (>„ > = 1 , you have 



Wole that if i y " M > is an eigenstate of W. list second-order correction equals 



zero. 

Okay, so E‘ M «■ a:o n i w Flow can you make that simpler? Well, 
(*J W ft. 

E'".-E 

R n = ( p -i ’ v ' , } gives you 



From using V> " L ) = X 






. Substituting that equation into 




w K"-)=W w IS: 



!*., W ft 






Bt 1 ! _ p 
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hrrr|W|*fll 

Now you have £".={$, w|d.) and F- 5 ,=£ ' - — z^j- 

™ h - w 1 

t.Ti-e total energy with the first- and second •order corrections:: 

£,=E''\ + AE\ + V&-'\+... (i«l) 



. Hires 



So from this equation, you can say 

I:'* | W I* ) 

E.=E„ + ^*,|W|i.).i : Xy- kf-*- 

■ltfT !_■ _ L- _ 






That gives you the first- and second-order corrections to the energy', accord- 
ing to perturbation theory 1 ’, 

Note that tor this equation to converge, the term in the summation must 
be small. And note in particular what happens to the expansion Term it the 
energy love Is arc degenerate: 

|{p,|wk)f 



III that case, you're going to end up with an E 1 *-^ that equals an E'" 1 ,,,. which 
means that the energy-corrections equation blows up. and this approach 
to perturbation theory is no good — which is to say Hurt you need a dif- 
ferent approach to perturbaiioa theory (coming up later in "Working with 
Pertubations to Degenerate Hamiltonians") to handle systems with degener- 
ate energy states. 

In the next section,, I show you an example to make the idea ot perturbing 
nondegen grate Hamiltonians more real. 

Perturbation Theory to the Test: 

Harmonic Oscillators in Electric Fields 

Consider the case in which you have a small particle oscillating in a harmonic 
potential, back and forth, as Figure 1 1-1 shows. 
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Figure 11-1; 
A la'npn: 
oscillator. 




Here h s the Hamiltonian for that particle, where the particle’s mass is m, its 
location isx and the angular frequency of the motion la n>, 



H = 



-y d 2 , 1 

2rri ilx- 2 






Now assume that (he particle is Charged with charge q, and that you apply A 
weak electric field, t. as Figure L 1-2 shows.. 






Figure 11-2; 

Applying 
an elHclric 
field to a 
harmonic 
oecillgtpr. 




> 



Z 






Tilt iurce due to the electric held Id this case is the perturbation, and the 
Hamilton ian becomes 



JI "fm^ 4 i mC?jV+ ^ 



In thi» section, you find the energy and wave functions of the perturbed 
system and compare them to the evact solutions. 
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That makes figuring out the second-order energy a tittle easier. 

Also* the following expressions turn out to hold for a harmonic oscillator: 



(o-l r v|rr}= n 13 - 



fi'* 



| 'Anuti i 



„ (ujIjt n\ 

With these four equations, you re ready to tackle q e 2, ■*. +..., the 



! p |B? _ cPl 



second-order correction to the energy, Omitting higher-power terms* the 
summation in this equation becomes 



< 4 c e 



E ,n, _ E |n;i 



fl'ff — 



E - H ;| 



And substituting in the fur K . -L'- n: n . and P'\ - gives, ynn 



l 


b+i \ x n >r 


<s € ~ 


- tli Vl 


fj : r 




tm 



Now, substituting in for -n t 1 Lv!n> and <n - 1 l.vl.n-- gives you 

,v , l^ 1 )* , * 

1 [-&«)[ 2 rma) 



nr 






I**-' 



( fT+1 ) , 

2 ht-r 






2m at 1 
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Note also that* E .. i 3 tu ' J and u’.- Ei fun 

These lour equations mean that 

k.) H fl ) * — ^w( rr 1 rr 0“ (n * L )} 



Nolo what this equation moans: Adding an electric field to a quantum har- 
monic oscillator spreads (Ere wave function of the harmonic oscillator. 



Originally, the harmonic oscillator’s wave function kg Just the standard har- 
monic oscillator wave function, lip > = ln>. Applying an electric field spreads 
the wave lunction. adding a component ot In - l >. which is proportional to 
the electric held, c. and the charge of the oscillator, if, like this: 



w.: = n\+-r—- 
■' .'i . -i - 



Q£ 



| 'Inna ) 



■ i 



kl' 3 - 1 }--) 



And the wave function also spreads to the other adjacent state. < n * 
like this: 



W-l") • I" - 1)- 1- 1 1 )> + 1)) 




You end up mixing states. That biending between states means that the per- 
turbation you apply must be small with respect to (fie separation between 
unperturbed energy states, or you risk blurring the whole system to the point 
that you can’t make any predictions about what's going to happen. 



In any case, that's a nice result — h lending the states in proportion to the 
strength of the electric lield you apply — mill it's typical of the result you gut 
with perturbation theory. 



Okay, Ural's how non degenerate perturbation theory works. As you can see. 
it’s strongly dependent on having the energy stales separate so itial your 



J py 
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Now multiplying that equation by #_ 1 gives you 

1 n 

x! f I E ■ , ) + { K I H , I#., ) | t. - e . X ". : K V, } 

Using | ho (act that l^> = l and «• 0 if m is not equal to n gives you 

a ,E,=.,E",H&.(j., H, i. ) 



Physicists often write that equation as 

{p=i,2X- j} 



where H-M^. >. And people also write that equation as 

i^ s -^=0 (£= 1,2,3,. 



where E' - E, - V"\. That's a system of linear equations, and the solution 
ex ists on ly iii 1 hen the dele? nninant to thtg array is nonvanishins; 



I'*,-*" 1 . IU I'... 

H„. 11.., . H,„ 



■■V 






H 






H 



fi/z 



H 



^3 






The determinant of this array is an ft h degree equation in Ef , and it has f 
dLFFereur roots. £0' . Those f dlFFerent roots arc the first-order corrections 

to rtiK Hamiltonian. Usually, Those roots are different because oE the applied 
perturbation. In other words, (Esc perturbation typically gets riel oJ the 
degeneracy. 

So here's the way you find the eigenvalues to the first order — you set up an 
fby-f matrls of the perturbation Hamiltonian, where H - <* I H lo >: 
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If,.- H,„ H,„ 

H^ e . 1 . H,,, 



H ,n, H,, H,„ ... H,,„ 



Then diagonalize this matrbc and determine the f eigenvalues E 1 and the 
marching e Igen vectors: 









/?" Ii 2. 3, j../ 






U.' 



TJien you get the energy eigenvalues to tirst order this way: 
K -E" ,J „ l E'"„_ (u = [,2.3,...,./) 

And the eigenvectors are 

^ 1 



in the next section, you look at an example to clarify thi? idea- 

Testing Regenerate Perturbation Theorg: 
Hydrogen in Electric Fields 

In this section, you see whether degenerate pert urbation theory can handle 
the hydrogen atom, which has energy states degenerate in dilterent angular 
momentum quantum numbers, when you remove that degeneracy by apply- 
ing an electric field.- This setup is called rise Sla/± affect. 
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]I, B 


1 1 .jii 






H,, 




H.,11 


H,., 









Drang The math gives you this remarkably simple result: 



H ,u 


H,,. 


1J,„ 


H,.u 


0 


f) 


1 


0 


H,:, 




1',, 




-a™. ° 


0 




0 


El,, 


E E J| j: 


ee, b 


EE, U 


1 


(t 


□ 


b 




H -V 


H,,, 




fj 


0 


0 


0 



Diagonalizing litis matrix gives you these eigenvalues — the first-order cor- 
rections to the unperturbed energies: 



v* E" 1 , - -3iM ri 
^E'VO 
w £' " 3 - 3efa„ 
F' ", - 0 

4 



where L (l , is Hid First-order correction to the energy of the I i > eigenfunction, 
EF'. Is the first-order correction to Uic encigy of the 1 2» eigenfunction, and 
so on. Adding these corrections to I he unperturbed energy lor the a = 2 state 
gives you the final energy levels: 



i> E *“ 

p-, = + 3fi*3, 










where ft is the Rydberg constant. Note this result: The Stark effect removes 
the energy degeneracy in 1 2Q0> and 1210s- (the l b and 13* eigenfunctions), 
hut the degeneracy in i 21 L -- and 21 - U (the 2 > and 1 4 ■> eigenlunct ions) 
remains. 
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Chapter 12 

Wham-Blam! Scattering Theory 



h This Chapter 

Switching between lib and cenler-of-mass frames 
Su I vlriL; flip Sc hrodinger eriuatlcn 
^ Finding the wavi- Eynctirjn 

Putting the Born approximation to work 



1 # out National Science Foundation grant finally came through, and you 
jb built your new synchrotron — a particle accelerator. Electrons and 
anti-electrons accelerate at near the speed of light along a giant circular track 
enclosed In a vacuum chamber and collide, letting you probe tire structure of 
t he high-energy particles ymi create. Ymi're sitting at the console ni ynnr giant 
new experiment, watching the lights Hashing and the signals on the screens 
approvingly. Millions of watts of power course through, the thick cables, and 
the radiation monitors are beeping, indicating that things are working. Cacti. 

You're accelerating particles and smashing them against each other to observe 
how they scatter. But this i& slightly more complex than observing how poo l 
bails collide. Classically, you can predict the exact angle at which, colliding 
objects will bounce off each other II l he collision Is eftrsfir (that Is. momentum 
arid kinetic energy are both conserved). Quantum mechanically* however* you 
can only assign probabilities to the angles at which things scatter, 

Physicists lssh Large particle accelerators to discover more shout the .struc- 
ture of matter, and that study is central In modern physics. This chapter 
serves as an introduction to that held oS study. You get to take a took at par- 
tick! scattering on the subatomic level. 

Introducing Particle Scattering 
and Cross Sections 

Think of a scattering experiment in terms of port ides in and panictes out- 
Luukal Figure 12-1, for example. In the figure, particles arc being scut in a 
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Figure 12-1: 
Scattering 
fram a 
targeE. 




stream from the left and Interacting with a target: most: of them continue on 
unscattered, but some particles interact with the target and scatter. 




Those part teles that do scaLler do so at a particular aitgle in three dimensions — 
that is. you give the scattering angle as a solid angle, dti, which equals sins? At do. 
where o ami u are the spherical angles I Introduce in Chapter U, 

The number of particles scattered into a specific dSi per unit time is propor- 
tional to a very important quantity in scattering theory: the differential cross 
section. 

.iivf n. 6 ) 

The diffemuuil emits seditm is given by \ ~ , and it's a measure ol the 

fai 

number of particles per second scattered into dtl per incoming Flux. The inci- 
dent flux, .1 (also called the current density), is the numherol inoidenl panicles 

dcr(^,tf) 

per unit area per un i 1 1 hue. So 5 s 

uli 

j d$[M) 

cO " J dQ 



where #) is due tiuinber of particles at angles 4> and 6. 



The differential cross section 



do j \ 

i&l 



has the dimensions of area, so calling; 





Chapter M\ Wtoani-BlamE Scattering Theory 277 



It a cross section is appropriate. The cross section is sort of like the site of 
the hull's eye when you’re aiming to scatter Incident particles through a 
specific solid angle. 




The differential caws section is the cross section for scattering to a specific 
solid angle, The tola ! rmu section, o, is the cross section for scattering of any 
kind, through any angle. So if the differential cross section for scattering to a 
particular solid angle is like the bull's eve, the total cross seel ion corresponds 
to I he whole target. 



You can relate the total cross section to the differential cross section by inte- 
grating the lolloping: 



ff= J 










da | 0 hint 

<K l 



Translating between the Center-of'Mass 
and Lab Frames 

Now you can start gelling Into the details of scattering, beginning with a dis- 
cussion of the center-of-mass frame versus the lab frame. Experiments take 
place In the tat? frame, but you do scattering calculations in the center-of-mass 
fmme. so you have to know hnw to translate between the two frames. This 
section explains how the frames differ and shows you how to relate the scat- 
tering angles and cross sections when you change frames. 



Framing the scattering discussion 

Look at Figure 12-2 — that’s scattering in the lab frame. One parlicle, travel- 
ing at i’li;,...., is incident on another particle that’s at rest - 0) and hits it. 

After the collision, the first particle Is scattered at arsgle w._ traveling at 
and the other particle Is scattered at angle and velocity iv,,,. 
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if 



V 



Y 


l lab 




2 lab 

# 








tr? 






-■rt 



I 2 

Figure 12-2: 

Scattering 
in Lhs lab 
Iranie. 



A 




Now in (he center-of-mass frame, the center of mass is 5 tat ionary ami th? par- 
ticles head toward each other. After they collide, they head away from each 
other at angles ft and ji-O. 

You have to move back and forth between these two l names — the lab frame 
and the cenler-ot-mass frame — so you need to relate the velocities and 
angles (in a noitrelativlslic way). 



Relating the scattering angles 
between frames 

To relate the angles ft, and 6, you start by noting that you nan connect j. : ,„i, 
and i’|.. lining the velocity ol the center of mass, r. 1 . L ,., tins way: 

^Lli* “ ®t> ■+ p inl 

In addition, here's what can say about the velocity of particle 1 after 51 col- 
lides with particle 2: 

" ' i :■■.'■ " '' lc * p .vi| 

Now you can find the components of these velocities;: 

I t-" f'lft* cose, = if,, cosft - if, 
v* p 'i, u p sinE'i =v' h sirs0 

Dividing the equation in the second bullet, by the one In the first gives you 
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Note also that lan(0 v ) = cot(y.). or tan^ft,} = tais( s /4- 

You know that ft : = */». anti tan{6 ; ) = tau(V* tells you that the following is 
trim: 



- V- - V: 



So substituting o, - */i into the preceding equation gives you. 
ft, - •/, - n. 
ft. + ft L = V- 



TJmrdurtL ft, and 9,. the angles of the particles iti the lab irame alter the colli- 
sion, add up in V- — which means ft. and 8, are at riy:lit angles with respect to 
each other. Cool. 



In this case, von can use the relations you've already derived to get these 
relations in tine special case where m, ■- m :j ; 



Arloi'l 


= 4cos(0, 1 

U 


dfl, 


rftr(o.O) 

dQ. 


=4c »! ! ) 
u 


dcr(^.0l 


= 4cos(^)- 

u 






=4 m] 



(K1 









Tracking the Scattering Amplitude 
Of Spinless Particles 

In the earlier section "Translating between the Cenler-of-Mass and Lab 
Frames,” you see how to translate from the lab frame to the center ohmass 
frame and back again, and those translations work classically as well as in 
quantum physics (as long as the speeds involved are nonrelativist Ic). Mow 
you look «t the elastic scattering of two spintess non relativistic particles 
From [he time-independent qmintoni physics point of view. 



Assume that the interaction between the particles depends cmlynn their 
relative distance, ! r, — >% I . You can reduce problems oE this kind to two 



J py 
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decoupled problems (see Chapter S for details). The first decoupled equation 
treats the center of mass of the two particles as a free particle, and the 

second cqnatimi is. lor a liclitious imrl iclcs ol mass 

1 m, r- m . 

Thelirst decoupled eq u.jikm, tire free-particle equaiiO]i oi the center of inass, 
is ol no in teres' to you ic scattering discussions. The second equation is the 

one to concern rate on, where u = m _ m : 






¥>!>■)+ V|>-j^|>-)=h>i>) 



You can use the preceding equation to solve for the probability that a par- 
ticle Is scattered into a solid angle ail — and you give this probability by l be 

different iat cross section, ^ „ 

ml 




In quantum physics, wave packets represent particles. In terms of scattering, 
(hn.se wave packets must he wide enough so that the spreading that occurs 
during the scattering process is negligible (however, Hie wave packet can’t be 
Sir spread that it encompasses the whult kb, including the particle detectors). 
Here's the crux: After the scattering, the wave function breaks up into two 
parts — an unscaltered part and a scattered pari. That's how scattering works 
in the quantum physics world. 



The incident Watfe fmetian 

Assume that the scattering potential ¥(r) has a very finite range, a Outside 
that range, the wave functions involved act like free particles- So the incident 
particle's wave Function, outside the limit of V[»> — that is, outside the range 
a from the other particle — Ls given by ihis equation, because V(r) is zero: 

r)+ M JV H> 

where A. = -^r 1 - - 

The form ' [r] ■+ , = n ls the equation for a plane wave, so <f M¥ .(c) is 

$ MV .(r) - Ae*' r , where A Ls a normalization Factor and k . ■ r is the dot product 
between the incident waves wave vector and r. in other words, you're treat- 
ing the incident particle ns a particle oF muhtentviiL A 1 ? ■ fik. 
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The scattered wade function 

Alter the scattering of the spinless particles, the nonscaltered wave function 
Isn't of much Interest to you, but the scattered wave function Is. Although the 
Incident wave function has the Form = Af' k ' ;r , the scattered wave func- 
tion has a R-ll^htly different form: 



The 1(4', 0) part is called the Sl vile ring innpiUixic?, and your job Is- to find ft. 
Here, A Ls a, normaifiiatioii [actor and 




where E is the energy ot the scattered particle. 



Relating the scattering amplitude 
and differential cross section 

The scattering amplitude of 5 pin Less particles turns out to he crucial to 
understanding scattering Erum l he tjuantuiu physics point uJ view. Tu see 
that, take a look at the (lux densities, (the llux density of the incident par- 
ade) and J (tlie Huh density lor the scattered particle): 

Inserting your expressions (or^,, and o„ into these equations gives you the 
following, where /('K d) is the scattering amplitude: 




Now in terms of the flux, the number ol particles dNOh n) scattered into <&i 
and passing through an area <#A = rrfti L$ 

dNfa9) = i^itn 
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And the homogeneous solution Is a plane wave — that Is, it corresponds to 
the Incident plane wave:: 



A:- v 



To lake a look at the scattering that happens, you have to (ind the particular 
solution, You can do that in terms of Gw n 's functions, so the solution to 
fv 1 eft I )^(rJ = ^+V(r) V r(r) 



is 



*¥ f<Kr-r’)V(/)»{r')rfV 






This intep/ttf breaks down to 



G{ r T ) ■ — r” 1 J <J 7— T dtt 

An ~i r-r 1 q -k 



. V ’ rl 



You can solve the preceding equation in terms of incoming and/or outgoing 
waves. Because the scattered particle is an outgoing wave, the Green's func- 
tion lakes this form: 





You already know that 

*?{ r ]| = Ac* ■+ f G( r - r' ] V{ r‘ ) V ( r \i V 



So subs t inning ^ ( r * r ) - 







.|„ p . r '\ into the preceding etpialLon gives you 

^y\mryy 



Wow, that's an integral equation tor qr(r) , ihc wave equation — tiuw do you 
go a Lout solving this whopper? Why, you use 8 tic Born approximation. «i 
course, 



-W 



riontf 
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The Bom Approximation: Rescuing 
the Wave Equation 

Okav, vrnir dilemma is to solve the Jul lowing equation for y(r), where = 

'"t r ) - 2 «F 1 v M v 



Ymi can rlfl- That with 3 series rtf successive approximations, called The Horn 
approximation (this is a laments result). To start, the zeroth order Bom 
approximation is just ^ (r) = fc lc (/). And substituting this zeroih-drikr Term, 
4T„(r). into the first equation in this section gives you the first-order term; 

>«,{r) = *„ - j&J |f ~j V(r, !*<r, M’r, 



which, tisinni|fg(r} = o M1 , <r) .gives you 



You get the second-order term by substituting this equation into 

^ l» “ P- = 2^7 J r^Ti V ( r ' M r ' ) d v ' : 



A r-rj 



M') = ^ - 2^/ 7I^i V ( r M' r '- 



And snhstiititing^,(r) = p„ ~ 9 ^ J J | tf _ v ( r V-i( f : ) d r '- in to the preceding 
equation skives you jT l r P| 



vAr) 



.i,p* 
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Tlit- pattern continues for the higher terms, which you can tied by plugging 
lower-order terns into higher ones. 



Exploring the far limits 
of the ufaOe function 



Now that you Vu used the Bom approximation (sue tike preceding section), 
take a look at the case where ris Large — in scattering experiments, r » r, 
where r& the distance from the target to the detector anti r is the size ol 

the detector. What happens to ^ ■ r I ■ 4i| « V{r')tfr(r r )tTr , the 

exact integral equation tor the wave function, when r r? Here's tire answer: 






U f e* - 
2xh : f\r r| 



Vfr'Mr'JdV 



Because r » f, yon can say that frlr-^l ® fer- ft • r. where b - r 1 is the dot 
product of it and r (k is The wave vector of the scattered particle). Anri 




Using the last two equations in 
gives you 



f—jt V [r }v [ r- ] d V 



yf(r )= Ae *' ' + ^-f{d>,0) 



r — ? ™ 



And htre 

'<**>■ 5§M* - V(/)v{r')<fV = ^p-:>|V|v-) 

dc | $ ,ff S | . p 

The differential cross section is given by — — ■ ![$.#} \ , which in this 



case becomes 
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Using the first Bom approximation 

If the potential Is weak, the Incident plane wave Is only a little distorted and 
the scattered wave Is also a plane wave, That’s the assumption behind the 
first Horn approximation, which you take a look al here, So If vou make the 
assumption that the potential is weak, you can determine from the equation 

V. l>> = f j~“j V ( h ) V, { r, ) tf'r, that 






u : e+ r 
2xty f r-r r \ 



vU'K 



Okay, so what is fl u. ';■)? Well 

•"Vir'Mr’Kr’ 



And this equals fhe follow ins', where <i = it,,, -k: 



'M-i&J* ■‘■WM’W-zfci* ^KP'- 



dui 6>,8 | | . : 

And because ^ , you have 



n* 

d£i -hr ti ' 




When the scalterinp Ls elastic. The magnitude of ft is equal to the magnitude 
of and you have 

ej = I Jt,j - h I = 2ft £in('j'0 

where 9 is the angle between it* and ft. 

In addition, if you say that V(j-) is spherically symmetric, and you can choose 
the a axis along tv, then q -t - w cosfl'. so 

'{*'*)= ^Orl* * r 'V [/)</>' = ; ^Ljr'V(r sinfl'dW'l^' 

That equals 

u JTirl £ 5 I Sr* ii 
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Because ' 



— - 1 *{*■& )| . you know that — ^ ^ ■ \ f ( «V »)| ’ , so 



drs\ p.M J -i jj 



Al 



<fh' 



J /V( r ' ) sin ( qr ' jifr' 



You've come far jn this chapter — from the SchriJdinger equation all the 
way through the Bom approximation, and now to the preceding equation lot- 
weak. spherically symmetric potentials. How about you pul this to work with 
5 i>me con crctc mi mhers';' 



Putting the Bom approximation to Work 

Iri This section, you find the differentia! cross sect ion fr»r two electrically 
charged particles of charge Z t e and '£&. Here, the potential looks like this: 

s/irt ^ 



So here's what the dil IcmdiaL cross section looks like in the Firs! Born 
approximation: 

4Z, J Z/i>y 

q>!\' 

And because J sir ' itr q . you know that 




j 

Ai 



4Z, : Z re* it : 
Ai ~ tfV 



And because <} - 2 fcs In {"/:), the tallowing Is true: 
j AX , % V V * Z,'Z,V , ( # ■] 

" ?v ” m 2 sn U J 



where E is die kinetic energy of the incoming particle: E - 



fyk 1 

- 



Now get more specific; say that you're smashing an alpha particle. Z, p 4, 
against a gold nucleus. Z : « 79. If the scattering angle In the lab frame Is 60 . 
what Is H in the center-ol-mass Frame? 
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The Part of Tens 



The 5 th Wave By Rich Tennant 









qxf 1 ^ m * V-*v- vj *,.st 

(<r 2 iT*i / ^ ^ *«IT> 

' * jt4 * n ' : i ^ e„ If- *■ If" Tfflifr*'" 

V?~ %V ■fr + #-- /l % 5 15 *tnrS 3*3 
4 [7---^-^ - y SoSsj'+!!jJ^i nft'TiHfl 
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^ 44 ^ 1^ - J * - J 
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Chapter 13 

Ten Quantum Physics Tutorials 



h This Chapter 

Understanding basic concepts and equations 
Viewing iJlustrarinns jmrS animations 



MMfivm scientists start mixing talk oE dice, billiard balls, and a possibly 
TT Lindetatl cat-in-a-box. you know you're dealing with a challenging sub- 
ject. Luckily, yon can firui pLnnly of online tutorials, some ol Ihmn featuring 
attiitialLmi, uj help you wrap your brain around quantum physics. This chap- 
ter presents a good starter list. 



An Introduction to Quantum Mechanics 

www . cb a:ni s ti y .aMu-atat*. ■.■■.- /bet ha / iff! 

What is a itme function? What is an orbital?: An fntmcluctitm to Quantum 
Mechanics comes Iron] Neal McDonald. Midori Kitagnwa-JVI.r-on. Anna 
Timasheva, Heath Hanlin, Zlt Lilas. anti SherwLn J, Singer at The Ohio Stale 
University. This site Includes tutorials on probability, particles versus waves, 
wave Junctions, and more.. Including Shockwave-based sound ('though If you 
don't have Shockwave installed, that's not a problem). 



Quantum Mechanics Tutorial 

www . q L 1 t£t v . com / tu tor i a 1 s / quan t urn . h tml 

This cool tutorial is one of the Flash-animated Modern Physics Tutorials by 
Giles Hogben. Extensively illustrated, this tutorial probes questions such 
as wave-panicle duality and offers a good general introduction to quantum 
physics. 
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Grains of Mystique: Quantum 
Physics for the Layman 

www . f»qa . org/doca/qp 

This site provides good historical and experimental background info— and 
they've documented their sources and made some attempts at .peer review. 



Quantum Physics On tine Version 2.0 

www . cf jarir cm - phys ic s . po ly t tebn ique . f e f ir.d^x . ht m 1 

Tli is is a cool set of programs that run in your browser, giving simulations 
of various quantum physics experiments, it's by Manuel joffre. Jean-Louis 
Basdevant, and Jean Dalibard of the Eeole Polyiechnique in France. Look for 
Information on wave mechanics, quantisation, quantum superposition. and 

spin 



T odd K. Timhertake's Tutorial 

f a col t yweb . be rry . edu 1 1 inter lake/ qeha oa /qp . Lit m 1 

This tutorial is by Todd ft. Tlmberlake, assistant professor of the Department 
of Physics. Astronomy. & Geology of Berry College In Georgia. It’s a fairly 
brief but well-written Introduction to the Ideas of quantum mechanics. 



Physics 2b/7’s Tutorial 

www . pbysics247 . com/physi cs tutorial/ quantum physics- 
bil 1 iards . shtml 

This is a text- based tutorial from Physics 24/7. ft Includes material on quanta, 
the uncertainty principle, and quantum tunneling (as well as some ads). 
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Start ZochoWski’s PDF Tutorials 

www . e-mccp , lit L . <a c . tik / - swz j c our se a / SBO 5 5 .■' SU *55.h tm! 

Sian Zocho wiki, from the department of Physics &. Astronomy at University 
College London, pul together these PDF-based tutorials on quantum phys- 
ics. These are tutorial handouts for a Quantum Mechanics course at the 
University College, and they serve as an excellent Introduction to quantum 
physics. 



Quantum Atom Tutorial 

ww . co lorado . e do / physics / 2 000 / quantunizanfi f index . hi m 1 



This is a fun. cartoon-centric tutor ial on the quantum nature of the atom from 
the University of Colorado Physics 2000 project. 



College of St. Benedict’s Tutorial 

wvw . phys ics. c sbs j u . edu / QX ' i ndex .lit ml 

This la a comprehensive quantum physics tutorial from the College of St. 
Benedict. It's a good, more serious, lest and equal lons-based tutorial with 
plenty of Illustrations. 



A Web-Based Quantum 
Mechanics Course 

elec trend . phys . u t V , edu ! cirri / Modul ee . h tm 



This one’s from the University of Tennessee, and it's on extensive online 
course In quantum physics. It Includes modules on square potentials, har- 
monic oscillators, angular momentum, spin, and so on, 
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Chapter 14 

Ten Quantum Physics Triumphs 



In This Chapter 

Explaining unexpected results 
Identifying characteristics of the quantum world 
^ Developing new models- 



f Quantum physks has been very successful In explaining ninny physical 
phenomena, such a$ w.ivc-parlicle duality. In fact, quantum physics 
was created to explain physical measurements that classical physics couldn't 
explain. Thi$ chapter is about ten triumphs of quantum physics, and it points 
you lo resources on the Web that examine those triumphs lor further 
information. 



WaPe~ Particle Duality 

Is I hat particle a wave? Or is that wave a particle? That's one of the questions 
that quantum physics was created to solve, because particles exhibited wave- 
like properties in the lab, whereas waves exhibited particle-like properties. 

These Web sites offer more insight: 

I s*" vmw . gilestv. com / fa tor i a 1 s / quail t uni . h t ml 

f" www . phys ' ■ ■ : , com/phys ics- fa tori a 1 /quart turn- phys acs- 
b i 1 1 i ^ rdrs . -sh hm 1 



The Photoelectric Effect 

Another founding pillar of quantum physics was explaining the pftoiQektcuic 
effect, in which experimenters shone light on a metal. No matter how strong 
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the light, the energy of ejected electrons from the metal didn’t rise. U turns 
out that the energy of electrons goes up with the frequency of the light, not 
Its Intensity — which gives support to the light as a stream ol discrete pho- 
tons theory. 

For more info on the photoelectric effect, check out ww.g i lestv . com/ 
i li tor lals / qxtar. fjcr. .html. 



Postulating Spin 

The Stern-Gerlach experiment results couldn't be explained without pos- 
tulating spin, another triumph of quantum physics- This experiment sent 
electrons through a magnetic field, and the classical prediction is that the 
electron stream would create one spot of electrons on a screen — hut there 
were two (corresponding to the two spins, up and down). 

TlliS Web Site has more iillo: electrons . pliys . Jtk . edjycanl /cr.ednles/ 
in , h t rr.. 



Differences between Newton's 
Laws and Quantum Physics 

In classical physics, hound particles nan have any energy or speed, hut that's 
not true in qiisnSum physics. And in classical physics, you can determine 
both the position and. mo men turn of particles exactly, which isn'L true in 
quantum physics (thanks to the I leisenberg uncertainly principle). And in 
quantum physics, you can superimpose states on each other, and have par- 
ticles tunnel into areas that would be classically impossible. 

You can find a nice discussion of the differences between classical and 
quantum physics at f aeul tyi^b - berry -edu/t timber lake /gthao*/ 

. hLnl r 



Heisenberg Uncertainty Principle 

One off the triumphs of quantum physics is the Heisenberg uncertainty princi- 
ple; Heisenberg theorized that you can't simultaneously measure a particle's 
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Square Welts 

Like trarmoinc Oscillators, quantizing pirtides bound in square wells at the 
micro level was another triumph lur quantum physics. Classically, particles 
in square wells can have any energy, but quantum physics says you can only 
have certain allowed energies. 

There's, plenty on the Web about U, including these two good treatments: 

I ts* ww . physi e* r c sb* j u - *clu /QH/ Indese , httal 
v* el ectrone -johy * - ntk - edu / qmi /nodul e* /ntadnlei - h Lm 



Schwdinqer's Cat 

SchrGdinger's Cat is a thought experiment that details some problems that 
arise 1 in the macro world from thinking; of the spin of electrons as completely 
tion-determined until you measure them. For example, If you know the spin of 
one of a pair of newly-created electrons, you know the other has to have the 
opposite spin. , So if you separate two electrons hy tight years and then mea- 
sure the spin of one electron, does the other electron's spin suddenly snap 
to rhe opposite value — even at a distance that would take a signal Emm the 
first electron years lo cover? Tricky stulil 



For more, lake a look al www r sri I e S t v . com/ tutorials/ cpm nt cur. . hi m L . 




Glossary 



ere's a glossary of common quantum physics terms; 



amplitude: The maximum amount ol displacement uFan oscillating particle. 

angular momentum: The product of the distance a particle is Irum a Certain 
point and its momentum measured with respect to the point. 

Mini!) Nation operator: An operator that lowers the energy level of an eigen' 
state hy one level, 

anii'Kannitian: The value you get when you take the Hermltlan adjoint of an 
expression and get the same thing back with a negative sign in front of it. 

black body: A body that absorbs all radiation and radiates it all away. 

flehr radius: The average radius of an electron 's orbit in a hydrogen atom, 
about L0 _|: meters. 

bound Atate: A state in which a particle isn't free to travel to Infinity. 

hoftMUK Particles with integer spins, including photons, pi mesons, and so on, 

bru-kel notation: Abbreviating the matrix Form of a state vector as a kef. 
or I (jr-, and abbreviating the kefs complex conjugate, or fern, as «qr . 

cenlcr-af-mass frame: In scattering theory, the frajue in which the center til 
mass is stationary and the particles head toward each other and collide. See 
also lab I ramie . 

central potential; A spherically symmetrical potential. 

commute: Two operators commute with each other if their commutator is 
equal to zero. The wnumtjitotofot operators A and H Is [ A, Hj = AH - BA. 

complex conjugate: The number yon get hy negating the imaginary pari oi a 
Complex number. The * symbol indicates a complex conjugate. 
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Conipion effect: An Increase of wavelength, depending on the scattering 
angle, that occurs after Incident light hits, an electron at rest. 

conservation of energy: The Saw of physics that says the energy oF a closed 
system doesn't change unless external influences act on the system, 

creation operator; An operator that raises the energy level oS an eigenstate 
by one level. 

current density: See incident Ibis. 

electron volts (tV): The amount of energy one electron gains falling Hi rough 
i volt. 

diagonalize: Writing a matrix so that the only nonzero elements appear along 
the matrix's diagonal, 

differential cross section.: In scattering theory, the cross section For scatter- 
ing a particle to a specific solid angle; it's like a bull's-eye. 

Dirac's constant: Planck's constant {h - 6.626 x lo S1 Joule-seconds) divided 
by 2n„ It's represented by an h with a bar going through 11. 

Dirac ttvlnlbti; See bra-kel notation. 

eigenvalue: A complex constant that represents the change in magnitude of 
a vector, 

eigenvector: A vector that changes ha magnitude hut not direction after you 
apply an operator. 

rlwlic collision; A collision in which kinetic energy is conserved. 

electric field: The force on a positive test charge per Coulomb due lo other 
electrical charges. 

electron: A negatively charged particle with half-integer spin, 
emisslvfly: A property nt a substance showing how well it radiates, 
energy: The ability oi a system To do work. 

energy degeneracy; The number of slates that have the same energy, 
energy well: See potential well. 

expectation value: The most probable value an operator will return. 
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fermion*: Particles with half-integer spin, including electrons, protons, neu- 
trons, quarks, and so on. 

frequency: The number cd cycles of a periodic occurrence per second. 

Hamiltonian: An operator for the total energy of a particle, both kinetic and 
potential. 

Heisenberg uncertainty principle: Jtef uncertainly principle. 

Hermilinn odjcsinb A value, represented as A 1 , that replaces complex num- 
bers with iheir complex conjugates, swaps bras and kels, and replaces opera- 
tors with their Hermltlan operators., 

Hennethn operator Operators that are equal to their Hermltlan ad joints: in 
Other words, an operator 5s Her m Ilian IF A' = A. 

incident flux: The number of Incident particles per unit area per unit time. 

inelastic collision: A collision, in which kinetic energy isn't conserved. 

intensity (wave): The time-averaged rate of energy transmitted by a wave per 
unit of area. 

.Me: The MKS unit of energy — one Kew ton-meter, 
fcet; Se.fi bra-ket notation. 

kinelk energy: The energy of an object due to its motion. 

lab frame: in scattering theory, the frame in which one particle is incident on 
a particle at rest and hits il. See also centeieaf-mass Ira me. 

Lnpkdoh: Aji operator, represented by A. that you use to find the Hamiltonian. 

magnetic field: The lorce on a moving positive test charge, per Coulomb. 
From magnets or moving charges. 

magnitude: The size or length associated with a vector (vectors are made up 
of a direction and a magnitude-). 

mass: The property that makes matter resist being accelerated, 
momentum; The product ol mass times velocity, a vector. 



MKS system: The measurement system that uses meters, kilograms, and 
seconds. 
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Newten: The MKS unit of force — one kilogram-meter per second-. 

normal] rad function: A function In which the probability adds up to 1. 

arfoliata: Different angular momentum stales of an electron, represented ns 
subs hells in atomic structure. 

orthogonal: Two lefts. I qr* and b--. lor which :tfdo=- = f). 

orthonormal: Two kets, I qr> artel that meet the tallowing conditions: 
'-iy o> = 0; 41 1 i+f-- = 1 ; and Id-- = L 



oMlUne: To move or swing side to side regularly, 

pair annihilation: The conversion of an electron and positron Into pure light. 

pair production.: The conversion of a high-powered photon info an electron 
and positron. 

particle; A discrete piece 0! matter. 

Pauli exclusion principle: The idea that no two electrons can occupy the 
same state in a single atom. 

period: The fLme It takes lor one complete cycle of a repealing event, 

perturbation: A stimulus mild enough that you can calculate the resulting 
energy levels and wave Fisnctions as correct ions tn the fundamental energy 
Levels and wave functions of the unperturbed system. 

photoelectric effect: A result in which the kinetic energy of electrons emit Led 
from a piece of melal depends only on Lite frequency ■= not the intensity — 
of llbe Incident light. 

photon: A quantum of electromagnetic radiation. An elementary particle that 
Is Its owe antlpartlcle. 

pi meson: A suhatortiLr particle that helps hold the nuclei is nJ an atom together. 

Flanek’s constant: A universal constant, h, that describes the relationship 
between the energy and Frequency of a photon. It equals 6.62ti x 10 Joule- 
secuuds. 

positron; A positively charged anti-electron. 

potential hairier: A potential step of limited extent; an electron may be able 
to tunnel through the barrier and come out the other side. 
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potential energy; An object's energy because of its position when a force is 
acting on It or Its Internal configuration. 

potential step: A region in which the energy potential Forms a stair shape; 
a particle striking the step may he reflected or transmitted. 

potential well: A region in which, there's a dip in the energy potential thresh- 
old; panicles without enough energy to overcome the barrier can become 
trapped in the well, unable to convert I he potential energy to kinetic. 

power The rate of change in a system's energy. 

probability amplitude The square root off the probability that a particle will 
occupy a certain state 

probability density: The likelihood that a particle will occupy a particular 
position or have a particular momentum, 

quantized: Coming in discrete values. 

quark: Particles that combine with antiquaries lo form protons, neutrons, and 
s o on. 

radian; The MKS unit of angle — 2s radians are In a circle. 

radiation: A physical mechanism that transports heat and energy as electro- 
magnetic waves. 

scalar: A simple number (without a direction, which a vector has). 

Sell rodinger equation: A wave function that describes how energies and 
probable locations of electrons change over lime. 

simple hnrmuiilc molten; Repetitive motion where the restoring force is pro- 
portional to the displacement. 

spherical coordinates: Coordinates that indicate location using two angles 
and the Length oFa radius vector. 

spin: The intrinsic angular momentum ol an electron, class ilied as up or down. 

synchrotron: A type u! circular particle accelerator. 

state vector; A vector that gives the probability amplitude that particles will 
be In their various possible slates. 
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Ihrraliokl frequency: II you shine light below this frequency on metal, no 
electrons are emllted. 

total cKMHjwctbn: tri scattering theory,, the cross, section for any kind of par- 
ticle scattering, through any angle. 

tunneling: The phenomenon where particles tract yet through regions that 
they're classically Jorbidden to jjo. 

ultraviolet catastrophe: The failure ol the Raloigh-Jcans Law lo ex plain Ij lack- 
body radiation at high frequencies, 

uncertainty principle: A principle that says Its Impossible to know an 
object’s exact momentum and position, 

vector: A mathematical construct that has both a magnitude and a direction. 

velocity: The rate of change of an. object's position , expressed as a vector 
whoso magnitude is speed. 

volt: The MKS -unit of electrostatic potential — one Joule per Coulomb. 

wave; A traveling energy disturbance. 

wavelength: The distance between crests or troughs of a wave. 

wave-particle duality: The observation that light has properties of both 
waves and particles, rle pending on the experiment. 

wave packet A collection of wave fund ions such that the wave fund ion ■; 
.Interfere constructively at one location and interfere destructively (go Lo 
zero) at aid other locations. 

work; Force multiplied by the distance over which that force acts. 
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perturbations, 11511 
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free particles. See q(so bound states 
Sehroclinger equation, 85-89, 172-173 
spherical coordinates, 195-196 
In three dimensions,, 172-177 
time dependence, 175-177 
total energy equation. 174-175 
wave packet, 87-89 
jt, y ? and z equations, 173-174 
frequency, 305 

g state. 251 
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center of mass, 208—209 
in electric fields, 271-273 
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Identical noninteractlng particles, 247— 25i> 
Identity matrix, S2 
Identity operator, 33 
Incident current density, 71-72 
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normalizing wave Function. 64-65 
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The Hamiltonian 

One of the central problems ei quantum (Mthanies is to CBkulfite the enarsy lewis of & system The 
energy operator, called the Humffftcnlsr), H, gives yurt me total energy, Finding ihe energy fevefa of K 
system breaks (town io findingtha eigenvalues ef ihe problem 

H kMk) 

Here's the shim eg nation in matris terms: 



del 



H -E H 

1 1 it 



H 



I 



L-E I-L 



H, 

H, 



H, 

H. 



H„-E H, 



K. 



H., -E 



= 0 



H V[r ) ! = r ) + V {* ) W r )= 





The Heisenberg Uncertainty Principle 

The Heisenberg uncertainty principle says that Jhe better you know the position gl a particle, (he less you 
knew the momentum., and vice versa. In the i direction, to r example, chat looks like this: 

Aafyp, £ | 

where i* is the measurement uncertainty iin the particle's jv position, &p t is its measurement uncertainty in its 
momentum in the j direction, and h = 

This relation ho Ids far ell three dimensions: 



The Sehrodinger Equation 



The Schrbdinger equation describes the energias and p rot ah 'h Ic cations o' electrons. The quantum phys 
ics in this book is largely a bout solving this differential equation lore variety oF potentials, Vl/|: 

M'h^nn 
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